Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,7200TC66.7% 24.1 / 19.0 24intergenic (+91/+63)yjcE/yjcFputative transporter YjcE/pentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base T (6/2);  new base C (8/8);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCA  >  NC_000913/4281643‑4281799
                                                                             |                                                                               
tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGActct                                                                              <  1:804374/81‑1 (MQ=255)
   ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTaa                                                                            >  1:804403/1‑80 (MQ=255)
     aGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAAc                                                                         <  2:804376/81‑1 (MQ=255)
     aGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAAc                                                                         <  2:804375/81‑1 (MQ=255)
      ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACa                                                                        <  2:804377/81‑1 (MQ=255)
         ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACaaat                                                                     <  1:804378/81‑1 (MQ=255)
                 aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTAt                                                              <  1:804379/80‑1 (MQ=255)
                 aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTAt                                                              <  2:804380/80‑1 (MQ=255)
                    tGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTAtctct                                                          >  2:804421/1‑81 (MQ=255)
                         aaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCttt                                                      >  2:804414/1‑80 (MQ=255)
                               cAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTc                                               <  2:804381/81‑1 (MQ=255)
                                aaTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCCCTCTTTTAAAGTc                                               >  1:804418/1‑80 (MQ=255)
                                 aTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCaa                                             >  1:804404/1‑81 (MQ=255)
                                       tGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAtt                                       >  1:804413/1‑81 (MQ=255)
                                       tGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAtt                                       >  1:804400/1‑81 (MQ=255)
                                        gCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTg                                      >  2:804411/1‑81 (MQ=255)
                                                 aCAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttctt                             >  1:804412/1‑81 (MQ=255)
                                                 aCAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttctt                             >  2:804415/1‑81 (MQ=255)
                                                 aCAGGTTAGCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttctt                             >  1:804410/1‑81 (MQ=255)
                                                         gCTCGGGCTAACATTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTaa                     <  2:804382/81‑1 (MQ=255)
                                                                     aTTCATGACTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa         >  1:804420/1‑81 (MQ=255)
                                                                      ttCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCAt        <  2:804383/81‑1 (MQ=255)
                                                                        cATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCt      >  2:804384/1‑81 (MQ=255)
                                                                           gATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTg    >  2:804424/1‑80 (MQ=255)
                                                                            aTTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCa  <  2:804385/81‑1 (MQ=255)
                                                                             |                                                                               
TAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCA  >  NC_000913/4281643‑4281799

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