Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,415,8060GT44.7% 59.1 / 16.6 38N527K (AAC→AAA) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base G (20/1);  new base T (0/17);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 6.25e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATCGCGAGTCTGGATCCCCTT  >  NC_000913/1415656‑1415949
                                                                                                                                                       |                                                                                                                                               
gcgATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAg                                                                                                                                                 <  2:125826/151‑1 (MQ=255)
     gCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATc                                                                                                                                            >  1:691031/1‑151 (MQ=255)
                  aGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTgg                                                                                                                               >  1:718050/1‑151 (MQ=255)
                     aTTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCg                                                                                                                            >  1:907534/1‑151 (MQ=255)
                          ctttGTGCGGTTGTTCAGCCCCTTCCTTTCCCCCCGGTTCTTTTCCCTTACCTCCACCCTGCCTCTTTTCATCGATTTTTCCCGGTTGGGTG‑GCTCCCCCCACCCCCCCCCCACAGTCAGGTCATTTACCCCCGCCAGTCTGGCCGCAg                                                                                                                         <  2:826970/147‑1 (MQ=255)
                                 cGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGc                                                                                                                  >  2:713331/1‑149 (MQ=255)
                                            ccccTCCCTTTCCCCCCGTATCTTTTTCCTTATCCCCACCCTCCCTCTCTTCACCAAATTTTCCCTGGTTTTTG‑GCCCCCCCCCCCCCCCCCCCACAGTCAGGGAATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTcc                                                                                                       <  1:320206/147‑1 (MQ=255)
                                             gcCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGt                                                                                                    >  1:876435/1‑151 (MQ=255)
                                                ttCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTcc                                                                                                 >  1:281859/1‑151 (MQ=255)
                                                        cccccGGATTTTTTTCTTTACCTCCACCCTGCCTCTTTCCATCGATTTTTCCCTGGTTTGTT‑GCGCCCCCCCCCCCCCCCCCACAGCCGGGCCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCac                                                                                         <  2:384228/148‑1 (MQ=255)
                                                           ccGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCactac                                                                                      >  1:500300/1‑151 (MQ=255)
                                                                 cATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTc                                                                                >  1:685016/1‑151 (MQ=255)
                                                                  aTTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTc                                                                                >  2:501227/1‑150 (MQ=255)
                                                                         ttAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGAt                                                                        >  2:746816/1‑151 (MQ=255)
                                                                         ttACCTCCACCCTGCCTTTTTTCATCG‑ATTTTTCCTTGTTTTTTGGCGCCCCCCTCCCCCCCCCCACGGGCAGGGAATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTg                                                                          <  1:499078/149‑1 (MQ=255)
                                                                                  ccctccctctctTCATCGATTTTTCCCTGTTTTGTG‑GCGCCCCCCACCCCCCCACCAGAGTCGGGGCTTTTCTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTg                                                               <  2:112107/145‑1 (MQ=255)
                                                                                          ctctTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGtt                                                       >  1:100833/1‑151 (MQ=255)
                                                                                             ttCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGttttt                                                    >  2:647132/1‑151 (MQ=255)
                                                                                               cagcgattttttCCTGGTATTTG‑GCGCCGCCCACCCCCGCCCCACGGTCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGttttt                                                    <  2:114546/140‑1 (MQ=255)
                                                                                               cATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTc                                                  >  1:147012/1‑151 (MQ=255)
                                                                                                   gAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGtt                                              >  1:769725/1‑151 (MQ=255)
                                                                                                   gAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGtt                                              >  2:118094/1‑151 (MQ=255)
                                                                                                              gttttgtg‑GCGCCCCCCACCCCCCCCCCACGGTCAGGGCTTTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGt                                    <  2:910090/146‑1 (MQ=255)
                                                                                                              ggTTTGTT‑GCGTCCCCCATCCCCGCACCACAGTCGGGCCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctggctg                                      <  1:518234/148‑1 (MQ=255)
                                                                                                               ggtttgtgccgccccccACCCCCGCCCCACGGTCGGGGCAGTTCTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGtt                                   <  2:519674/137‑1 (MQ=255)
                                                                                                                tATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCt                                 >  1:626321/1‑151 (MQ=255)
                                                                                                                 aTGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTg                                >  1:321127/1‑151 (MQ=255)
                                                                                                                     t‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTac                            >  1:522267/1‑151 (MQ=255)
                                                                                                                        cgccGCCCACCCCCCCACCACAGTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCt                                 <  2:466888/141‑1 (MQ=255)
                                                                                                                           cccccACCCCCGCACCCCAGTGGGGGCTTTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTac                            <  1:430297/144‑1 (MQ=255)
                                                                                                                           cGCCCACCCCCGCACCACAGTAAGGGCGTTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTacaca                         <  2:885417/149‑1 (MQ=255)
                                                                                                                            ccccACCCCCGCACCACAGTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACaa                      <  1:799449/150‑1 (MQ=255)
                                                                                                                             cccccccccccaccacAGTCAGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACaa                      <  2:767344/140‑1 (MQ=255)
                                                                                                                              ccaccccccccccacATTCGGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATc                    <  2:339768/140‑1 (MQ=255)
                                                                                                                                 cccccGCACCACAGCCGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATcgcg                 <  2:387798/148‑1 (MQ=255)
                                                                                                                                          cacaGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATCGCGAGTCTGGAt        >  2:465134/1‑151 (MQ=255)
                                                                                                                                             aGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATCGCGAGTCTGGATccc     >  1:915490/1‑151 (MQ=255)
                                                                                                                                               tcggggCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATCGCGAGTCTGGAt        <  1:550140/143‑1 (MQ=255)
                                                                                                                                                  gggCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATCGCGAGTCTGGATCCCCtt  <  1:177566/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
GCGATGCGGATCGCTGCCAGCATTCTCATTGTGCGGATGTTCAGCGCCTTCCATTTCCTCCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTT‑GCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTGGCTGGTTCTGGTACACACAATCGCGAGTCTGGATCCCCTT  >  NC_000913/1415656‑1415949

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