Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	43.3%	115.3 / 24.5	67	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (31/7);  new base C (0/29);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.08e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993764‑994063                                                                                                                                                         |                                                                                                                                                       cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGc                                                                                                                                                       <  1:953234/151‑1 (MQ=255)  aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGcc                                                                                                                                                      <  2:1012882/151‑1 (MQ=255)   cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGccc                                                                                                                                                     <  1:172740/151‑1 (MQ=255)    ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGccc                                                                                                                                                     >  1:92383/1‑150 (MQ=255)      tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCccacc                                                                                                                                                  >  1:412640/1‑151 (MQ=255)       aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCccacc                                                                                                                                                  >  2:576288/1‑150 (MQ=255)        gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                >  2:9849/1‑151 (MQ=255)                      tGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGa                                                                                                                                   >  1:465619/1‑150 (MQ=255)                       gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                 >  1:797071/1‑151 (MQ=255)                            atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                            >  1:339028/1‑151 (MQ=255)                                   gcgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                     >  2:1033411/1‑151 (MQ=255)                                    cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                    >  2:125659/1‑151 (MQ=255)                                    cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                    >  1:975869/1‑151 (MQ=255)                                       gtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                  >  1:503465/1‑150 (MQ=255)                                           tCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  2:535042/1‑150 (MQ=255)                                            ccTGAATTTCAGCCAGGGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCCCCC‑CCGCTGCGGCCCCCCACCGCCACAGACTGACCTGCCGGAATATGTAAATCCAg                                                                                                            <  2:799318/151‑1 (MQ=255)                                              tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          >  1:637385/1‑151 (MQ=255)                                                aaTTCCACCCGGGGGTGTGGTGCCCGCTAACACAGCGCTTGGGGTTGGCT‑TTTCCACCCCAGCAAGGGGGCCCACCGGGGTATTTTCG‑CCACCC‑CCGCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          <  2:302353/149‑1 (MQ=255)                                                 aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       >  1:1130113/1‑151 (MQ=255)                                                  tttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                      >  2:579556/1‑151 (MQ=255)                                                  ttccagccggcGGTGTGGTCCCCGCTAACACACCGCCGGGGGTTGGCT‑TTCCCAGCCCTGCCAGCAGGCGCAGC‑GGGTTCTTTTGCCCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                           <  2:428115/143‑1 (MQ=255)                                                      agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  >  1:1150554/1‑151 (MQ=255)                                                          agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                               >  2:593422/1‑150 (MQ=255)                                                             ggTGTGGTCCCCGTTAACACATCCCCTGCGTTGGGCT‑TTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCCCCC‑CCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                           <  2:92907/151‑1 (MQ=255)                                                               ggtggcccccccaaaCACCCCGCGTGCGGTTGGCT‑TTCCCAGCCCTGCCGGCGGCCGCAGAAGGGTATTTTGG‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  1:393914/138‑1 (MQ=255)                                                                gtggtgCCCGCAAACACACCCCCTGTGTTTGGCT‑TTCCCACCCCCCCCAGCGGGCGCAGCAGGTTTCTTTCC‑CCCCCC‑CCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                             <  1:271890/146‑1 (MQ=255)                                                                  ggtcCCCGCTAACCCACCGCCTGGGGTTGGCT‑TTTCCAGCCCTGCCAGCGGGCGCAGCGGGGTATTTTCC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  2:869372/148‑1 (MQ=255)                                                                  gggccccGCAAACCCACCGCTGGCGGTGGGCT‑TTCCCAGCCCGCCCCGCGGGCGGAGGGGGTTTTTTTTC‑CCCCCC‑CCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                      <  1:939940/147‑1 (MQ=255)                                                                         gCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                               >  2:394748/1‑151 (MQ=255)                                                                            accccaccGCCTGGGTTTGGCT‑TTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTCCTTTTC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            <  2:282288/144‑1 (MQ=255)                                                                            aaCACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            >  2:1148027/1‑151 (MQ=255)                                                                              caaatccccggggtttGGCT‑TTCCCAGCCCTCCCACCGGGCGCAGCGGGTTTTTTTTC‑CCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                           <  2:872484/136‑1 (MQ=255)                                                                               cctccGCCTGGGTTGGGCT‑TTCCCACCCCTGCCAGCGGGCGCAGCGGGGTCTTTTTC‑CCCCCA‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                         <  2:549917/147‑1 (MQ=255)                                                                               cacaccgcTGGGGTTTGCGTTTTCCACCCCGCCCAGCGGGCGCGGCGGGGTTTTTTTC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          <  1:1007683/144‑1 (MQ=255)                                                                                caTCGCCTGCGGGTGGCG‑TTTGCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTG‑CCACCA‑CCGCTTCGGCCCACCACCGTCACAAACTGACCTCCCGGACTCTGTAAATCCAGTTGGCTCAGGACGATATTTTCTGCGTAATGttt                                                                           >  1:1143523/1‑148 (MQ=255)                                                                                  tCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       >  1:700979/1‑150 (MQ=255)                                                                                        gCGGTGGGCT‑TTCCCAGCCCTGCCAGCGGGCGCAGCAGGTTTCTTTGC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                <  2:1061813/151‑1 (MQ=255)                                                                                         cGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                               >  1:347505/1‑151 (MQ=255)                                                                                            tgggCT‑TTCCCAGCCCTCCCAGCAGGCGCAGCAGGGTTCTTTGC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                 <  1:1032959/144‑1 (MQ=255)                                                                                             tgggctttcccaccccTCCCAGCGGGCGCAGGGGGTTACTTTTC‑CCCCCC‑CCGCTGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                             <  1:973802/136‑1 (MQ=255)                                                                                             tGGCT‑TTTCCAGCCCTGCCAGCAGCCGCAGCGGGGTCTTTTTG‑CCACCC‑CCGCTGCGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  2:685441/150‑1 (MQ=255)                                                                                               ggtttttCCAGCCCGGCCAGCGGCCGAAGCGGGTTACTTTCC‑CCCCCC‑CCCGTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  2:966618/145‑1 (MQ=255)                                                                                                cG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                        >  1:528780/1‑151 (MQ=255)                                                                                                   ttcccgccccTGCCGGCGGGCGCGGAAGGTTCTTTTTC‑CCACCC‑CGGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       <  1:390247/126‑1 (MQ=255)                                                                                                     tCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    >  1:429929/1‑151 (MQ=255)                                                                                                      ccagcccgcccAGCGGCCGCAGCAGGGTCTTTTCG‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:492491/142‑1 (MQ=255)                                                                                                          aCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               >  1:871985/1‑151 (MQ=255)                                                                                                           ccgcccAGCGGGCGAAGCGGGGTATTTTCC‑CCCCCC‑CCTCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               <  2:742374/146‑1 (MQ=255)                                                                                                            cTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  2:268051/1‑151 (MQ=255)                                                                                                                     ggCGCAGCGGGTTACTTTTG‑CCACCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  2:241935/149‑1 (MQ=255)                                                                                                                         cggagtgttttttttC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  1:464725/134‑1 (MQ=255)                                                                                                                         cagcagtttttttttG‑CCCCCC‑CCGCCGCGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  1:636496/138‑1 (MQ=255)                                                                                                                              gggtttttttC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAgaa                             <  2:105207/143‑3 (MQ=255)                                                                                                                              gggTTCTTTTG‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                             <  2:118681/149‑1 (MQ=255)                                                                                                                              gggTTCTTTTC‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  2:772851/146‑1 (MQ=255)                                                                                                                               ggTACTTTTG‑CCCCCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                          <  2:548707/151‑1 (MQ=255)                                                                                                                               ggTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  2:679765/1‑150 (MQ=255)                                                                                                                                gTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                         >  2:992414/1‑151 (MQ=255)                                                                                                                                    ttttG‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                        <  2:953408/148‑1 (MQ=255)                                                                                                                                    ttttC‑CCCCCC‑CCGTCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:647667/149‑1 (MQ=255)                                                                                                                                     tttgcccaccc‑CCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:1115882/144‑1 (MQ=255)                                                                                                                                     tttC‑CCCCCC‑CCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     <  2:359744/150‑1 (MQ=255)                                                                                                                                      tttccccccccccGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                     <  1:364205/124‑1 (MQ=255)                                                                                                                                           cccccgcgggggcgcCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:1042043/133‑1 (MQ=255)                                                                                                                                           ccccccccGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                  <  2:919417/143‑1 (MQ=255)                                                                                                                                              ca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc            >  2:353604/1‑151 (MQ=255)                                                                                                                                                  cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc         >  1:640528/1‑151 (MQ=255)                                                                                                                                                    cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt       >  1:1152738/1‑151 (MQ=255)                                                                                                                                                         gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  1:520348/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993764‑994063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑