Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0400CT35.1% 30.1 / 35.5 37F273F (TTC→TTT) ‡ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base C (5/19);  new base T (0/13);  total (5/32)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATA  >  NC_000913/1373976‑1374102
                                                                |                                                              
tctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGGTCTtc                                                                >  1:1791214/1‑65 (MQ=255)
tctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGCTCTtc                                                                >  1:813493/1‑65 (MQ=255)
gctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtc                                                                <  1:2140293/64‑1 (MQ=255)
   aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtcatc                                                             >  1:3504823/1‑65 (MQ=255)
   aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGCTCTtcatc                                                             >  1:1977462/1‑65 (MQ=255)
               ttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCt                                                 <  1:840686/65‑1 (MQ=255)
                 gcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGt                                               <  1:4350017/65‑1 (MQ=255)
                 gcctttAATGATTGAGGGGGAGGCGTGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTc                                              <  1:2522784/65‑1 (MQ=255)
                  gctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                              <  1:1522866/64‑1 (MQ=255)
                  cctttAATGCTTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                              <  1:1112886/65‑1 (MQ=255)
                  cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                              <  1:186830/65‑1 (MQ=255)
                  cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                              <  1:4062751/65‑1 (MQ=255)
                  cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                              <  1:3178950/65‑1 (MQ=255)
                    tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCAt                                            <  1:2445134/65‑1 (MQ=255)
                      tAATGATTAGGGAAAAGGGGCGGCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTct                                         <  1:3600184/65‑1 (MQ=255)
                       aaTGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct                                         <  1:1276180/65‑1 (MQ=255)
                        aTGATTTAGGAAGAGGGGCGGCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTctgc                                       <  1:2506311/65‑1 (MQ=255)
                           aTTTAGAAAAAGCGGCGGCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGa                                    <  1:4260352/65‑1 (MQ=255)
                             ttAGGAGAAGGCGCGTCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGAcc                                  <  1:943094/65‑1 (MQ=255)
                                 gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                               <  1:907045/65‑1 (MQ=255)
                                 gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                               <  1:2521374/65‑1 (MQ=255)
                                 gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                               <  1:3212496/65‑1 (MQ=255)
                                      gcggcggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACg                          <  1:2581001/65‑1 (MQ=255)
                                            tagATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAc                   <  1:1373450/63‑1 (MQ=255)
                                            gcgATCTCGGTCGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAGCGGGTGaa                    <  1:3741038/65‑1 (MQ=255)
                                             cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:1258388/65‑1 (MQ=255)
                                             agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:1995763/64‑1 (MQ=255)
                                             agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:3809222/64‑1 (MQ=255)
                                             agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:2948070/64‑1 (MQ=255)
                                             agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:528010/64‑1 (MQ=255)
                                             agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAAAGGGTGaaaa                  <  1:3062859/64‑3 (MQ=255)
                                             agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGAACAGGAAAAGGGTGaaaa                  <  1:3012492/64‑3 (MQ=255)
                                              gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:2982532/65‑1 (MQ=255)
                                              gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:2418223/65‑1 (MQ=255)
                                              gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                  <  1:2681582/65‑1 (MQ=255)
                                                   cGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAAt             >  1:1968269/1‑65 (MQ=255)
                                                       ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca         >  1:4300517/1‑65 (MQ=255)
                                                         tgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaac       <  1:51684/65‑1 (MQ=255)
                                                         tgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaac       >  1:4185407/1‑65 (MQ=255)
                                                          gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaaca      >  1:2508423/1‑65 (MQ=255)
                                                           ctctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaacaa     <  1:457000/65‑1 (MQ=255)
                                                            tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa    >  1:3618812/1‑65 (MQ=255)
                                                              ttCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATa  <  1:3307927/65‑1 (MQ=255)
                                                                |                                                              
TCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATA  >  NC_000913/1373976‑1374102

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