Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,304	0	G	A	44.1%	88.2 / 38.6	68	intergenic (+313/‑100)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (30/8);  new base A (0/30);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.67e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTT  >  NC_000913/345157‑345447                                                                                                                                                       |                                                                                                                                                 gTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGaa                                                                                                                                                >  1:218825/1‑150 (MQ=255)      gCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATg                                                                                                                                          <  1:516439/151‑1 (MQ=255)      gCAGCCATGCGGTGCATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATg                                                                                                                                          >  1:924017/1‑151 (MQ=255)          ccATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGt                                                                                                                                      <  1:694378/151‑1 (MQ=255)           cATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGt                                                                                                                                      >  2:1124839/1‑150 (MQ=255)              gCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAgg                                                                                                                                  <  1:550733/151‑1 (MQ=255)               cGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGc                                                                                                                                 <  2:448690/151‑1 (MQ=255)                        aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                         >  1:1133342/1‑150 (MQ=255)                        aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                         >  1:1158376/1‑150 (MQ=255)                        aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCa                                                                                                                         >  1:149982/1‑150 (MQ=255)                         aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaa                                                                                                                       <  2:67693/151‑1 (MQ=255)                         aaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaa                                                                                                                        >  2:424418/1‑150 (MQ=255)                          aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaaa                                                                                                                      <  1:730320/151‑1 (MQ=255)                                  tcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCAt                                                                                                             >  2:22937/1‑151 (MQ=255)                                    aggCAGTTTTTTTGTAAATGTGTTGGTTGACATTCATATAAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGT‑TTTTTTCTGTACCTTACGCATCCCACTAAGGCCTACCAGAAAATATGCTGTAAGGCTCTTATCAAAAA‑CCGCCAt                                                                                                             <  2:722363/147‑1 (MQ=255)                                        aGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCtc                                                                                                       >  2:1114950/1‑151 (MQ=255)                                           tGTTTTGTAAATGTGTTGTTTGACATTCATATGAAAAAAATCATAATC‑CCACTAT‑TTTTGTGGGGGT‑TTTTTTCTGATTCTTCCCCTCCCCACCAAAGCTTAACAAAAAATATGCTGTAAGGCTAATATCAAAAA‑CCGCCATTAGCtca                                                                                                      <  2:808709/149‑1 (MQ=255)                                            gTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCtcatc                                                                                                    >  2:417328/1‑150 (MQ=255)                                                 gtgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGa                                                                                                >  2:463005/1‑149 (MQ=255)                                                        ttttgtttgACATTCATATAAAAAAAATCATAATC‑CCATCTT‑TTTTGTGGGGGG‑TCTTTTCTGTATCTTCCGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAg                                                                                         <  1:836350/147‑1 (MQ=255)                                                          ttggttgACATTCATATAAAAAAAATCAAAATC‑CCATCAT‑TTTTGTGGGGGG‑TTTTTTCTGTATCTTACGCACCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCa                                                                                          <  1:387276/146‑1 (MQ=255)                                                          ttggttgACATTCATATAAAAAAAATAAAAATC‑CCATCAT‑GTTGGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAgacg                                                                                      <  1:112149/150‑1 (MQ=255)                                                          tgggtagACATTCACAGGAAAAAAATAACAACC‑CCA‑CATGGTTTGTGAGGGG‑TTTTTCGCGTATCTGCCGCACCCCACTCAAGCCTGACAAAAAATTTGCTGTAAGTCACA‑ACCAAAAACCCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                       <  1:1121673/147‑1 (MQ=255)                                                              ttgACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                 >  2:916183/1‑151 (MQ=255)                                                              ttgACATTCATATAAAAAAAATCATAATC‑CCATCAT‑GTTGGTGGGGGG‑TCTTTTCTGTATCTTACGCATCCCACTCAAGCCTGACAAAAAATATGCTGTAGGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACaa                                                                                  <  1:234847/150‑1 (MQ=255)                                                                     tCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGGCAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                          >  1:671460/1‑151 (MQ=255)                                                                     tCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt                                                                          >  1:144136/1‑151 (MQ=255)                                                                      cATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttg                                                                         >  2:619941/1‑151 (MQ=255)                                                                       atatGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                        >  2:1119444/1‑151 (MQ=255)                                                                          aaaaaaaaaaCAAAATTCCCAATAT‑GTTTGTGTGGGT‑TTTTTTCTGTTTCTTACGCACCCCCCTCAACCCTAAAAAAAAATAGGCTGAAAGGCTAAAAACAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt                                                                           <  2:778313/144‑1 (MQ=255)                                                                            aaaaaaaTCATAATC‑CCATCAT‑TTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa                                                                    <  1:1030147/150‑1 (MQ=255)                                                                            aaaaaaaTAAAAATC‑CCATCAT‑TTTTGTGGGGGT‑TTTTTTCTGTATCTTACGCACCCACCCCACGCCGGAAAAAAAATATGCTGAAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgtt                                                                       <  2:317978/147‑1 (MQ=255)                                                                                 aaTCATAATC‑CCATCTT‑GTTTGTGGGGGT‑TCTTTTCTTTTTCTTACGCATCGCCCTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                              <  1:1104075/151‑1 (MQ=255)                                                                                  aTCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                             >  1:849802/1‑151 (MQ=255)                                                                                     aaaaTC‑CCATCTT‑GTTTGTGGGGGG‑TTTTTTCTGTTTTTTCCGCACCGCACTCAAGCCGGAAAAAAAATATGCTGAAAGGCTAAAACAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                             <  2:853387/146‑1 (MQ=255)                                                                                       aaTC‑CCATCTT‑GTTTGTGGGGGT‑TCTTTTCTTTATCTTCCGCATCCCACTCACGCCTGACAGAAAATATGCTGTAAGCCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                        <  2:1012190/151‑1 (MQ=255)                                                                                       aaTC‑CCATCAT‑TTTTGGGTGGGT‑TTTTTTCTGTATCTTCCGCACCGCACTCAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                          <  1:956376/149‑1 (MQ=255)                                                                                        aTT‑CCAT‑AT‑GTTTGTGGGGGGTTTTTTTCTTTTTTTTCCGCACCCCACCCAAGCCTGAAAAAAAATATGCTGAAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                        <  2:832038/150‑1 (MQ=255)                                                                                            ccATCAT‑TTTTGTGGGGGG‑TTTTTTCTGTATCTCACGCATCGCACCCAAGCCGCAAAAAAAATATGCTGAAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                        <  1:261550/147‑1 (MQ=255)                                                                                            ccATCAT‑TTTGGTGGGGGT‑TTTTTTCTGTATCTTACGCACCGCACCCAACCCTGACAAAAAATAGGCTGAAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                        <  2:791141/147‑1 (MQ=255)                                                                                            ccATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                      >  2:662106/1‑149 (MQ=255)                                                                                            ccATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg                                                      >  1:83701/1‑149 (MQ=255)                                                                                            ccATCAT‑GTTTGTGGGGGT‑TTTTTTCTTTATCTTACGCACCCCACGAAGGCCGGACAAAAAATATGCTGTAAGGCTAATACCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                    <  1:753842/151‑1 (MQ=255)                                                                                            ccATCAT‑GTTGGTGGGGGT‑TCTTTTCTGTATCTTCCCCATCGCACTCAAGCCTGAAAAAAAATATGCTGAAAGGCTAAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                    <  2:366929/151‑1 (MQ=255)                                                                                             catcat‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                   >  2:1142992/1‑151 (MQ=255)                                                                                                        gtgtgGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGc                                         >  2:114873/1‑151 (MQ=255)                                                                                                         tgtgGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCg                                        >  2:543995/1‑151 (MQ=255)                                                                                                          gtgGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCgg                                       >  1:1089670/1‑151 (MQ=255)                                                                                                            gggg‑TCTTTTCTTTATCTTCCCCACCCCACTCAACCCTAACAAAAAATATGCTGTAGGCCTAAAAACAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                          <  2:986583/146‑1 (MQ=255)                                                                                                               gttCTTTTTTGTATCTTCCGCACCGCACAAAAGCCTGACAAAAAATATGCTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                   <  1:122009/149‑1 (MQ=255)                                                                                                               g‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                   >  1:116854/1‑150 (MQ=255)                                                                                                                 tcttttctttATCTTACGCACCCCACTCAACCCTGAAAAAAAATATGCTGTAAGGCTAATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                   <  2:726694/149‑1 (MQ=255)                                                                                                                 tCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                 >  1:1106422/1‑151 (MQ=255)                                                                                                                   ttttCTGTATCTTACGCATCGCACTCAAGCCTGAAAAAAAATATGCTGTAAGGCTAATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACt                                <  1:1025338/150‑1 (MQ=255)                                                                                                                      tCTGTATCTTCCGCACCCCACTCAAGCCTGAAAAAAAATATGTTGTAAGGCTCATATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                             <  2:828434/150‑1 (MQ=255)                                                                                                                            tCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                      >  1:386326/1‑151 (MQ=255)                                                                                                                              ttACGCACCGCACCCAAGCCTGACAAAAAATATTCTGAAAGCCCCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                     <  1:414825/150‑1 (MQ=255)                                                                                                                                 cGCATCGCACTCAAGCCTAACAAAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtctt                   <  1:450432/149‑1 (MQ=255)                                                                                                                                 cGCATCGCACCCAAGCCTAAAAAAAAATATGCTGTAAGGCTCAAATAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAActc                     <  1:225118/147‑1 (MQ=255)                                                                                                                                      cGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc            >  2:609688/1‑151 (MQ=255)                                                                                                                                         aCTCAAGCCTGAAAAAAAATATGGTGAAAGACACATACAAAAAC‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAg          <  2:154125/150‑1 (MQ=255)                                                                                                                                          cTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc            >  1:861743/1‑147 (MQ=255)                                                                                                                                            cAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc      >  2:279035/1‑151 (MQ=255)                                                                                                                                             aaGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc     >  1:788044/1‑151 (MQ=255)                                                                                                                                             aaGCCTGAAAAAAAATATGCTGAAAGGCTAATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt        <  1:44574/148‑1 (MQ=255)                                                                                                                                             aaGCCTAAAAAAAAATAGGCTGAAAGGCTAAAATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATc       <  2:265284/149‑1 (MQ=255)                                                                                                                                                ccTGAAAAAAAATATGCTGTAGGCCTCATAAAAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt  <  2:135272/151‑1 (MQ=255)                                                                                                                                                ccTGAAAAAAAATATGCTGAAGGGCTCATATAAAAGA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt  <  2:1053169/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 GTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATT‑CCATCAT‑GTTTGTGTGGGG‑TCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAA‑CCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTT  >  NC_000913/345157‑345447

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