Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,328	0	C	A	45.2%	64.3 / 46.9	62	intergenic (+337/‑76)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (31/3);  new base A (0/28);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.57e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGG  >  NC_000913/345178‑345468                                                                                                                                                          |                                                                                                                                              taAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATc                                                                                                                                               >  2:732305/1‑151 (MQ=255)  aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCa                                                                                                                                              <  1:754129/151‑1 (MQ=255)  aaaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCa                                                                                                                                              <  1:1123926/151‑1 (MQ=255)             tcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCAt                                                                                                                                  >  2:498602/1‑151 (MQ=255)               aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCAtt                                                                                                                                 >  1:1239117/1‑150 (MQ=255)                tGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCAtt                                                                                                                                 >  2:668208/1‑149 (MQ=255)                   aGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCtc                                                                                                                            >  2:101489/1‑151 (MQ=255)                        ttttGTGTATGTGTTGTTAGCCATTTAAATAAAAAAATTCATAATCCCATCATTTTGGTGGGGGGTCTTTT‑‑TTTTATCTGAGGCAGGGCCCT‑AAACCCGAAAAAAAATAAGGTTAAAAGCCTCAGATAAAAAA‑CCGCCATTAGCtcatca                                                                                                                        <  2:64479/150‑1 (MQ=255)                         tttGTGAATGTTTTGTTTGACTTTCATAAAAAAAAAATCATATTTCCATCTTTTTTGTGGGGGTTTTTTT‑‑CTGTTTTTTCCCCACCCCACCAAAGCCTGAAAAAAAAT‑ATGTTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAg                                                                                                                       <  2:303787/150‑1 (MQ=255)                             tgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGa                                                                                                                     >  1:618319/1‑148 (MQ=255)                               attgtgTTGTTTAAAATTAATAAAAAAAAAATCATAATCCCATCTTTTTTGTGTGGGTTTTTTT‑‑TTTTATCTCCCCCCCCCCACCAAACCCTAAAAAAAAAT‑ATGTTGAAAGGCTAATAAAAAAAA‑CCGCCATTAGCTCATCAGGAaga                                                                                                                  <  1:1019213/147‑1 (MQ=255)                                   ttttgtttgATATTAATAGAAAAAAAATAATAATCCCATCATGTTTGTGTGGGTTTTTTT‑‑CTGTATCTTACCCACCCCCCTCAAGCCTAAAAAAAAAT‑ATGCTGTAAGCCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCa                                                                                                               <  1:96733/146‑1 (MQ=255)                                   tgtTGGTTGCCTTTAATAAAAAAAAAATCAAAATCCCATAGTTTTGGGGTGGGTTTTTTC‑‑TTTTATCTTCCCCACCCCCCCAACCCCAAAAAAAAAAA‑ATGGTAAAGGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAg                                                                                                              <  1:66239/149‑1 (MQ=255)                                     ttggttaACATTAAAAAAAAAAAAAACAAAATACCATCATGTTTGTGGGGTTTCTTTT‑‑CTTTTTTTTACGCACCGCACTAAGGCTTAAAAAAAAAT‑ATGTGGAGAGACTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAg                                                                                                              <  1:1149276/147‑1 (MQ=255)                                       ggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                         >  2:787998/1‑150 (MQ=255)                                          tgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                     >  1:407259/1‑151 (MQ=255)                                            aCATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                     >  2:293459/1‑149 (MQ=255)                                              aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                  >  1:414982/1‑150 (MQ=255)                                               ttaatatAAAAAAAATCATAATCCCATCATGTTTGTGTGGGGTTTTTT‑‑CTTTATCTTACGCATCCCACCCAAGCCAAAAAAAAAAT‑ATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGa                                                                                                 <  1:1116968/147‑1 (MQ=255)                                               ttaaaaaaaaaaaaaTAAAACTCCCTTCTTGTTTGGGGGGGTTTTTTT‑‑CTGTATCTTCCCAATCGCCCCAAACCCAAAAAGAAAAA‑AGTTTGTAAGCCTCAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAg                                                                                                   <  1:127871/140‑1 (MQ=255)                                               ttCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGGGCAGACGACAACCAGGAt                                                                                                >  1:516277/1‑151 (MQ=255)                                               ttCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt                                                                                                >  1:651459/1‑151 (MQ=255)                                                  atatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                                             >  1:899102/1‑151 (MQ=255)                                                      aaaaaaaaTCATAATCCCATCTTTTTTGGGTGGGTTTTTTT‑‑TTTTATCTTCCCCACCCCACCAAACCCTAAAAAAAAAT‑GTGGTGAAAGCCAAAAATAAAAAT‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                                             <  1:138249/146‑1 (MQ=255)                                                       aaaaaacACATAATTCCAT‑AT‑TTTTTGTGGGATCTTTTTCCTTTATCTAACCCCCCGCCCAAACCCCTGACAAAAAAT‑ATGTTGAAAGTTTCAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                          <  2:622083/149‑1 (MQ=255)                                                       aaaaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                        >  1:608223/1‑151 (MQ=255)                                                       aaaaaaaCCAAAATCCCATCATTTTTGGGTGGGTTTTTTT‑‑CTTTATCTTCCACATCCCACACACCCCGAAAAAAAAAA‑ATGGTTGAAGGCACATATCAAAAC‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATtgttgt                                                                                          <  1:177156/149‑1 (MQ=255)                                                          aaaatcaaaatcCCATCATTTTTGGGTGGGTTTTTTT‑‑TTGTTTTTTACCCACCGCACCAAACCCGAAAAAAAAAT‑TTGTTGAAAGCCAAAAAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                     <  1:1194169/151‑1 (MQ=255)                                                                   aTTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCCTCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTc                                                                            >  1:1232048/1‑151 (MQ=255)                                                                   aTTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTc                                                                            >  1:1017636/1‑151 (MQ=255)                                                                      ccATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                                           >  2:900856/1‑149 (MQ=255)                                                                        atcatGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc                                                                       >  1:702857/1‑151 (MQ=255)                                                                        atcatGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc                                                                       >  2:525393/1‑151 (MQ=255)                                                                        atcatGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc                                                                       >  1:194446/1‑151 (MQ=255)                                                                         tcatTTTTGGGGGGGTTTTTTT‑‑CTGTATCTTACACACCCCACTAAGGCCGAACAAAAAAT‑ATGCTGAAGGGCTAATAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                                        <  2:299142/149‑1 (MQ=255)                                                                          catGTTTGTGGGGGGTTTTTT‑‑CTGTATTTTCCGCATCGCACTCAAGCCTGACAAAAAAT‑ATGCTGAAAGCCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                                        <  2:302562/148‑1 (MQ=255)                                                                           atGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                    >  2:626875/1‑151 (MQ=255)                                                                            tGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                   >  2:255872/1‑151 (MQ=255)                                                                            tGTTTGTGGGGGTTCTTTT‑‑CTTTATCTCACCCATCGCACTCAGCCCTAACAAAAAAT‑ATGCTGTAAGGCAAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                    <  2:1249611/150‑1 (MQ=255)                                                                             gTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                   >  1:328526/1‑150 (MQ=255)                                                                              tttGTGTGGGTTTTTTT‑‑CTGTTTCTTCCGCATCCCACCAAACCCTAAAAAAAAAT‑ATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGACTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                    <  1:393848/148‑1 (MQ=255)                                                                                  tgtgGGTTCTTTT‑‑CTATATCTACCACGCCCCACACGAGCCAGATAAAAAAT‑ATGCATTAGGGCT‑ATATAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                               <  1:519424/149‑1 (MQ=255)                                                                                   gtgGGTTTTTTT‑‑CTGTATCTTCCCCACCGCACTCAAGCCGGAAAAAAAAT‑ATGCTGTAGGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                                              <  2:135423/149‑1 (MQ=255)                                                                                     gggtTCTTTT‑‑CTGTATTTTCCGCACCCCACTCAACCCAGACAAAAAAT‑ATGCTGTAAGGCTCAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                           <  1:408655/150‑1 (MQ=255)                                                                                           tttt‑‑CTGTATCTTCCGCATCGCACTCAAGCCTGAAAAAAAAT‑ATGCTGAAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                                     <  2:105580/150‑1 (MQ=255)                                                                                           tttt‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                    >  2:1192012/1‑151 (MQ=255)                                                                                           tttt‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt                                                    >  1:728065/1‑151 (MQ=255)                                                                                                              aTCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                    >  1:914628/1‑150 (MQ=255)                                                                                                                cccacacaACCCTGAAAAAAAAT‑ATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                 <  2:582457/149‑1 (MQ=255)                                                                                                                   aCTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                 >  1:246293/1‑148 (MQ=255)                                                                                                                   aCTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                 >  2:884538/1‑148 (MQ=255)                                                                                                                    cTCAAGCCTAAAAAAAAAT‑ATGCTGAAAGGCTCATAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAg                               <  2:497849/149‑1 (MQ=255)                                                                                                                       agcccaaaaaaaaaaT‑ATGCTGAAAGCCTAATAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                <  2:94712/134‑1 (MQ=255)                                                                                                                       aaccctaaaaaaaaaT‑ATGCTGAAAGCCTCATAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCa                                <  2:632113/134‑1 (MQ=255)                                                                                                                       aaGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                          >  1:1089546/1‑151 (MQ=255)                                                                                                                            tGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttat                     >  2:451750/1‑151 (MQ=255)                                                                                                                               aaaaaaaT‑ATGTTGTAAGATTTATAAAAAAAG‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt                    <  2:507875/146‑1 (MQ=255)                                                                                                                                aaaaaaT‑ATGCTGTAGGTCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                 <  1:635266/149‑1 (MQ=255)                                                                                                                                 gAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatc                 >  1:1239412/1‑150 (MQ=255)                                                                                                                                              tAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt    >  2:958167/1‑151 (MQ=255)                                                                                                                                              aaaGGTTCAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg          <  2:806071/144‑1 (MQ=255)                                                                                                                                               aaGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt    <  2:317387/150‑1 (MQ=255)                                                                                                                                                      atatAAAAAG‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg  <  1:424370/145‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                              TAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTT‑‑CTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAAT‑ATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGG  >  NC_000913/345178‑345468

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