Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	45.8%	46.2 / 33.9	60	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base A (0/27);  total (17/43)
Fisher's exact test for biased strand distribution p-value = 3.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039                                                                                                                                                           |                                                                                                                                                           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                           <  1:71206/150‑1 (MQ=255)    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                        <  4:797434/150‑1 (MQ=255)    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                       <  4:1101845/151‑1 (MQ=255)     aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                      >  2:690695/1‑151 (MQ=255)     aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                      <  1:354709/151‑1 (MQ=255)        gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                                   <  1:473467/151‑1 (MQ=255)                      aatcccgcccGCAGCCTTGGGTCAGCGCCTCCGGTTTTTCATTCCCCCCAACACCCGCGCTTTTCC‑CCACCCCCTTGTAAAAAAAT‑C‑CTTG‑C‑GGGCGGGC‑TCTCCGCCAGTTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCagta                                                                                                                                      <  4:735337/142‑1 (MQ=255)                         cccAGCCGAAGCCTTGTGTCAGCGGCTACCGTGATTCTTTCCCCCCAACACCCGCGCTTTTAC‑CCAACGCCATGTAAAAAAAT‑G‑CTTT‑C‑GGGGCGGT‑TTTCTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                                  <  4:1483118/150‑1 (MQ=16)                          cccgCCGCACCCCTGTTCCAGCGTCTACCGTAATTCATCCCCGCCAACACCCGCCCTTTTCC‑CCACCCCATTTTGAAAAAAT‑G‑CCTC‑C‑GCGGCGCC‑TTTCTGCCATTTATTTTCCCCCCCCAACAAAACCGCCGTTCTCTCCagta                                                                                                                                      <  1:124427/143‑1 (MQ=255)                          ccACCCGCAGCCCTGTGTCACCGCCTACCGTAATTCATTCCCGCCAACAACCGCGCATTTAT‑CCAACGCCAGGTGCAAAAAT‑G‑CTTT‑C‑GCGGCGGC‑TGTCTGCCAGCTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                                 <  2:965151/151‑1 (MQ=21)                                gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                            >  1:368808/1‑149 (MQ=255)                                   gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                            >  4:171573/1‑146 (MQ=255)                                          gtCGGCGGCTCCCGTGTTTTTTTCCCCCAAACACCCGCGCATTTAC‑CAACCCCCATGT‑AAAAAAA‑G‑CCTTCC‑GGGGGGGC‑TTTCTGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                             <  4:1221425/139‑1 (MQ=255)                                           tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑GAGA‑G‑g‑aaaa                                                                                                              >  2:180535/1‑147 (MQ=255)                                             aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgcgaggaaaagagt                                                                                                              >  4:1552096/1‑137 (MQ=255)                                             aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑GGGA‑G‑a‑aaaaaaat                                                                                                          >  1:379985/1‑145 (MQ=255)                                             aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑AAGA‑G‑a‑aaaaaatt                                                                                                          >  4:1309151/1‑145 (MQ=255)                                               cggttcccGTGTTTCTTTCCCGCAAACACCCGCGCTTTTAT‑CCAACGCCATGTGAAAAAAG‑C‑CCTC‑C‑GGGGGGGC‑TTTCGGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAg                                                                                                             <  3:176502/140‑1 (MQ=255)                                                ggcaccggtgtTTCATTCCCGCCAACACCCGCGCTTTAAC‑CCAACCCCATGTGAAAAAAT‑C‑CCTC‑C‑GCGGCGGC‑TGTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTaa                                                                                                      <  4:720700/144‑1 (MQ=255)                                                ggCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTAC‑CCACCGCCATGTGAAAAAAT‑C‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATa                                                                                                              <  4:1630490/144‑1 (MQ=255)                                                ggCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑Aaca‑g                                                                                                                     >  2:351393/1‑139 (MQ=255)                                                   tcccGTGATTCATTCCCCCCAACACCCGCGCATTTAC‑CCACCCCCATGTGAAAAAAT‑G‑CCTT‑C‑GCGGCGGT‑TGTCTGCCGGTTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGc                                                                                                            <  1:107682/141‑1 (MQ=255)                                                   tACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑CAGA‑G‑g‑aaaaaagttta                                                                                                       >  2:1275378/1‑139 (MQ=25)                                                   tACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑Gaggggaaaatatttaa                                                                                                          >  4:940038/1‑134 (MQ=25)                                                        tGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑GAGAGG‑A‑AAAAAT‑TTTAACCCCgcc                                                                                               >  2:928461/1‑149 (MQ=25)                                                               ttccccccAACACCCGCGCTTTTAC‑CCACCCCCTTGTGAAAAAAT‑G‑CCTC‑C‑GCGGCGGC‑TGTCTCCCAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCa                                                                                              <  4:736073/138‑1 (MQ=255)                                                                  ccGCAAACACCCGCGCGTTTAC‑CAACCCCCATGT‑CAAAAATGG‑CCTT‑C‑GCGGCGGT‑TGTCCGCCTGTTTTTTTCCCGCGCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTgg                                                                                    <  3:238826/151‑1 (MQ=255)                                                                      caacaaCCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGGG‑G‑A‑AAAAGG‑TTAAAGcccccatgtttcgtgg                                                                                     >  3:63802/1‑134 (MQ=25)                                                                           cccGCGCTTTTAC‑CCACCCCCTTGTCAAAAAAT‑G‑CCTG‑C‑GGGGCGGC‑TGCCGCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                            <  4:1273795/149‑1 (MQ=255)                                                                           accccccATTTATCCCAACCCCATG‑ACCAAAAT‑G‑CCCT‑CGGCGGCGCC‑TTCCTCCCCGTTTTTTCTCCCCGGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                            <  4:1047779/145‑1 (MQ=255)                                                                            ccGCGCTTTTAC‑CAACCCCCATTTGAAAAATT‑C‑CCTC‑G‑GGGGCGCC‑TGTCTGCCACTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                          <  4:862945/151‑1 (MQ=255)                                                                             cgcgtttttAC‑CAAACGCCGTGTGCAAAAAT‑G‑CCTT‑C‑GGGGCGGC‑TTTCTGCCATCTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCaa                                                                                 <  4:875199/137‑1 (MQ=255)                                                                                   tttAC‑CCAACGCAATGTGAAAAAAT‑G‑CTTC‑C‑GGGGCGGC‑TGTCTGCCAGTTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                                    <  3:163790/150‑1 (MQ=255)                                                                                   tttAC‑CCAACCCCTTGTCAAAAAAG‑G‑CCTC‑C‑GCGGGGGC‑TGTCTGCCATTTTTTTCCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                                       <  3:433355/147‑1 (MQ=255)                                                                                    ttAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑CAGA‑GAA‑AAAAAC‑TTTAAACCC‑CCCTGTTTTGTGGGGATCGCTGCTggggga                                                                   >  4:451410/1‑146 (MQ=25)                                                                                      aT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑GAGG‑G‑A‑AAAAAA‑TTTTAAACCGCCATGTTGTGTGGGGATCGCGGCTGGGAaaaa                                                                 >  4:838358/1‑150 (MQ=25)                                                                                           aacccattgggaaaaatT‑G‑CCTT‑C‑GGGGGGGT‑TTTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                               <  1:333630/136‑1 (MQ=255)                                                                                            ccGCCTTGTG‑AAAAAT‑GCCCTG‑C‑GCGGGGGC‑GTTCTGCAAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                           <  2:1627257/150‑1 (MQ=255)                                                                                               cagggtgAAAAAAT‑G‑CCTT‑C‑GCGGCGGT‑TTTCTGCCATCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                     <  3:37191/134‑1 (MQ=255)                                                                                                   gtgAAAAAAT‑G‑CCTC‑C‑GCGGCGGC‑TGTCTCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                      <  2:974201/149‑1 (MQ=255)                                                                                                   gtcaaaaagT‑C‑CCTT‑C‑GCGGCGGT‑TGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                     <  1:184918/146‑1 (MQ=255)                                                                                                    taaaaaaaT‑C‑CCTG‑C‑GCGGCGGC‑TGTCTCCCATTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                                    <  2:1097316/147‑1 (MQ=255)                                                                                                       aaaaaC‑G‑CCGT‑C‑CCGGGGGCTTGTTTGCCGTTTTTTTTTCCCCGGAAAGAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                                 <  4:1523477/151‑1 (MQ=255)                                                                                                           aT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AAAAAA‑ATTTTACACCCCATGTTGTGTGGGGATCGCTGCTGGTAAAAAAAAGGAAAGGGGAGAgcgggg                                            >  3:36359/1‑149 (MQ=255)                                                                                                                       gcggcgGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑GAGG‑G‑A‑AAAAACATTTAACCCCCCCTCGTT‑TTTGGGGATCGCTGCTGGTGACAAAATGGGAAGGGGAgggggggggggcagc                                     >  3:280658/1‑136 (MQ=255)                                                                                                                        gggggggttGTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                                    <  3:450455/143‑1 (MQ=255)                                                                                                                         ggcggt‑TGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                       <  1:197012/146‑1 (MQ=255)                                                                                                                           ggggttGTCTGCCGGCTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  1:85819/146‑1 (MQ=255)                                                                                                                                tgccTGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  2:74987/143‑1 (MQ=255)                                                                                                                                     gctatttATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                       <  4:1240296/145‑1 (MQ=255)                                                                                                                                       cAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑AAGA‑G‑A‑AAAAAC‑ATTAAAACCCCCCTGTTTTGTGGGAAACGCTGCTGGGAAAAATGTGGGAAGGGGAGGGGGGCGGAACAACTCCCGAAACgcgc                        >  4:1542662/1‑148 (MQ=255)                                                                                                                                        agtttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                    <  4:539399/146‑1 (MQ=255)                                                                                                                                           tttttttGCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg                  <  4:126292/148‑1 (MQ=255)                                                                                                                                           tttttccccgccaCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                     <  2:1247191/140‑1 (MQ=255)                                                                                                                                             ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt                   <  3:186823/147‑1 (MQ=255)                                                                                                                                              tttcCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt                   <  1:276567/146‑1 (MQ=255)                                                                                                                                              tttcCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt               <  4:827153/150‑1 (MQ=255)                                                                                                                                               tccccgccgCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt              <  2:570095/147‑1 (MQ=255)                                                                                                                                                  ccgccgCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt              <  2:1077959/147‑1 (MQ=255)                                                                                                                                                  ccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc          <  2:761940/151‑1 (MQ=255)                                                                                                                                                  cccccccaaaaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg            <  4:1686692/143‑1 (MQ=255)                                                                                                                                                     ccgaAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt         <  4:1217243/149‑1 (MQ=255)                                                                                                                                                           caaAACCGCCGTTCTCTCCAGTAGTGGCGGGGG‑GAGA‑G‑AGAAAAAA‑TTTAATCCC‑CGATGTTGGGTGGGAGAAGCTGCTGGGTAAAATGTGGGAAGGGGAGTGCGGGGGGACAGCCCCCAAACCCCGCCAATTTTGTTCgtccccatttc  >  1:467369/1‑140 (MQ=255)                                                                                                                                                           |                                                                                                                                                           GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAAT‑G‑CCTT‑C‑GCGGCGGC‑TGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGC‑AAGA‑G‑A‑AATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039

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