Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,947,850	0	C	T	40.7%	76.5 / 36.1	59	E40K (GAA→AAA)	hdfR	flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (10/25);  new base T (24/0);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 9.62e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAT  >  NC_000913/3947701‑3948000                                                                                                                                                        |                                                                                                                                                        gTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTcc                                                                                                                                                        <  2:222571/151‑1 (MQ=255)       tCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTCACACCCCGTTGA‑TTTTTCAGtttt                                                                                                                                                  >  2:913369/1‑149 (MQ=255)        cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTc                                                                                                                                                 <  2:627325/151‑1 (MQ=255)           gAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAAGGTCTCACCCAGTTGA‑TTTTTCAATTGTCt                                                                                                                                                >  2:818578/1‑149 (MQ=255)            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGAt                                                                                                                                             <  1:380543/151‑1 (MQ=255)            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGTTTACACCCCGTTGT‑TTTTTCACTTTTCt                                                                                                                                                >  1:448884/1‑148 (MQ=255)             ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGAtt                                                                                                                                            <  1:343433/151‑1 (MQ=255)                  gcgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTTACACCCCGTTGG‑TTTTTCAGTTGTTTGATTTGaa                                                                                                                                        >  1:128124/1‑150 (MQ=255)                     ccACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCt                                                                                                                                    <  1:906108/151‑1 (MQ=255)                     ccACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAAGGTCACACCCCGTTGA‑TTTTTCAGTTGTCTGATTCGaa                                                                                                                                        >  2:37408/1‑147 (MQ=255)                          tcctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAAGTTTACAACCCGTTGA‑TTTTTCAGTTGTCTGATTCGAAAGATCAct                                                                                                                                >  2:736840/1‑150 (MQ=255)                           cctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGGTCACACCCCGTTGG‑TTTTTCCGTTTTTTGAttttaaaaatcac                                                                                                                                 >  1:504171/1‑137 (MQ=255)                            ctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGGAAAAGGTCACACCCCGGTGG‑TTTTTCAGTTGTCTTGTTTGAAAACGCActtc                                                                                                                              >  1:542513/1‑148 (MQ=255)                               aCGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGAc                                                                                                                         <  1:623937/151‑1 (MQ=255)                                cGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGGAAAAGGTCACACCCCGTTGA‑TTTTTCAGTTGTTTGATTTGAAAGCTCACTG‑CTGAc                                                                                                                         >  1:236817/1‑150 (MQ=255)                                      gcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTTGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCTCTG‑CTGACTGGGTc                                                                                                                   <  2:33153/150‑1 (MQ=255)                                         tGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTTACAACCCCTTGA‑ATTTTCCGGTGTTTGAtttcaaaagccac                                                                                                                                 >  1:223055/1‑123 (MQ=255)                                            cAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATaga                                                                                                            <  2:739758/151‑1 (MQ=255)                                                tGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAACTCACTG‑CTGACTGGGGCAGATAGAGCg                                                                                                         >  2:709594/1‑150 (MQ=255)                                                       gAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACGCCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGcc                                                                                                 <  2:33668/151‑1 (MQ=255)                                                       gAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGc                                                                                                  <  1:11590/150‑1 (MQ=255)                                                           gTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAAGGTTACACCCCGTTTG‑TTTTTCCTTTTTTTGATTTGAAAACTCACTG‑CTGACTGGGGCAGGTAGAGCGATTTAGcccccc                                                                                             >  1:538591/1‑147 (MQ=255)                                                           gTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTCACACCCAGTTGA‑TTTTTCAGTTTTTTGATTTGAAAAATCACTG‑CTGTCTGGGTCAGATAGAGCGATTCAAcc                                                                                                 >  2:73267/1‑147 (MQ=255)                                                            tttCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGGGAAAAGGTT‑ACACCCACTTGA‑TTTTTCAGTTGTGTCATTTGAAAGCTTTCTG‑TTTACTGGGTCAGATAGAGCGGTTTAGGCGCCCg                                                                                            >  2:434780/1‑151 (MQ=255)                                                              tCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGTAAAAGGTCACACCCCATTGA‑TTTTTCCCTTGTCTGATTTTAAAGATCACTG‑CTGACTGGGGCAGATAGAGCGATTCAGCCGCTc                                                                                             >  1:770745/1‑148 (MQ=255)                                                              tCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGa                                                                                           <  1:62071/150‑1 (MQ=255)                                                                 gCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCTCTCGACCa                                                                                        <  1:335127/150‑1 (MQ=255)                                                                 gCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGAc                                                                                          <  2:146185/148‑1 (MQ=255)                                                                  cATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGTGAAGATTCACACCCCAATGAGTTTTCCA‑TTTTGTGGTTTGAAAGCTCACTA‑CAGGGTGGGTCGGATAAAGCGATTTAGCCGCTCGGCCa                                                                                        >  1:656731/1‑149 (MQ=255)                                                                  cATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGCGGAAAAGTTCACACCCCGCTGA‑TTTTTCCGTTGTCTTATTTGGAAGCTCCCCGCCGGTCTGGG‑CAGAAAGAGCGGTTCAGGCGCTCGAcc                                                                                         >  2:477767/1‑148 (MQ=255)                                                                        gCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGGAAAAGGTCACAACCCGGTGA‑TTTTTCACTTTTTTGATTCGAAAGGTAACTG‑CCGACTGGGGTAGATAGAGCGGTTTAGCCGCCCGACCaaaaagac                                                                                 >  1:877542/1‑145 (MQ=255)                                                                         cAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGGAAAGGTTCACACCCCGGTGG‑TTTTTCAGTTGTTTTATTTGGAAAACAACCG‑CCGGCTGGGTGAGATAGAGCGATTGAACCGCGCGGACAAAgaga                                                                                  >  1:417266/1‑148 (MQ=255)                                                                                 tttCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGc                                                                        <  2:554159/150‑1 (MQ=255)                                                                                      ccGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACtt                                                                  <  2:378932/151‑1 (MQ=255)                                                                                      ccGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTCACACCCCGGTGA‑TTTTTCAGTTGTCTGATTCTAAAACTCCCTG‑CTGACTGGGGCCGATAGAGGGGTTCAGCCGCCCGACCCAAATGACGCGGTTGGCTAACtt                                                                  >  2:951033/1‑151 (MQ=255)                                                                                          cAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGGAAAAGGTCACACCCAGATGA‑TTTTTCCGTTGTTTTTTTTGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGGGGTTCAGCCGCCCGACCAAAAGGACGCGGTCGGGTAACtttc                                                                >  1:29584/1‑147 (MQ=255)                                                                                          cAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAg                                                              <  1:59568/151‑1 (MQ=255)                                                                                              ggTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAGGTTCACACCCAGTTGA‑TTTTTCCGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCCGATAGAGCGATTCAGCCGCTCGACCCAAGTGACGCGTTCGGCTAACTTCCAGGaa                                                           >  2:731085/1‑150 (MQ=255)                                                                                              ggTTAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAGGGTCACCACCACATGG‑TTTTTTATTTTTTTGGTTCGAAAACCCACTG‑CTGACTGGGGGAGGTAGAGCGATTTCGCCGCCCGCCCCAAAAGACGCGGTCGGCTAAATTCCCGGaa                                                           >  1:965900/1‑150 (MQ=255)                                                                                                ttAAACGGATATTGTTTCTGTGGA‑GGAGGGAAAAATTCAGACCCTGATGA‑AGTTTCTCTTTTCTGATGCGAAAACGCACTG‑CTGACTGGGTGAGGTTGAGAGAGTCCGGCGCCCGACAAAAGACAGGCGGTCGGCCCAATTCCCGGaa                                                           >  1:108637/1‑148 (MQ=255)                                                                                                 tAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGt                                                        <  1:928174/150‑1 (MQ=255)                                                                                                 tAAACGGATATTGTTTCTGTGGC‑GGGGGAAAAAGGTCACCCCCCGTTGG‑TTTTTCCGTTTTTTGTTTCGGAAACTCACTG‑CTGACTGGGGCAGATAGAGCGATTTAGCCGCCCGACCAAAATGACGCGTTTGGCTAACTTCCAGGaaa                                                          >  1:930394/1‑148 (MQ=255)                                                                                                  aaaCGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGttt                                                      <  1:399503/151‑1 (MQ=255)                                                                                                   aaCGGATATTGTTTCTGTGGC‑GGGTGAAAAAGGTTTCACCCCGTTGG‑TTTTTCAGTTTTCTCTTTTTAAAGCACACTG‑CTGACTGGGGCAGAAAGAGCGATTCAGACGCTCGACCCAAGTGGCGCGTTCGGCTAACTTCCAAGaaaagttt                                                     >  2:774962/1‑146 (MQ=255)                                                                                                       gATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAAc                                                 <  1:852439/151‑1 (MQ=255)                                                                                                          atTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAAt                                              <  1:426432/151‑1 (MQ=255)                                                                                                                    tgGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATcc                                     >  2:238343/1‑150 (MQ=255)                                                                                                                     ggC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCcac                                   >  1:763889/1‑151 (MQ=255)                                                                                                                      gC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCcaca                                  <  1:731058/151‑1 (MQ=255)                                                                                                                                aaGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTcc                         <  2:86560/151‑1 (MQ=255)                                                                                                                                     tCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCaa                     >  2:380448/1‑150 (MQ=255)                                                                                                                                     tCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCaa                     >  2:797476/1‑150 (MQ=255)                                                                                                                                      cacaCCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCaaa                    >  2:157644/1‑150 (MQ=255)                                                                                                                                      cacaCCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCaaa                    <  1:73267/150‑1 (MQ=255)                                                                                                                                        cacCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAatt                 >  2:399567/1‑151 (MQ=255)                                                                                                                                            cAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTAttt             <  1:268418/151‑1 (MQ=255)                                                                                                                                             aGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTg            <  2:224945/151‑1 (MQ=255)                                                                                                                                                tGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCg         >  1:562478/1‑151 (MQ=255)                                                                                                                                                  a‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGGCGt        <  2:519762/150‑1 (MQ=255)                                                                                                                                                      ttCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATg    >  1:492763/1‑151 (MQ=255)                                                                                                                                                        ccAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAt  >  1:525088/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        GTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGAAAAGGTTCACACCCAGTTGA‑TTTTCCAGTTGTCTGATTCGAAAGCTCACTG‑CTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAT  >  NC_000913/3947701‑3948000

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