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

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027                                                                                                                                                           |                                                                                                                                                     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTAttttt                                                                                                                                                      <  2:306167/151‑1 (MQ=255)   tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑cc                                                                                                                                                   >  1:663205/1‑151 (MQ=255)     gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑ccgc                                                                                                                                                 <  1:1060771/151‑1 (MQ=255)     gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑ccgc                                                                                                                                                 <  2:1089116/151‑1 (MQ=255)       aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑ccgccg                                                                                                                                               >  1:1808632/1‑151 (MQ=255)        ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑Ccgccgc                                                                                                                                              >  1:1434841/1‑151 (MQ=255)         gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGca                                                                                                                                             <  1:1575829/151‑1 (MQ=255)         gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGca                                                                                                                                             <  2:1604390/151‑1 (MQ=255)          ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGcaa                                                                                                                                            <  1:1075133/151‑1 (MQ=255)               tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACaaa                                                                                                                                        >  1:365355/1‑150 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACaaaa                                                                                                                                       <  2:235129/150‑1 (MQ=255)                 aGCGTCGCCCTTCAGGCACTGCCACCCGAAGACCTGTGCCAGGGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCATCGCCATGTGCAAAAAT‑GCTTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAaaaaaa                                                                                                                                       <  1:387973/149‑1 (MQ=255)                  gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAAccg                                                                                                                                    >  1:912497/1‑151 (MQ=255)                   cttgccgcCTCAGGCATTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGcaaca                                                                                                                                          <  1:295295/139‑1 (MQ=255)                    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                  <  2:441451/151‑1 (MQ=255)                    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                  <  1:1105219/151‑1 (MQ=255)                          cATCAGCCAATTCCAGCCCCAGACCTGTGCCAGCGCCTCCCGTGA‑TTCATTCCCG‑CCAACAACCGCGCATTTATCCAACCCCATGTGCAAAAAT‑GCCTTCGCGGCGGCGGTCTCCC‑AGCTTTTTTT‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                              <  2:242101/149‑1 (MQ=18)                            tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtctctc                                                                                                                          >  2:1153426/1‑151 (MQ=255)                              aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                         >  1:989938/1‑150 (MQ=255)                                    tcccgcccGCAGACCTGTGTCAGGGGCTCCCGTTT‑TTCTTTCCCG‑CCAACACCCGCGCTTTTCCCCACCGCCGTGTGAAAAAAC‑CCCTCGGGGGGGGCTGTCTCCC‑ACTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                   <  2:1134357/144‑1 (MQ=255)                                    tcccAGCGGCACCCCTTTGGCAGGGCCAACCGTAA‑TTCATTCCCC‑CAAACAACGGGTCTTTTCCCCATCCCTTTGTAAAAAAAT‑GCCTTCGGGGCGGCTGTCTGCC‑AGTTTTTTTT‑CCGCCCAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                     <  2:283535/146‑1 (MQ=255)                                    tcccAGCCGCAGACTTGTGTCAGCGGCTCCCGTGA‑TTCTTTCCCC‑CCAACACCCGCGCTTTTATCCAACGCCTTGTGAAAAAAT‑GCCTCCGCGGCGGCTGTCTGCC‑AGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCagta                                                                                                                     <  2:1252093/146‑1 (MQ=255)                                          ccGCAGACTTGTGCCGGCGGCTACC‑TGATTTCTTCCCCC‑CCACCACCCGCACTTTAATCCACCGCCTTGTCAAAAAAT‑CCTTCCGCGGGGCTTGTCTGCC‑GTTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                             <  1:899741/150‑1 (MQ=255)                                           cGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  1:1612955/1‑150 (MQ=255)                                              aGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                           >  2:1696826/1‑147 (MQ=255)                                               gACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  2:601809/1‑146 (MQ=255)                                                 ccTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  2:1615529/1‑144 (MQ=255)                                                    gtgtCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                            >  1:1060482/1‑141 (MQ=255)                                                        cAGCGGCTACCGTGA‑TTCATTCCCG‑CCAACAACCGCGCATTTATCCAACCCCATGTGCAAAAAT‑GCCTTCGGGGCGGCTGTCTGCC‑AGCTTTTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                              <  2:1449651/151‑1 (MQ=255)                                                         aGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaa                                                                                                 >  2:883993/1‑145 (MQ=255)                                                         aGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGaaaaaaaat                                                                                             >  2:1317644/1‑145 (MQ=255)                                                             gccaCCGTGA‑TTC‑TTTCCGCCCAACACCCGCCCTTTTACCCAACACCATTCGCAAAAGC‑GCCTTCGCGGGGGCTGTCTGCC‑ACCTTTTTTC‑CCGCGACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                           <  2:700538/146‑1 (MQ=255)                                                                 ccGTGA‑TTCATTCCCG‑CAAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAAT‑CCCTTCGCGGGGGCTGTCTCCC‑AGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                       <  2:244258/149‑1 (MQ=255)                                                                  cggtatttaaTTCCCC‑CCAACACCCGCGCTTTTACCAAACCCCATGGGAAAAAAT‑GCTTGCGGGGCGGCTGTCGGCC‑AGCTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                     <  2:741391/142‑1 (MQ=255)                                                                   gTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattttatacc                                                                                      >  2:370371/1‑135 (MQ=25)                                                                   gTGA‑TTCATCCCCG‑CCAACACCCGCGCTTTTACCCACCCCCAGTGGAAAAAAT‑GCCTCGGCGGGGGCTGTCGCCC‑ACTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                     <  1:1625161/149‑1 (MQ=255)                                                                        ttCTTTCCCG‑CCAACACCCGCGCTTTTACCCACCGCCTTTTCAAAAAAG‑GCCTTCGCGGCGGCTTTCGCCC‑AGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                         <  2:1289583/141‑1 (MQ=255)                                                                        ttCTTTCCCG‑CCAACACCCGCGCTTTTACCAAACGCCAGGTGAAAAATT‑GCCTTCGCGGCGGCTGTCTGCC‑ACTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                                   <  2:373263/147‑1 (MQ=255)                                                                               ccG‑CCAACACCCGCGTATTTACCCACCGCCATGTGAAAAAAT‑GCCTTCGGGGCGGCTGTCTGCC‑AGCTATTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                         <  2:35757/150‑1 (MQ=255)                                                                                cccccAACAACCGCCAATTAACCCACCGCCGTTTGAAAAAAT‑CCCTTCGGGGCGCCTGTCTGCC‑ACCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                          <  2:1362836/146‑1 (MQ=255)                                                                                     aacaccCGCGCTTTTACCCAACGCCATGTGAAAAAAT‑GCTTCCGCGGCGGCTGTCTCCC‑AGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                    <  2:1504090/150‑1 (MQ=255)                                                                                      acaccCGCGCTTTTACCCACCCCCGGGT‑CAAAAATGGCCTCCGCGGCGGCTGTCGGCC‑AGCTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                    <  2:1300496/149‑1 (MQ=255)                                                                                       cacccGCGCATTTATCCACCCCCATGTGCAAAAAT‑CCCTTCGCGGCGGCTGTCGGCC‑AGCTATTTTT‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                  <  2:904266/147‑1 (MQ=255)                                                                                             cccttttACCCAACGCCTTGTGAAAAAAT‑GCCTTCGGGGCGGCTTCCTGCC‑‑GCTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                           <  2:224402/147‑1 (MQ=255)                                                                                                 tttATCCACCGCCAGGTGAAAAAAT‑CCTTGCGCGGGGGCTTTCTGCC‑GGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                          <  1:1044850/148‑1 (MQ=255)                                                                                                    acccACCGCCATGTGAAAAAAT‑GCTTTCGCGGCGGCTGTCTGCC‑AGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                      <  2:269373/147‑1 (MQ=255)                                                                                                      ccacccccATGTGAAAAAAG‑CCCTTCGCGGGGGCTGTCTGCC‑AGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                     <  2:1158093/142‑1 (MQ=255)                                                                                                      ccAACGCCATGGGAAAAATC‑GCTGCCGCGGCGCCTGT‑TGCCAAGCTTTTTTC‑CCCCCGAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                     <  1:1510393/148‑1 (MQ=255)                                                                                                      ccAACCCCTTGTCAAAAAAT‑CCCGCGGCGGGGCCTGTCTGCC‑AGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                    <  1:612403/149‑1 (MQ=255)                                                                                                                   aaaaaaT‑CCCTCCGCGGCGGCTGTCTGCC‑AGCTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                      <  2:1830968/149‑1 (MQ=255)                                                                                                                       aaT‑GCCTCCGGGGGGCCTGTCGCCC‑AGTTTTTTTC‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                 <  1:267318/151‑1 (MQ=255)                                                                                                                           tccTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACTTTTAGCCCGCCCTGTTTTGTGGGGAATGCCGCCGGGGACAAAGTGggaagggaagggggggg                                >  2:1731767/2‑142 (MQ=255)                                                                                                                                 gcggggGCTGTCTGCC‑AGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                           <  1:969707/148‑1 (MQ=255)                                                                                                                                 gcggcgGCTGTCTGCC‑AGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                         <  2:1013417/150‑1 (MQ=255)                                                                                                                                        cTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGAGCACGAGAAAAAAATTTAAAACCCCCACGTTGTGTGGGGAACGCTGCTGGTGGAAAAGGGGGAAGGGAAGTGCGGCGGATCaacccccga                   >  2:506690/1‑141 (MQ=255)                                                                                                                                               cc‑AGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg             <  2:240159/148‑1 (MQ=255)                                                                                                                                                ccgccttttttc‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                <  2:1109994/138‑1 (MQ=255)                                                                                                                                                ccattttttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa          <  1:1374165/144‑1 (MQ=255)                                                                                                                                                c‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGATTTTAGGCCGGCATTTTTTGTGGGGAACGCTGCCGGGAGCAAAGTGGGAAAGGGAGGGGGGGGGATCAAATCCAGAACCccgccc          >  2:1231087/1‑148 (MQ=255)                                                                                                                                                     ttttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt       <  1:1507442/147‑1 (MQ=255)                                                                                                                                                     tATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACTTTAAAACCCTCCTTTTGTTTGGGGATCGCTGCTGGGGAAAATATGGGAAGGGGAGGGGGGGGGATAAACTCCCGAAACGCGCtcaat       >  2:1781331/1‑145 (MQ=255)                                                                                                                                                         ttt‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:1273264/150‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                     ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027

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: ‑