Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	44.6%	79.9 / 21.9	65	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (23/13);  new base C (0/29);  total (23/42)
Fisher's exact test for biased strand distribution p-value = 1.04e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCT  >  NC_000913/357730‑358018                                                                                                                                                         |                                                                                                                                               aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                >  2:358006/1‑151 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                               <  2:366583/149‑1 (MQ=255)     gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                            >  2:154941/1‑151 (MQ=255)      tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                           >  1:1511654/1‑151 (MQ=255)       aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                          <  1:1172872/151‑1 (MQ=255)         gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                         >  1:295409/1‑150 (MQ=255)         gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                        <  2:318939/151‑1 (MQ=255)         gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                        <  2:155836/151‑1 (MQ=255)            tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                      >  2:121487/1‑150 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                   <  1:651382/151‑1 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                  <  2:988260/150‑1 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                  <  2:274497/150‑1 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                 <  2:938242/151‑1 (MQ=255)                 aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                <  1:54201/151‑1 (MQ=255)                   cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                              <  2:1076950/151‑1 (MQ=255)                      aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                           >  2:714464/1‑151 (MQ=255)                        cgcaccggtcattcccACCCGCAGACCTGTGTCAGCGGCTACCGTGATTCTTTCCCGCCAACACCCGCGCATTTATCCACCGCCATGTGAAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGCAAAAAAACCGCCGTtc                                                                                                                         <  1:494094/137‑1 (MQ=18)                                gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  1:398266/1‑151 (MQ=255)                                    cccccgCCGCAGTCCGTTGCCGGCGGCTCCCGTGTTTCTTTCCCCCAAACACCCGCGCTTTTCCCAACCGCCATGTCAAAAAAGCCCTC‑C‑GCGGCGGCTTTCTCCCATTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                              <  1:1109668/145‑1 (MQ=255)                                      cccgCCGCAGACCTGTGTCGGCGGCTACCGTTTTTCTTTCCCCCCAACACCCGCGCTTTCACCCACCGCCATGTGAAAAAACGCCTC‑C‑GGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                           <  1:1553474/148‑1 (MQ=255)                                        aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                         >  2:1318137/1‑151 (MQ=255)                                          ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                       >  2:338129/1‑151 (MQ=255)                                           cacggACTTGTGTGGGGGGCTCCCGTTTTTCTTTCCCCCCACCACCCGCGCATTTACCCACCGCCATGTAAAAAAACGCTTT‑C‑GGGGCGGTTGTCGGCCAGTTTTTTTTCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                      <  1:498559/147‑1 (MQ=255)                                             cAGCCCTGTGTCGGCGGCTCCCGTGATTCATTCCCCCAAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCCTT‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                    <  2:998994/151‑1 (MQ=255)                                             cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGTCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggc                                                                                                     >  2:209387/1‑148 (MQ=255)                                             cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                     >  2:746458/1‑148 (MQ=255)                                              agccttgttTCAGCGGCTACCGTGATTCTTTCCCCCAAACACCCGCGCTTTTACCCAACCCCATGTGAAAAAATGCCTT‑G‑GCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                   <  2:875453/146‑1 (MQ=255)                                              aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggg                                                                                                   >  1:426398/1‑147 (MQ=255)                                                cccTGTGGCAGGGCCTACGGTGATTCATTCCCGCCAAAAACCGCCATTTTACCCAACGCAAGGGAAAAAAAGGCCTG‑G‑GGGGGGGGGGTGGGCCAGTTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                 <  2:108088/150‑1 (MQ=255)                                                    gtgtagggggCTCCCGGGTTTCATTCCCGCCAACACCCGCCTTTTCCCCCACCGCCATGGAAAAAAATCCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt                                                                                            <  2:534562/143‑1 (MQ=255)                                                    gtgtCGGCGGCTACCGTTTTTCTTTCCCCCAAACACCCGCGCTTTTACCAACCCCCATGTGAAAAAATGCTTC‑C‑GGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                             <  2:635760/150‑1 (MQ=255)                                                    gtgcCAGCGCCTCCGGTGTTTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCAATGTGAAAAAATCCCTG‑C‑GCGGCGGCGGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt                                                                                            <  1:463897/151‑1 (MQ=255)                                                    gtgcCAGCGCCTACCGTTTTTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCATGGTGAAAAAATCCCTC‑C‑GCGGCGGCTGTCGGCCAGTTATTTGCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                             <  1:187299/150‑1 (MQ=255)                                                     tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaaa                                                                                            >  2:928907/1‑149 (MQ=255)                                                         ggcgcCTCCCGTGTTTCTTTCCCGCAAACACCCGCGCTTTTACCCAACCCCATGTGAAAAAATGCTTT‑C‑GCGGCGGCTGTCTCCCACCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑tt                                                                                      <  2:290296/150‑1 (MQ=255)                                                           gggCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCTTC‑C‑GGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTa                                                                                    <  2:76224/150‑1 (MQ=255)                                                            ggCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑aaaaaaaattt                                                                                     >  2:1276094/1‑142 (MQ=255)                                                             cctcccgtttttCATTCCCCCCAACACCCGCCCTTTTACCCACCGCCATG‑GAAAAAAAGCCCT‑CGGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAg                                                                                  <  1:1127150/141‑1 (MQ=255)                                                              cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑aaaaaaaatttaa                                                                                   >  2:1319992/1‑140 (MQ=25)                                                               tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaaaaaagttaaaac                                                                                 >  2:511171/1‑139 (MQ=25)                                                               tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaacatttaaac                                                                                 >  2:1337555/1‑139 (MQ=25)                                                                 ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAGG‑TTTAAGccg                                                                               >  2:405848/1‑150 (MQ=25)                                                                              ccGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑ATTAAAACCGCCCTGTTTTGTgg                                                                 >  1:1530347/1‑151 (MQ=25)                                                                               cccaaccacccGCGCTTTTACCCACCCCAGTGTGAAAAATTCCTTT‑C‑GCGGGGGTGTTCCGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGc                                                                <  2:1357122/142‑1 (MQ=255)                                                                                gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAC‑TTTAAAACCCCCCTGTTGTGTgggg                                                               >  2:485882/1‑149 (MQ=25)                                                                                ccaaacaccCGCGTTTTTACCCACCGCCATGTAAAAAAATCCCTG‑C‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCa                                                               <  2:655805/148‑1 (MQ=255)                                                                                     cacccGCGCTTTTACCCACCCCCATGGAAAAAAATGCCTG‑C‑GGGGCGGTTGTGTCCCAGGTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                          <  1:693755/148‑1 (MQ=255)                                                                                       cccGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTC‑C‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                        <  2:276949/150‑1 (MQ=255)                                                                                         cccgcTTTTACCCAACCCGATGTG‑AAAAATTCCTTCC‑GCGGGGCCTGCCTCCCAATTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                      <  2:767523/149‑1 (MQ=255)                                                                                          gcgcTTTTACCCAACCCCATGTAAAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGctgct‑g                                                    <  2:180363/151‑1 (MQ=255)                                                                                           cgcATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑AAAAAAA‑TTTAAACCCCCCCTGTTTTGTTGGGATGGCTCCT‑gg                                                   >  2:1321993/1‑151 (MQ=25)                                                                                                 taccaacccccATGTGAAAAATGGCCTC‑C‑GCGGCGGCTGTCTCCCAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACa                                              <  2:1060254/141‑1 (MQ=255)                                                                                                    ccAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAACTTTTAACCCCCCCAT‑TTTTGTGGGGAACGCTGCT‑GGGGACAAagt                                          >  1:411651/1‑151 (MQ=21)                                                                                                            atgtaaaaaaaTGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAggg                                 <  2:595460/145‑1 (MQ=255)                                                                                                              gtgAAAAAATCCTTC‑C‑GGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGaa                               <  2:1323473/151‑1 (MQ=255)                                                                                                                 aaaaaaTCCCTC‑C‑GGGGGGGCTGCCTGCCGTTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTg                            <  1:217770/150‑1 (MQ=255)                                                                                                                    aaaTGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTATCCCCCCCTGTTTTGGGGGGATCGCTGCTGGGTAAAAAGAT‑GGAAGGGGagggggg                          >  2:505945/1‑145 (MQ=255)                                                                                                                     aaTGCCTC‑C‑GCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTgcggc                        <  2:295425/151‑1 (MQ=255)                                                                                                                      aTGCCTC‑C‑GCGGCGGTTGTCTCCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTgcggcg                       <  1:447723/151‑1 (MQ=255)                                                                                                                         ccTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑AAAAAAA‑ATTAAAACCCTCCTGTTTTGTGGGGAATGCTGCT‑GGTGAAAAAGTGGGAAAGGGaggggggggg                       >  2:1267298/1‑140 (MQ=255)                                                                                                                              ccggggcggcTTTCGGCCGCCTTTTTTCCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATc                   <  2:348460/145‑1 (MQ=255)                                                                                                                                 gggcggcTGTCTGCCAGCTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTc              <  2:154407/150‑1 (MQ=255)                                                                                                                                     ggcTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGa          <  2:185732/151‑1 (MQ=255)                                                                                                                                     ggcTGTCTGCCAGTTTTTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGa          <  1:275165/151‑1 (MQ=255)                                                                                                                                          cctcccggtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg      <  1:813443/139‑1 (MQ=255)                                                                                                                                             gCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg    <  2:120736/149‑1 (MQ=255)                                                                                                                                             cccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt  <  2:1240830/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                               ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCT  >  NC_000913/357730‑358018

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