Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	50.9%	31.7 / 80.9	55	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (16/11);  new base A (0/28);  total (16/39)
Fisher's exact test for biased strand distribution p-value = 4.38e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-05

TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357744‑358034                                                                                                                                                        |                                                                                                                                                 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                                                                                                                                         <  1:1105219/151‑1 (MQ=255)      gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                         <  2:441451/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)                                                                                                                                            ttcccgcccAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt         <  1:169715/147‑1 (MQ=255)                                                                                                                                               ccgcccCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc        <  1:1448310/148‑1 (MQ=255)                                                                                                                                               cccccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc     <  1:922980/148‑1 (MQ=255)                                                                                                                                                  ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa  <  1:819990/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                 TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCG‑CCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTCGCGGCGGCTGTCTGCC‑AGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357744‑358034

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