Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 856,402 A→C E147A (GAA→GCA)  ybiT → ABC‑F family putative regulatory ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913856,4020AC93.8% 51.0 / ‑7.0 16E147A (GAA→GCA) ybiTABC‑F family putative regulatory ATPase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (4/11);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCGCATCTATGCTTTGCCGGAAATGAGTGAAGAAGACGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGAACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGATATTCTCCTGCTCGACGAACCGACCAACAACCTCGACATCGACACCATTCGCTGGCTGGAACAGG  >  NC_000913/856259‑856551
                                                                                                                                               |                                                                                                                                                     
cgACCGCATCTATGCTTTGCCGGAAATGAGTGAAGAAGACGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCa                                                                                                                                                 <  1:130848/150‑1 (MQ=255)
   ccGCATCTATGCTTTGCCGGAAATGAGTGAAGAAGACGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTAc                                                                                                                                             <  1:127644/151‑1 (MQ=255)
                    cGGAAATGAGTGAAGAAGACGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGt                                                                                                                            <  1:198552/151‑1 (MQ=255)
                                   aagaCGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGaa                                                                                                             <  2:308663/151‑1 (MQ=255)
                                      aCGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGc                                                                                                           <  1:76804/150‑1 (MQ=255)
                                                           tGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGgc                                                                                     <  2:99249/151‑1 (MQ=255)
                                                                gTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAGTTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGgcg                                                                                    >  1:181701/1‑147 (MQ=255)
                                                                gTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGgcg                                                                                    >  1:183356/1‑147 (MQ=255)
                                                                     aTACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCt                                                                           <  1:260255/151‑1 (MQ=255)
                                                                            gAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCgg                                                                    >  1:61635/1‑151 (MQ=255)
                                                                              aaTGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGat                                                                  <  2:37311/151‑1 (MQ=255)
                                                                                    cGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGCACACCACTACGGCCCGCTGAGTGAAGTTGTTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGAtcttct                                                             >  2:266737/1‑150 (MQ=255)
                                                                                                            cGGTGAACTGTTGCTTGGCGTGGGAGTTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGATATTCTCCTGCTCGACGAACCGACCaacaac                                    <  2:181700/151‑1 (MQ=255)
                                                                                                                  aCTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGATATTCTCCTGCTCGACGAACCGACCAACAACCtcgac                              <  1:279094/151‑1 (MQ=255)
                                                                                                                   cTGTTGCTTGGCGTGGGAATTCCAGTGGCACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGATATTCTCCTGCTCGACGAACCGACCAACAACCtcgaca                             <  1:229262/151‑1 (MQ=255)
                                                                                                                                              gAACAGCACTACGGCCCAATGAGTGAAGATGCTCCTGGATGGAACCTGCGGGTGCTTCTGGCGAAGGCCCTGTTTGCTGATCCGGATATTCTCCTGCACGACGAACCGACCACCCACCTCGACATGGACACCAGTCGCTGGCTGCAGCAgg  >  2:228422/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CGACCGCATCTATGCTTTGCCGGAAATGAGTGAAGAAGACGGCTATAAAGTGGCCGATCTGGAAGTTAAATACGGCGAAATGGACGGTTACTCTGCGGAAGCTCGCGCCGGTGAACTGTTGCTTGGCGTGGGAATTCCAGTGGAACAGCACTACGGCCCGATGAGTGAAGTTGCTCCTGGCTGGAAGCTGCGTGTGCTTCTGGCGCAGGCGCTGTTTGCTGATCCGGATATTCTCCTGCTCGACGAACCGACCAACAACCTCGACATCGACACCATTCGCTGGCTGGAACAGG  >  NC_000913/856259‑856551

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