Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	40.6%	126.3 / 11.5	69	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (5/36);  new base G (28/0);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.41e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAC  >  NC_000913/3133156‑3133447                                                                                                                                                   |                                                                                                                                                       gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTc                                                                                                                                                     <  2:83955/150‑1 (MQ=255) gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTc                                                                                                                                                     <  1:1346166/150‑1 (MQ=255)        ggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGAGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCacgc                                                                                                                                             >  1:614956/1‑148 (MQ=255)          cTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATcc                                                                                                                                             >  1:364530/1‑149 (MQ=255)                cgcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGGCCCTCC‑CGt                                                                                                                                         >  2:1451453/1‑146 (MQ=255)                 gcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAgg                                                                                                                                   <  1:26535/151‑1 (MQ=255)                  cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑AGTTCaa                                                                                                                                     >  1:243313/1‑148 (MQ=255)                   ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGAt                                                                                                                                 <  1:1257763/151‑1 (MQ=255)                   ggTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑ACTTCaaag                                                                                                                                   >  2:889449/1‑147 (MQ=255)                    gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGGCCATCC‑CGGTTAAGGAtt                                                                                                                                >  2:778192/1‑150 (MQ=255)                    gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCATCC‑CGTTCaaaggt                                                                                                                                 >  2:1127766/1‑146 (MQ=255)                     ttGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑CGTTCAAAGATCt                                                                                                                               >  2:22285/1‑151 (MQ=255)                      tGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑CGTTCAAGGATCtt                                                                                                                              >  1:252599/1‑151 (MQ=255)                       gTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCCGGGGGTCGCCCATCC‑CGTTCAAAGAACttt                                                                                                                             >  1:806567/1‑151 (MQ=255)                         aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTg                                                                                                                            <  1:999173/150‑1 (MQ=255)                         aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCAGCCCCTAGCCAGTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTgc                                                                                                                           <  2:348884/151‑1 (MQ=255)                              tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCTGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCAt                                                                                                                      >  2:961589/1‑151 (MQ=255)                              tAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑CGTTTAAGGATCTTTTCGccct                                                                                                                      >  1:478952/1‑149 (MQ=255)                               aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCAtt                                                                                                                     <  2:781100/151‑1 (MQ=255)                               aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCAtt                                                                                                                     <  1:1322981/151‑1 (MQ=255)                               aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCAtt                                                                                                                     <  1:916988/151‑1 (MQ=255)                                aCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCC‑AGCTTAAAGATGTTTTTGCCAtt                                                                                                                     >  1:432656/1‑150 (MQ=255)                                  agtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑CGTTCAAGGATCTTTGCGCCAtt                                                                                                                     >  1:198152/1‑148 (MQ=255)                                     ttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATg                                                                                                                <  2:497272/150‑1 (MQ=255)                                      tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑AGTTCAAGGAACTTTTCGCCATTgta                                                                                                                  >  1:780145/1‑145 (MQ=255)                                       gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCa                                                                                                             <  1:370777/151‑1 (MQ=255)                                               aGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcgg                                                                                                      <  1:35081/150‑1 (MQ=255)                                                  gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑AGTTTATAGAGTTTTGCGCCCTTGAATGTCAGTTgcgg                                                                                                      >  2:1318714/1‑147 (MQ=255)                                                      ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACgg                                                                                              <  1:1139647/151‑1 (MQ=255)                                                           cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑Att                                                                                       <  1:870313/151‑1 (MQ=255)                                                            gTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑Attt                                                                                      <  2:1202178/151‑1 (MQ=255)                                                             ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑AAGTTAAAGATGTTTTTGGCAATTTATTTCAGTTGCGGCGGCAAGGCG‑‑Gtttt                                                                                     >  2:662002/1‑151 (MQ=255)                                                                aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGACCATTC‑CGCTCAAAGGTCTTTGCGCCATTTAATGTCCGTTTCGGCGTAACAGCG‑‑Attttttt                                                                                  >  1:1408017/1‑150 (MQ=255)                                                                   catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCt                                                                                <  1:187691/150‑1 (MQ=255)                                                                    atcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGt                                                                              <  1:1016527/151‑1 (MQ=255)                                                                        ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑CGTTCAAGGGTCTTTGCGCCATTTAATGTCAGTTGCGGCGGAACGGCG‑‑GTTTTTTTCtttta                                                                            >  1:907292/1‑145 (MQ=255)                                                                         cTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑CGTTTAAGGATCTTTTCGCCCTTTAAAGTCAGTTGCGGCGGAACGGCG‑‑ATTTTTTTCTTTAACGc                                                                         >  1:1200288/1‑151 (MQ=255)                                                                           aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTTTAACGCt                                                                        <  1:1157416/150‑1 (MQ=255)                                                                           aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCt                                                                        <  1:380091/150‑1 (MQ=255)                                                                           aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCC‑CGTTTCAGGATCTTTTCGCCCTTTAATGTTAGTTGCGGGGTGACGGCGC‑ATTTTTTCCTTTaagc                                                                          >  2:237491/1‑147 (MQ=255)                                                                              cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCa                                                                    <  1:1185319/151‑1 (MQ=255)                                                                                       gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTGGCCCCTCC‑ACTTCAAGGATCTTTGCGCCATTTAATGGTAGTTGCGGGGGGACGGGG‑‑GTTTTTTTCTGTAACGATGCCC‑GCtgggggg                                                          >  2:1338359/1‑150 (MQ=255)                                                                                          cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGGCCGAAC‑ATGTCCAGGTTCTTTGTGGCAATGATTGACATGTGGGGGGTGACGGCGGTATTTTTTTCTG‑AGGGCCGCCA‑GTtgggg                                                            >  2:1127480/1‑145 (MQ=255)                                                                                             gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGtt                                                    <  1:1258542/151‑1 (MQ=255)                                                                                              aaaGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGttt                                                   <  1:1453095/151‑1 (MQ=255)                                                                                                aGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCCTCC‑AGGTCAAGGGTCTTTTCGGCACTGCAAGTTAGATTCGGCGTTACGGGG‑‑GGTTTTTTCTGTAACACCGCCC‑GTttgtgggtcggggt                                                   >  1:1102683/1‑141 (MQ=255)                                                                                                 gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTAcg                                                <  2:1423232/151‑1 (MQ=255)                                                                                                  cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCACC‑CGTTCAAAGATCTTTTCGCCCTTTAATGTCAGTTGCGGGGGAAAGGCG‑‑GTTTTTTTCTGTTAAGCTGGCAGGTTGGGGGCCGGGtt                                                    >  2:1141346/1‑147 (MQ=255)                                                                                                     ggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCa                                           <  1:644714/151‑1 (MQ=255)                                                                                                      gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑CGTTCAAAGATCTTTGCGCCATTTAATGTCAGTTTCGGGGGAAAGGCG‑‑GGTTTTTTCTGTAAACCTGCCC‑GTTGGTGGGCAgggt                                                    >  1:544838/1‑140 (MQ=255)                                                                                                      gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCaa                                          <  2:1402673/151‑1 (MQ=255)                                                                                                           aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑AGTTTCAAGATATTTTCGCCATTTTATGTCCAGTGCGGCGCAACGGCG‑‑GTTTTTTCCTTTAAATCTGGCA‑GTTGGTGGGCAGGGTTTCG‑CGCaaa                                         >  1:512421/1‑147 (MQ=255)                                                                                                               gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGttt                                 <  1:1450864/151‑1 (MQ=255)                                                                                                               gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGttt                                 <  2:589958/151‑1 (MQ=255)                                                                                                                cgcgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTg                                <  2:1427811/151‑1 (MQ=255)                                                                                                                   gTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTtc                             <  1:96133/151‑1 (MQ=255)                                                                                                                     aGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCC‑AGTTTAAGGATCTTTGCGCCATTTAATGTCAGTTGCGGCGGAACGGCG‑‑GTTTTTTTCTGTAACGCTGCCC‑GGTGGTGGGCAGGGTTTCG‑CGCAAACGCGGGTTGTttat                           >  1:1211740/1‑151 (MQ=255)                                                                                                                       tGACGCTTGCAGAATTTCATCGGCGGGGGGTCTCCCATCCTTGTTTAA‑GATTTTTGCGCCCTTTAATGTTTGTTGCGGCGCTACCGGG‑‑GTTTTTTTCTTTTATACTGCCA‑GTTTTTGGAGAGGTTTATT‑CGGAAAATCtgtttgtt                              >  2:1246652/1‑146 (MQ=255)                                                                                                                        gACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTtcatca                         <  2:720159/150‑1 (MQ=255)                                                                                                                        gACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCAACC‑CGCTCCAGGAGCTTTGCGCCCTTGAATGTCAGGTGGGGGGCGACGGGG‑‑AACTTTTTCTTTGACGCTGCCA‑GGTGGGGGGGAGGGTTTCG‑CGCAAACGCGGGTTGTttatt                          >  2:1427268/1‑148 (MQ=255)                                                                                                                         aCGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAAt                       <  1:1325569/151‑1 (MQ=255)                                                                                                                           gCTTGCAGAATTTCATCGGCGGGGGGTTGGCCCTCC‑AGATCAAAGGTCTTTTCGCCACTGGATGTTAGTTGCGGGGGAACGGCG‑‑GTTTTTTTCTTTTACGCTGCCA‑CGTGGGGGTTAGGGTTTCGCCCCAAACCCGGttgtttat                            >  1:1409263/1‑139 (MQ=255)                                                                                                                              tGCAGAATTTCATCGGCGGGGGGTCGCCCCTCC‑ACTTTCAGGATCTTTTCGCCCTTGTATGTCAGTTGCGGCGGAACGGCG‑‑GTTTTTTTCTTTTACGCCGCCA‑GTTGGGGGTGAGGGTTTCG‑CGCAAAAGCGGGTTGTTTTTCaaaggc                    >  2:74776/1‑145 (MQ=255)                                                                                                                               gCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAg                 <  2:1299203/151‑1 (MQ=255)                                                                                                                                 aGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAGAc               <  2:757158/151‑1 (MQ=255)                                                                                                                                      ttCAGCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTTTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCATGACGAAt           <  1:1201718/150‑1 (MQ=255)                                                                                                                                         aTCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaac       <  2:243310/151‑1 (MQ=255)                                                                                                                                          tCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGTG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaac       >  1:1359964/1‑150 (MQ=255)                                                                                                                                              cggGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc  <  2:1243832/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCC‑AGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCG‑‑ATTTTTTCCTGTAACGCTGCCA‑GTTGGTGGTCAGGTTTACG‑CGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAC  >  NC_000913/3133156‑3133447

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