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

AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTG  >  NC_000913/345180‑345475                                                                                                                                                       |                                                                                                                                                    aaaaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCaa                                                                                                                                                   <  2:215846/151‑1 (MQ=255)   aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCaaaa                                                                                                                                                 <  2:767321/151‑1 (MQ=255)    aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCaaaaa                                                                                                                                                <  2:228154/151‑1 (MQ=255)       gATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAAcc                                                                                                                                              <  2:200283/150‑1 (MQ=255)        aTCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGc                                                                                                                                            >  2:723429/1‑151 (MQ=255)           tcATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCAt                                                                                                                                         >  2:568271/1‑151 (MQ=255)             aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                                        >  2:326601/1‑150 (MQ=255)             aTGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCAtt                                                                                                                                        >  1:349601/1‑150 (MQ=255)              tGCAGATGTTTTGTGATTGTGTTGGTTGACTTTCATATGAAAAAAACCATAATCCCATCT‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCACCCCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATAAAAAACCGCCAt                                                                                                                                         <  2:54365/148‑1 (MQ=255)                cAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCt                                                                                                                                    >  2:548123/1‑151 (MQ=255)                cAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCt                                                                                                                                    >  2:81201/1‑151 (MQ=255)                    ttttttGTGATTGTTTTGGTTGACATTCATATAAAAAAAATCATAATCCCATCA‑TGTTTGTGTGGGGTTTTTTCTGTATCTTCCGCATCGCACTCAACCCT‑AACAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCtca                                                                                                                                  <  1:615490/147‑1 (MQ=255)                    ttttttGTGAATGTGTTGGTTGACATTCATATAAAAAAAATCATAATCCCATCA‑TGTTTGTGGGGGTTTTTTTCTGTATCTTACGCACCGCCCTCAAGCCT‑GACAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCt                                                                                                                                    <  1:503256/145‑1 (MQ=255)                     gTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                >  1:80924/1‑150 (MQ=255)                             attgttTTGGTTGACTTTAATAAAAAAAAAAGCATATCCCCAATA‑TTTTGGTGGGGGTTTTTTTTTTTATCTAACGCGCCGCCATCAAGCTG‑AAAAAAAAATATGGTGAAAGCCTAACACAAAAAGCCGCCATTAGCTCATCAGGAagag                                                                                                                        <  2:622626/148‑1 (MQ=255)                                 tgtTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                                   >  2:124135/1‑151 (MQ=255)                                   ttgtttgACTTTTATATAAAAAAAATCAAAATCCCATCA‑TGTTTGTGGGGGGTTTTTTTTGTATCTTACGCACCCCACCCACCCTT‑AAAAAAAAATATGTTGTAGGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgac                                                                                                                   <  2:140218/149‑1 (MQ=255)                                   tattgtcatttttttataAAAAAAAAAATAA‑TCCGTCATTTTTTGGGGGGGTTTTTTTCGGTTTCTTCCGCACCCACCGAAAGCCT‑AACAAAAAATATGCTGAGGGGCGAATATAAAAAGCCGCCATTAGCTCATCAGGAaga                                                                                                                         <  2:685842/129‑1 (MQ=255)                                     gtttaaCATTAATATAAAAAAAACCATAATCCCATAT‑TTTTGGTGGGGGGTTTTTTCTGTATCTCACGCACCCCACCCAGCCCT‑GAAAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                               <  1:579534/146‑1 (MQ=255)                                     ggttgACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac                                                                                                                >  1:925139/1‑150 (MQ=255)                                       ttgACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                             >  1:535126/1‑151 (MQ=255)                                       ttgACATTAATATAAAAAAAATAATAATCCCATCA‑TGTTTGTGGGGGTTTTTTTCTGTTTTTTACGCACCCCACTCAAGCCT‑GACAAAAAATATGCTGTAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACa                                                                                                               <  1:977094/149‑1 (MQ=255)                                         gACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                            >  1:596576/1‑150 (MQ=255)                                                aaataaaaaaaaTCACAGCCCTAACA‑TTTTTGTGTGGTGTCTTTTTTGAATCTTCCCCTTCTCTCCCAAGACT‑TAAAAAAAATATTCAGAAATGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttg                                                                                                     <  2:698418/145‑1 (MQ=255)                                                 tatGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgtt                                                                                                   >  2:3214/1‑151 (MQ=255)                                                 tatGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgtt                                                                                                   >  1:133143/1‑151 (MQ=255)                                                     aaaaaaaTCATAATCCCATCA‑TTTTTGTGTGGGGTTTTTTCTTTATCTTACGCATCGCACTCAAGCCG‑GAAAAAAAATATGCTGTAGGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                               <  1:506603/151‑1 (MQ=255)                                                         aaaCCATATTCCCATCA‑TGTTTGTGGGGGTTTTTTTCTGTATCTTACGCACCGCACCCACGCCT‑GAAAGAAAATATGCTGTAAGCCTCATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTa                                                                                           <  2:898086/151‑1 (MQ=255)                                                          aaTCATAATCCCATCT‑TTTTTGTGGGGGTTTTTTTCTGTATCTTACCCACCGCACCCAAGCCA‑AAAAAAAAATATGCTGTAAGGCACATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc                                                                                          <  2:642221/151‑1 (MQ=255)                                                                aaTTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                    >  2:553660/1‑151 (MQ=255)                                                                    ccATCA‑TGTTTTTGGGGGGTTTTTTCTGTATCTTACGCACCCCACTAAGGCCT‑AAAAAAAAATATGCTGTAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa                                                                                 <  1:234737/150‑1 (MQ=255)                                                                       tca‑tTTTGGTGGGGGTTTTTTTCTGTATTTTCCGCACCGCACTCAAGCCT‑GAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                                                <  2:349662/148‑1 (MQ=255)                                                                             ttggtgggggtgtttttCTATTTCTTCCGCCCCGCAC‑CAAACCTAGAAAAAAAATATGCTGTAAGGCTTAAAACAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGa                                                                         <  2:317816/137‑1 (MQ=255)                                                                             ttggggtgGGTTTTTTTTTTTATCTTCCGCACCCCACTCAACCCT‑GAAAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                          <  2:202779/144‑1 (MQ=255)                                                                               tgtgtgGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                      >  1:325864/1‑151 (MQ=255)                                                                                gtgggGGTTTTTTTCTGTTTCTTACGCATCCCACTAAGGCCT‑GAAAAAAAATTTGCTGAAGGCCTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATg                                                                       <  2:79182/149‑1 (MQ=255)                                                                                gtgggGGTTCTTTTCTGTATCTTACGCACCGCACTCGAGCCG‑AACAAAAAATATGCTGTAGGGCTATTATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGAt                                                                        <  2:246768/148‑1 (MQ=255)                                                                                      gtttttttCTGTATCTTCCGCACCGCCCTAAAGCCT‑GAAAAAAAATATGCTGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa                                                               <  1:902808/147‑1 (MQ=255)                                                                                      ggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGATGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                                 >  2:886417/1‑149 (MQ=255)                                                                                        ttttttCTGTATCTTACGCTTCGCACCCAAGCCT‑GAAAAAAAATAGGCTGTAAGGCTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCaa                                                              <  1:266254/148‑1 (MQ=255)                                                                                            ttCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                         >  2:42190/1‑151 (MQ=255)                                                                                               tGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCa                                                       >  1:748983/1‑150 (MQ=255)                                                                                                gTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaa                                                      >  1:425988/1‑150 (MQ=255)                                                                                                   tCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct                                                  >  1:493580/1‑151 (MQ=255)                                                                                                          caccccccTCAAGTCC‑AACAAAAAATATGTTTTAAGGCTCATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                           <  1:258669/144‑1 (MQ=255)                                                                                                          caccGCACTAAAGCCT‑GACAAAAAATATGTTGTAAGCCTAATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACAACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                           <  1:150941/148‑1 (MQ=255)                                                                                                             cGCACACAAGCCT‑GAAAGAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttct                                             <  1:703850/146‑1 (MQ=255)                                                                                                                  gaaggCCT‑GACAAAAAATATGCTGTGGTGCTAAAAACAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                                    <  1:497194/146‑1 (MQ=255)                                                                                                                   cAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                  >  2:153918/1‑151 (MQ=255)                                                                                                                   cAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcc                                  >  2:797045/1‑151 (MQ=255)                                                                                                                    aaGCCT‑GAAAAAAATTATGCTGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                            <  1:696293/140‑1 (MQ=255)                                                                                                                    aaGCCT‑GAAAAAAAATATGCTGAAAGGTTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                                    <  1:125843/148‑1 (MQ=255)                                                                                                                    aaGCCT‑AAAAGAAAATATGTTGTAAGTCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACtcttc                                              <  2:320947/138‑1 (MQ=255)                                                                                                                       ccT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                              >  1:114509/1‑151 (MQ=255)                                                                                                                       ccT‑GAAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                              <  2:655751/151‑1 (MQ=255)                                                                                                                               aaaaaTATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatca                       <  2:149661/150‑1 (MQ=255)                                                                                                                                aaaaTATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcat                      <  1:908485/151‑1 (MQ=255)                                                                                                                                        cTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGCACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa               >  2:92126/1‑150 (MQ=255)                                                                                                                                         tgaagggCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGc                <  1:182417/143‑1 (MQ=255)                                                                                                                                         tGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa              <  2:456249/150‑1 (MQ=255)                                                                                                                                           tAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGt           >  1:704417/1‑151 (MQ=255)                                                                                                                                            aaGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAg            >  2:497153/1‑149 (MQ=255)                                                                                                                                              ggCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTgg         >  2:494391/1‑150 (MQ=255)                                                                                                                                              ggCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa        <  2:565956/151‑1 (MQ=255)                                                                                                                                               gCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCTTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa        >  2:621656/1‑150 (MQ=255)                                                                                                                                               gCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa        >  2:384627/1‑150 (MQ=255)                                                                                                                                                cTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAAc      >  2:728793/1‑151 (MQ=255)                                                                                                                                                    tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTg  >  2:756820/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    AAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCA‑TGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCT‑GACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTG  >  NC_000913/345180‑345475

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