Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,603,342	0	A	C	40.0%	88.3 / 56.1	65	T509P (ACC→CCC)	hyfB	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (36/3);  new base C (0/26);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 3.64e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCG  >  NC_000913/2603193‑2603491                                                                                                                                                         |                                                                                                                                                       cgcgCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑c                                                                                                                                                      <  2:233863/151‑1 (MQ=255)   cgcAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑ccc                                                                                                                                                    >  2:1196762/1‑151 (MQ=255)      aGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTc                                                                                                                                                 <  1:230201/151‑1 (MQ=255)          cacacGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGc                                                                                                                                             <  2:786355/151‑1 (MQ=255)             acGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCAc                                                                                                                                          >  1:329463/1‑151 (MQ=255)               gCTGTAGACGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACt                                                                                                                                         >  1:1195394/1‑150 (MQ=255)                  gAAGGGGCGCAGGAAGTGCAAAATACGATGATCGTCGCCATGCTACTGCGCGCGGCCCTCTGCTTATTAATTGCGCTTAGTGCC‑AGTTGGCTGCCACCGAAGATAATGCATTTTGCCCTTGCGTTTCCC‑AA‑TC‑CCCCTCCCGCCACTGTc                                                                                                                                      <  1:81010/150‑1 (MQ=255)                      agGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGcc                                                                                                                                   >  2:1073699/1‑149 (MQ=255)                           caGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGa                                                                                                                             >  1:673404/1‑150 (MQ=255)                                 gtccaaaaTACGATGGCCGCCGCCTTGCTCCTGCCCGCGCCACTCTGTTTATTATTTGCGCTTGGTGCC‑GGTGGGCGGCCACAAAAAAAATTCCATTTTCCCCATGCTTTTCCA‑AA‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCAc                                                                                                                      <  2:1334238/146‑1 (MQ=255)                                        aTACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAcc                                                                                                               >  2:804485/1‑151 (MQ=255)                                          aCGATGTTCGTCGCCATGCTCCTGCCCGCGGCCCTCTTCGTATTAATTGCGTTTGGTGCC‑TGTTGGCTGGCACCGAAAAAAATGCATTTTGCCCATGCTTTTACA‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAccc                                                                                                              <  2:105901/150‑1 (MQ=255)                                            gatgatCGTCGCCATGCTCCTGCTCGCGCCCCTCTGCTTATTATTTGCGTTTGGTGCG‑TGTTGGCGGGCACAGAAAATAATGCATATGCCCCTTGCTTTTACC‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAccc                                                                                                              <  2:1187931/148‑1 (MQ=255)                                              tgatCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCAc                                                                                                         >  2:694605/1‑151 (MQ=255)                                                   gtcgCCATGCTACTGCTCGCGGCACCCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTc                                                                                                    >  2:733115/1‑151 (MQ=255)                                                     cgcctttctccggCTCGCGCCCCTCTGTGTTTAATTTGCGCTTTGGGCC‑GGTGGGCGGGCACCAAAAATAATGTATTTCCCCCATGC‑TTTCCC‑AACCC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCa                                                                                                   <  1:127675/138‑1 (MQ=255)                                                     ccccttGCTCCTGCCCGCGGCCCTTTGTTTATTATTTGCGTTTAGTGCC‑TGTGGGCGGCCACCAAAAATAATGCATTTTCCCCTTGTTTTTACC‑AC‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCAt                                                                                                  <  2:78484/146‑1 (MQ=255)                                                      gCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCAt                                                                                                  >  2:565686/1‑150 (MQ=255)                                                       ccATGCTACTGCCCGCGCCCCTCTGGTTATTATTTGCGCTTGGGGCC‑TGTGGGCTGGCACCGAAGATAATGCATATTGCCCTTGCGTTTACA‑AC‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCa                                                                                                   <  1:1291296/148‑1 (MQ=255)                                                           gCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATacac                                                                                              >  2:1131049/1‑149 (MQ=255)                                                              aCTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTc                                                                                         >  1:281452/1‑151 (MQ=255)                                                               cTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTc                                                                                         >  1:1375074/1‑150 (MQ=255)                                                                  ccggcgcccctctGCTTATTATTTGCG‑TTGGTGCCTGGTGGGCGGGCCCCAAAAAAAATCCTTTTCCCCCTTGC‑TTTACC‑AA‑CCCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCa                                                                                        <  1:714122/139‑1 (MQ=255)                                                                  cTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCAccc                                                                                     >  2:338885/1‑151 (MQ=255)                                                                  cTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCAccc                                                                                     >  1:77795/1‑151 (MQ=255)                                                                  cTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCAccc                                                                                     >  2:902823/1‑151 (MQ=255)                                                                    cgcgGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCAccccc                                                                                   >  1:197691/1‑151 (MQ=255)                                                                         ccctctGCTTATTAATTGCGTTTGGTGCC‑TGTTGGCGGGCCCAGAAGAAAATACATATCGCCCATGCTTTCACC‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTc                                                                                 <  1:1068015/146‑1 (MQ=255)                                                                          aCTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCAtt                                                                              >  2:454659/1‑150 (MQ=255)                                                                          aCTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCAtt                                                                              >  2:860713/1‑150 (MQ=255)                                                                            tctgtttATTAATTGCGCTTGGTGCC‑GGTGGGCGGGCACGGAAAATAATGCATATCCCCCTTGCTTTTACC‑AC‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTa                                                                             <  1:522078/143‑1 (MQ=255)                                                                             ctGCTTTTTAATGGCGCTTAGTGCG‑GGTTGGCGGACAGAAAAAAAAATCAATTTTCCCTGTGCTTTTCCC‑AA‑TC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCAtt                                                                              <  2:846538/147‑1 (MQ=255)                                                                              tgcttttaaaTTGCGCTTGGGCCG‑GGTTGGCTGACACAGAAAAAAATCAATTTGCCCCATGTTTTTCCA‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTAc                                                                            <  2:117890/140‑1 (MQ=255)                                                                              tGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGt                                                                         >  2:248777/1‑151 (MQ=255)                                                                               gtatattatTTGCGCTTAGTGCG‑AGGTGGCGGCCACAAAAAATAATGAATATCCCCCATGCTTTTACC‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACt                                                                           <  2:559275/145‑1 (MQ=255)                                                                               gcttattatTTGCGCTTGGGCCC‑GGTGGGCGGCCACAAAAAAAAATGAAT‑TTTCCCCTTCGTTTTCCAAA‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGt                                                                         <  1:770326/147‑1 (MQ=255)                                                                                        ttGCGCTTGGTGCC‑GGTGGGCTGGCACGAAAAAAAATGATTATCCCCCTTGCTTTAACC‑AA‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGt                                                                   <  1:151630/147‑1 (MQ=255)                                                                                        ttGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtact                                                               >  2:1150558/1‑151 (MQ=255)                                                                                        ttGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtact                                                               >  1:292960/1‑151 (MQ=255)                                                                                           cgcTTGGTGCT‑AGTGGGCGGGCACAAAAGATAATGCATATCCCCCATGCTTTTACA‑AA‑TC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtact                                                               <  2:465519/148‑1 (MQ=255)                                                                                              ttggTGCC‑GGTGGGCTGCCACCGAAAAAAATGCATTTTCCCCATGCTTTTACA‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTta                                                                 <  2:277981/140‑1 (MQ=255)                                                                                              tgaGTGCG‑ATTGGGCTGGCACCAAAAATAATGCATATCCCCCATGCGTTCCCC‑AA‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTg                                                           <  2:71373/147‑1 (MQ=255)                                                                                                     c‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCt                                                  >  1:957888/1‑151 (MQ=255)                                                                                                           ggctggcACCGAAGATGATGCATAGTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCCCTGACGCCAGCGGAATAGCACTTGTACCCGCCAGGTTTCATACACATGTCACCCCCCCATGACTGTCGCGGTAACTACTGACGAAGCCTTtgc                                              >  1:242590/1‑150 (MQ=255)                                                                                                           ggcgggcACGGAAAATATTGAATTTTCCCCTTGCTTTACCC‑AC‑TC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtg                                               <  2:664167/149‑1 (MQ=255)                                                                                                             ctggcACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtgctgc                                           >  2:858629/1‑151 (MQ=255)                                                                                                             cgggcACCGAAAATAATTTATTTTGCCCTTGCTTTCCCC‑AC‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtgct                                             <  1:241789/147‑1 (MQ=255)                                                                                                               ggcACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCt                                         >  2:424334/1‑151 (MQ=255)                                                                                                               gccccagAAAATAATGCATTTCCCCCTTGCTTTACCC‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCttt                                                <  1:848167/138‑1 (MQ=255)                                                                                                                 cACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTgg                                       >  1:64725/1‑151 (MQ=255)                                                                                                                 cACCAAAAATAATACATATTGCCCATGCTTTTCCC‑AA‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGATACTACTGGCGATGCCTTtgctgc                                           <  1:1073690/147‑1 (MQ=255)                                                                                                                     gAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCtt                                   >  2:871618/1‑151 (MQ=255)                                                                                                                      aaGATAATGAATATCGCCTTTTCTTTTCCC‑AA‑TC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtgctgc                                           <  1:296072/142‑1 (MQ=255)                                                                                                                           aaTGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTgg                             >  1:461527/1‑151 (MQ=255)                                                                                                                              gCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCtg                          >  2:396892/1‑151 (MQ=255)                                                                                                                                    tGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGgtgt                    >  1:379201/1‑151 (MQ=255)                                                                                                                                      cccTTGCTTTTCCC‑AA‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTgtggtg                     <  1:850773/148‑1 (MQ=255)                                                                                                                                      cccATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCg                  >  2:941237/1‑151 (MQ=255)                                                                                                                                              tttACC‑AC‑CC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGtt                <  2:1005736/145‑1 (MQ=255)                                                                                                                                              tttACC‑AA‑CC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgc             <  2:349702/148‑1 (MQ=255)                                                                                                                                                tACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGTTGTGGTGTCGTTcgcgccgc        >  1:1218346/1‑151 (MQ=255)                                                                                                                                                tACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgcgccgc        >  2:1201451/1‑151 (MQ=255)                                                                                                                                                tACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgcgccgc        >  2:734294/1‑151 (MQ=255)                                                                                                                                                  cc‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgcgccgcgc      >  2:310512/1‑151 (MQ=255)                                                                                                                                                     aa‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgcgccgcgc      >  1:226997/1‑149 (MQ=255)                                                                                                                                                         a‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCg  >  1:110342/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                       CGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCC‑AGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACC‑AA‑TA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCG  >  NC_000913/2603193‑2603491

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: ‑