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

GTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAT  >  NC_000913/3947701‑3948000                                                                                                                                                         |                                                                                                                                                       gTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTcc                                                                                                                                                       <  1:943599/151‑1 (MQ=255) gTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTcc                                                                                                                                                       <  2:530158/151‑1 (MQ=255)   tATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCa                                                                                                                                                      <  2:1091514/150‑1 (MQ=255)    aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAg                                                                                                                                                     <  1:925782/150‑1 (MQ=255)        cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTc                                                                                                                                                <  1:1095714/151‑1 (MQ=255)           gAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTg                                                                                                                                              <  2:388084/150‑1 (MQ=255)            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑G‑AAAAGGTT‑ACACCCAGTTGATTTTTCCGtt                                                                                                                                                   >  1:1195000/1‑144 (MQ=255)            aGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAAGTTCACACCCCGTTGATTTTCCAGTTTTCt                                                                                                                                               >  2:656439/1‑148 (MQ=255)               tATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGT‑A‑AAAAGGTTC‑CCCCCCGCTGGTTTTTCAGtttttttttt                                                                                                                                           >  1:579901/1‑141 (MQ=255)                 tGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTCCAACCCCCGTGGTTTTTCCattgtgttgtt                                                                                                                                            >  1:1262102/1‑137 (MQ=255)                  gcgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTTACACCCCGTTGATTTTTCAGTTGTCTGATTTGaa                                                                                                                                       >  2:733049/1‑150 (MQ=255)                     ccACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGc                                                                                                                                    <  1:769232/150‑1 (MQ=255)                      cACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTTACACCCCGTTGATTTTCCAGTTGTTTGATTCGaaa                                                                                                                                      >  2:1205328/1‑147 (MQ=255)                        cctcctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCAc                                                                                                                                <  2:384706/151‑1 (MQ=255)                           cctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCActgc                                                                                                                             <  2:1095027/151‑1 (MQ=255)                           cctTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCAct                                                                                                                               <  1:955591/149‑1 (MQ=255)                                cGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTTACACCCCGGTGATTTTTCAGTTTTCTGATTTGAAAAATCAc                                                                                                                                >  2:971952/1‑143 (MQ=255)                                   gcggcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTggg                                                                                                                     <  1:1136779/151‑1 (MQ=255)                                          gCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTTGGTCAGAt                                                                                                               <  1:830194/150‑1 (MQ=255)                                            cAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATaga                                                                                                            <  1:1118862/151‑1 (MQ=255)                                            cAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATaga                                                                                                            <  1:964261/151‑1 (MQ=255)                                            cAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGA‑G‑GAAAAGTTTACACCCCGATGATTTTTCAGATGTGTGTTTCGAAAGCTCACTGCTGCCTTGGGCTGAt                                                                                                               >  1:720906/1‑148 (MQ=255)                                                 gCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGAt                                                                                                       <  1:526281/151‑1 (MQ=255)                                                  cTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑GAAAAGGTCACACCCCGGTGGTTTTTCAGCTGTCTGCATCCGAAGCTCACTGCTGTCCGGGTGAGATAGAGCGGtt                                                                                                      >  2:194025/1‑151 (MQ=255)                                                   tCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑‑GAAGGGTCACACCCCGTTGATTTTTCAGTTTTCTGTATTGAAAGCTCACTGCTTACAGGGGTCGATAAAGCg                                                                                                         >  2:368428/1‑146 (MQ=255)                                                    cATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGAtt                                                                                                      <  2:817299/149‑1 (MQ=255)                                                          cGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAAGGTCACACCCCGGTGGTTTTTCCGTTTTTTGATTTGAAAAATCACTGCTGACTGGGGCAGATAAAGCGAtttatcccc                                                                                               >  1:66784/1‑143 (MQ=255)                                                           gTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAGGTTTACACCCCGGTTGTTTTTCAGTTGTCTGATTAGAAAACTCACTTTTGGCTGGGGCAGGAAAAGCGATTCa                                                                                                    >  2:699325/1‑144 (MQ=255)                                                                 gCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCaa                                                                                       <  1:752081/151‑1 (MQ=255)                                                                  cATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCaaa                                                                                      <  2:570973/151‑1 (MQ=255)                                                                   aTAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑G‑AAAAGGTT‑AAACCCAGTTGATTTTTCCGTTGTTTGATTTGAAAACTCACTGCTGACTGGGGTAGATAAAGGGATTTAGCCGCCCGCCCaaa                                                                                      >  2:757643/1‑150 (MQ=255)                                                                   aTAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAAGGTTAAACCCCGTTGATTTTTTAGTTGTCTGATTTGAAAAATCACTGCTGACTGGGTGCGATAGAGCGATTTAGCCCCCCGACCaaa                                                                                      >  1:973508/1‑150 (MQ=255)                                                                     aaGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAAGTTCACACCCAGTTGATTTTTCCGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGGCAGATAGAGCGATTCAGCCGCTCGACCAAAGTg                                                                                   >  2:865276/1‑151 (MQ=255)                                                                      aGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGG‑AAAAGGTT‑ACACCCCCGTGATTTTTTCGTTTTCTGATTCGAAAGATCCCTCCTTTCTGGGGTCGATAGAGGGATTTAGCCCCTCGACCcaaa                                                                                     >  1:897620/1‑147 (MQ=255)                                                                        gCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGGTCACACCCCGGTGATTTTTCAGTTGTCTGATTTGAAAGCTCACCGCTGACTGGGGCAGATAGAGCGATTTAGCCCCTCGACCCAAgaga                                                                                  >  2:1083493/1‑149 (MQ=255)                                                                        gCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTTACACCCCGTTGATTTTTTTGTTTTTTGTTTTGAAAGATCACTGCTGACTGGGGCAGATAAAGCGATTTAGccccccgccccaa                                                                                      >  1:605894/1‑134 (MQ=255)                                                                             agtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGGTCACACCCAGTTTATTTTTCAGTTGTCTGATTCGCAAACTCCCTGCTGCCTGGGTCAGATAGAGCGATTCAGCCGCTCCACCCAAAAGAcggggt                                                                            >  2:838586/1‑146 (MQ=255)                                                                               tttttCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑GAAAAGGTCACCCCCCGGTGGTTTTTCAGTTGTCTGGTTTGAAAGGACACTGCTGGCTGGGGGAGGAAAAGGGGTTTAGCCGCTCGGCCaaaaggaagcggttgg                                                                         >  2:779136/1‑138 (MQ=255)                                                                                ttttCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTCACACCCAGTTGGTTTTTCCGTTGTCTGGTTTGAAAACTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCgg                                                                         >  2:141872/1‑150 (MQ=255)                                                                                  ttCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑GAAAAGGTCACACCCACTTGGTTTTTCCGTTGTCTTTTTCTGAAACTCACTGCTGACTGGGGCTGATAGGGCGGTTCAACCCGTCGCCCAAAGGGAcgcg                                                                              >  2:256535/1‑143 (MQ=255)                                                                                   tCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGc                                                                        <  2:809707/148‑1 (MQ=255)                                                                                    cACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAAc                                                                    <  2:5809/151‑1 (MQ=255)                                                                                     aCCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTCACACCCCGGTGATTTTTCAGTTGTTTGATTTGAAAACTCACCGCTGACTGGGGCAGATAGAGGGATTTAGCCGCTCGACCAAAAAGAAGCGGTCGGCTaa                                                                     >  2:750779/1‑149 (MQ=255)                                                                                      ccGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑G‑AAAAAGGTCACACCCAGATGGTTTTTCTGTTGTTTGATTTGAAAACTCACTTCTTGCTGGGGCAGATAGAAGGAATCAACCGCCCCACCAAAATGACGCGGTCGGCTTACtt                                                                  >  1:285860/1‑151 (MQ=255)                                                                                       cGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTc                                                                 <  2:464024/151‑1 (MQ=255)                                                                                        ggCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑G‑AAAAGGTT‑ACACCCCGTTGATTTTTCAGTTGTTTGATTTGAAAGATCACTGCTGACTGGGGCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGGTAAATTcc                                                                >  2:663071/1‑151 (MQ=255)                                                                                           aGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGTGG‑AAAAGGTT‑ACACCCCGGTGGTTTTTCAGTTTTTTGGTTTGGAAGATCCCTGCTGAATGGGGCAGATAAAGGGATTTAGCCGCTCGACCCAAAGGAAGCGGTCGGCCAACtttccgg                                                             >  1:581318/1‑146 (MQ=255)                                                                                           aGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGC‑G‑AAAAAGGTTACCACCAGCTGATTTTTCAGTTGTCTGATTTGAAAAATCCCTGCTGACTGGGGCAGATAAAACGATTCAGCCCCTCGCCCCAAGAGATGCGGTTGGGTCACTTCCAgg                                                             >  1:284275/1‑151 (MQ=255)                                                                                               gTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAg                                                         <  2:1205567/151‑1 (MQ=255)                                                                                                    aCGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGttttt                                                    <  2:418493/151‑1 (MQ=255)                                                                                                    aCGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGttttt                                                    <  2:6958/151‑1 (MQ=255)                                                                                                          atTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACTCCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAAt                                              <  2:732315/151‑1 (MQ=255)                                                                                                          atTGTTTCTGTGGC‑CGGT‑G‑AAAAGGTTCACACCCAGTTGATTTACCAGTGGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATCGAGCAATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACGTCCAGGAAAGTTTTTAACa                                                <  1:161568/149‑1 (MQ=255)                                                                                                                tCTGTGGC‑GGGG‑G‑AAAAAGTTTACACCCAGGTGATTTTTCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGGCAGAAAAAGCGATTCAGCCGCTCGACCAAAGAGACGCGTTCGGCTAACTTCCAGGAAAGATTTTTACAATTCCGt                                         >  1:313904/1‑150 (MQ=255)                                                                                                                     ggC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCca                                    <  2:1233980/150‑1 (MQ=255)                                                                                                                       c‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCcac                                   >  2:958520/1‑149 (MQ=255)                                                                                                                       c‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCGCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCca                                    <  2:158862/148‑1 (MQ=255)                                                                                                                              g‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGct                             <  2:987213/150‑1 (MQ=255)                                                                                                                                 aaaGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCtctc                          >  2:778006/1‑151 (MQ=255)                                                                                                                                          cacCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCaaaa                   <  1:100589/149‑1 (MQ=255)                                                                                                                                          cacCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAatt                 >  1:879947/1‑151 (MQ=255)                                                                                                                                           acCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCAGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAatt                 >  2:1138485/1‑150 (MQ=255)                                                                                                                                               aGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTg            <  2:255326/151‑1 (MQ=255)                                                                                                                                                 ttGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGt           >  1:36550/1‑150 (MQ=255)                                                                                                                                                        tCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATg    <  2:219508/150‑1 (MQ=255)                                                                                                                                                        tCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATg    <  2:10692/150‑1 (MQ=255)                                                                                                                                                         ccAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAt  <  2:1252425/151‑1 (MQ=255)                                                                                                                                                         ccAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAt  <  2:901355/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                       GTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑G‑AAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATGAT  >  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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑