Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,328,195 (T)5→4 intergenic (+112/+44) yhbY → / ← greA ribosome assembly factor YhbY/transcription elongation factor GreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,328,1910T.96.4% 105.3 / ‑2.4 28intergenic (+108/+48)yhbY/greAribosome assembly factor YhbY/transcription elongation factor GreA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (12/15);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTG  >  NC_000913/3328060‑3328320
                                                                                                                                   |                                                                                                                                 
cGTAAAATCTCGCTGCCACGCTAAGATTATCCTGAAGTTTCACACATTGGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCcttt                                                                                                                          <  2:349038/140‑1 (MQ=255)
   aaaaTCTCGCTGCCACGCTAAGGTTATCCTAAAGTTAAACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCGGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCctttctt                                                                                                                       <  1:325236/140‑1 (MQ=255)
              gCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATc                                                                                                            <  2:624357/140‑1 (MQ=255)
                      aaGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCtt                                                                                                    <  1:90668/140‑1 (MQ=255)
                         aTTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                 >  2:467999/1‑140 (MQ=255)
                           tATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                               <  2:664133/140‑1 (MQ=255)
                           tATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                               <  1:467999/140‑1 (MQ=255)
                             tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                             >  2:771147/1‑140 (MQ=255)
                                   aGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTaa                                                                                       >  1:613634/1‑140 (MQ=255)
                                          cacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTAc                                                                                <  1:771147/140‑1 (MQ=255)
                                                 gCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                         <  2:197164/140‑1 (MQ=255)
                                                     tgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTCCAATACATTAACATGTCGAGTATTGTGTAACTCTTACAGGTATTCCAc                                                                     >  2:361580/1‑140 (MQ=255)
                                                     tgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTcc                                                                       <  1:414385/138‑1 (MQ=255)
                                                     tgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTcc                                                                       >  2:414385/1‑138 (MQ=255)
                                                       tAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAt                                                                          <  2:598745/133‑1 (MQ=255)
                                                       tAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAt                                                                          >  1:598745/1‑133 (MQ=255)
                                                       tAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCt                                                                   >  1:507323/1‑140 (MQ=255)
                                                               ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                         >  1:51158/1‑126 (MQ=255)
                                                               ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                         <  2:51158/126‑1 (MQ=255)
                                                                             gcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTc                                             >  1:456336/1‑140 (MQ=255)
                                                                             gcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTc                                             >  2:539471/1‑140 (MQ=255)
                                                                                       aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGc                                   <  2:613634/140‑1 (MQ=255)
                                                                                        aaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCg                                  <  1:61041/140‑1 (MQ=255)
                                                                                         aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTGCAGGTATTGCACGTTAGTTACTTCATATTCTACTTCGCCGCCCGGCGt                                 >  2:222835/1‑140 (MQ=255)
                                                                                               aCGCTCTATTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                           <  2:456336/140‑1 (MQ=255)
                                                                                                  ctctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGAc                        <  1:539471/140‑1 (MQ=255)
                                                                                                         cGTTGATAAAAGGCCACATAGCCGCTTTTCCCTTTCTCTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCacaaca                 <  1:222835/140‑1 (MQ=255)
                                                                                                                        gCATAGCGGCCTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg  >  2:586190/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                                 
CGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTG  >  NC_000913/3328060‑3328320

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: