Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	46.9%	74.3 / 18.0	64	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (24/10);  new base C (0/30);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 8.02e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357733‑358024                                                                                                                                                     |                                                                                                                                                  cGTTGGCCTGATAAGCGTAGCGCATCAGACAATTACAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATGTCCGCCAACAACCGCTCATTTACGCAACGGCATGTGCAAAAATGCCTTC‑GCGGCGGCTATCTGCCAGCTATGTTTc                                                                                                                                                  <  1:939385/149‑1 (MQ=14)  gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                               >  2:1450369/1‑151 (MQ=255)    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                             >  1:599413/1‑151 (MQ=255)      gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                           <  2:859446/151‑1 (MQ=255)      gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                           <  1:1314138/151‑1 (MQ=255)         tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                         >  1:506113/1‑150 (MQ=255)         tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                         >  1:921270/1‑150 (MQ=255)             aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                    <  1:1359615/151‑1 (MQ=255)                  tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                               >  2:935314/1‑151 (MQ=255)                   aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                              >  2:1326/1‑151 (MQ=255)                          cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                       >  2:523498/1‑151 (MQ=255)                             gaattcccACCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTGCAAAAATGCCTCC‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                    <  1:400853/145‑1 (MQ=18)                             caaaTCCCACCCGCAGACTTTTGTCAGCGGCTCCCGTGATTCTTTCCCCCCAACACCCGCGCATTTATCCACCGCCATGTCAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagt                                                                                                                    <  1:845604/149‑1 (MQ=255)                                aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                 >  1:869322/1‑151 (MQ=255)                                 ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCTGTAGTg                                                                                                                >  1:215660/1‑151 (MQ=255)                                  cccAGCCGCAGCCCTGTGCCAGCGGCTACCGTGATTAATTCCCGCCAACACCCGCGCATTTACCCAACGCCTTGTAAAAAAATGCCTCC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                               <  1:205169/150‑1 (MQ=255)                                   ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                              >  2:1276439/1‑151 (MQ=255)                                        cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                         >  2:438042/1‑150 (MQ=255)                                         gCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTCCCCCAAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                        <  2:497215/151‑1 (MQ=21)                                         gCAGACTTGTGTCAGCGCCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCACCGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                        <  2:306067/151‑1 (MQ=21)                                           aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggac                                                                                                       >  2:190467/1‑147 (MQ=255)                                            gcccTGTGTCAGCGGCTACCGTGATTCATTTCCGACAACAACCGCGCATTTCCCCAACGCCATGTACAAAAATGCCTTC‑GGGGCGGCTGTCTCCCGTTTATTTTTCCGCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                     <  2:1096285/149‑1 (MQ=21)                                             cccTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTTC‑GGGGGGGCTGTCTGCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                    <  1:721587/150‑1 (MQ=255)                                              ccTGTGTCAGCGACTCCCGTGATTCATTTCCGCCAACACCCGCCTATTTACCCACCGCCATGTGAAAAAATCCCTGC‑GCGGCGCCTGTCTGCCAGTTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                  <  2:193390/151‑1 (MQ=255)                                                 gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTGAGAG‑aaa                                                                                                >  1:175550/1‑150 (MQ=255)                                                 gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaa                                                                                                >  2:1262651/1‑150 (MQ=255)                                                   gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAa‑aaaaa                                                                                              >  2:582414/1‑148 (MQ=255)                                                    ccggCGGCTCCGGTGTTTCATTTCCCCCAACACCCGCGCATTTCCCCAACCCCATTTGAAAAAAGCCCCCG‑GCGGGGGCTCTCTCCCGCTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGc                                                                                            <  1:260247/148‑1 (MQ=255)                                                        cGGCTACCGGGTTTATTTCCCCCCAACACCCGCGCATTTACCAACCGCCATGTCAAAAAATGCCTTC‑GCGGCGCTTGTCGGCAAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCtttn                                                                                        <  1:313007/151‑2 (MQ=255)                                                          gCTCCCGTGATTCTTTTCCCCCAACACCCGCGCATTTACCCAACCCCATGGGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg                                                                                      <  2:326654/151‑1 (MQ=255)                                                            tccggTGATCCATTCCCGCAAACACCCGCGCTTTTCCCCAACG‑CATGGTGAAAAAAGCCTTCGGGGGCGGTGGCTTGCCAGCTATTTTTCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg                                                                                      <  1:1273520/145‑1 (MQ=255)                                                             aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaaaaaacattaaacc                                                                                   >  1:731275/1‑138 (MQ=25)                                                             aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaattttaaacc                                                                                   >  2:108562/1‑138 (MQ=25)                                                              cggTGATTCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTCC‑GCGGCGGCTGTCGCCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg                                                                                  <  1:131205/149‑1 (MQ=255)                                                              ccGTGATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTCCAAAAATGCTTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg                                                                                  <  2:755731/151‑1 (MQ=255)                                                                gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAAAAAATTTaaacgcggc                                                                                >  1:280143/1‑144 (MQ=25)                                                                  gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAACTTTTAAccccccat                                                                              >  2:538333/1‑146 (MQ=25)                                                                   attcattcCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  2:501814/151‑1 (MQ=255)                                                                   attcattcCCCCCAACACCCGCGCTTTTATCCAACCCCATGTGAAAAAATGCCTCC‑GGGGCGGCTGTCGGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  2:1322341/151‑1 (MQ=255)                                                                    ttCATTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAATGCCTTC‑GCGGCGGCTGTCTGCAAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  1:762902/150‑1 (MQ=255)                                                                            cGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAAGGCCTTC‑GGGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                    <  2:758031/151‑1 (MQ=255)                                                                                 ccacccGCGCATTTCCCCACCCCGATGCGAAAAAATCCCTCC‑GCGGCGCTTCTCTGCCAGCTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                               <  1:396006/147‑1 (MQ=255)                                                                                 acaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAG‑AAAAAGAGTTATACCCGCCATGTTTGGTGGGAATcc                                                               >  1:186325/1‑150 (MQ=25)                                                                                     ccGCTCTTTTATCCAACGCCAGGTGAAAAAATGCCTCC‑GCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                           <  1:779763/151‑1 (MQ=255)                                                                                             ttATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCCG‑CCTGTTTTGTGGGAAACGCCGCTgggtac                                                    >  2:705248/1‑146 (MQ=25)                                                                                               aTCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAACTTTAAACCCCTCCTGTTTTGTGGGGATCGCTGCTGGTAaaaa                                                  >  2:1225082/1‑150 (MQ=37)                                                                                                acccacccccgtatgaaaaaaGCCTTCGGCGGCGGCTGTCGGCCGGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                 <  1:1195386/136‑1 (MQ=255)                                                                                                   aaCGCCATGTGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                             <  2:128052/151‑1 (MQ=255)                                                                                                    ccGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                            <  2:516676/150‑1 (MQ=255)                                                                                                         aTGTGAAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                       <  2:512904/151‑1 (MQ=255)                                                                                                          tgtgCAAAAATGCCTGC‑GCGGCGGCTGTCTGCCGTCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                        <  2:255303/149‑1 (MQ=255)                                                                                                               aaaaaTGCCTTC‑GCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                 <  2:134841/151‑1 (MQ=255)                                                                                                                aaaaTGCCTTC‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                <  2:506793/151‑1 (MQ=255)                                                                                                                 aacggCCTTC‑GGGGCGGCTGTTGGCCAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                               <  2:281469/147‑1 (MQ=255)                                                                                                                     gCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑GAAAAAATTTTAACCCGCCCTTTTTTGTTGGGAACGCTGCGGGGGAAAAAGTGGGAAGGGGAgggggggggg                           >  1:1272346/1‑142 (MQ=255)                                                                                                                        ttC‑GGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                          <  1:864728/149‑1 (MQ=255)                                                                                                                          ccggggcggcTGTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                         <  1:119666/145‑1 (MQ=255)                                                                                                                                cggcTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                 <  2:599407/151‑1 (MQ=255)                                                                                                                                cggcTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                 <  2:295787/151‑1 (MQ=255)                                                                                                                                  gcTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa               <  1:183293/151‑1 (MQ=255)                                                                                                                                    tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAGATTTAAGCCCCCCCTGTTTTGTTGGGATCGCTGCTGGGGAAAATGTGGGAAGGGGAGGGGGGGGGAGCAGATCCAGaaac             >  2:861969/1‑149 (MQ=255)                                                                                                                                     gtcggcaggttATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg            <  1:1276440/141‑1 (MQ=255)                                                                                                                                     gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  2:588428/150‑1 (MQ=255)                                                                                                                                            agttATTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt     <  1:1262652/148‑1 (MQ=255)                                                                                                                                                tttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  2:921274/149‑1 (MQ=255)                                                                                                                                                tttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  1:935310/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357733‑358024

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