Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	55.6%	23.9 / 90.6	63	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (19/9);  new base A (0/35);  total (19/44)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-09

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357741‑358037                                                                                                                                                        |                                                                                                                                                        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)                                                                                                                                             tcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc           <  2:934835/148‑1 (MQ=255)                                                                                                                                              cccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg          <  1:427362/149‑1 (MQ=255)                                                                                                                                                ccccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt       <  1:1337306/149‑1 (MQ=255)                                                                                                                                                ccccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt       <  2:451068/149‑1 (MQ=255)                                                                                                                                                    gcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAt   <  1:797247/151‑1 (MQ=255)                                                                                                                                                      aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAA‑TTTTAACCCCCCCTGTTTTGTGGGGAAAGCTGCC‑GGGTACAAAGT‑GGGAAGGGGAGGGGGGGGGAACAACTCCCGAAACCCCCCAATTTAGTTTGCCTCCAta  >  2:1588193/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                        TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG  >  NC_000913/357741‑358037

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