Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,283,043 G→A F11F (TTC→TTT)  yjcF ← pentapeptide repeat‑containing protein YjcF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,283,0430GA88.9% 33.6 / ‑3.6 18F11F (TTC→TTT) yjcFpentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (8/8);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AAAATATTTATGTATTTTCCTGTCCGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCG  >  NC_000913/4282800‑4283286
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
aaaaTATTGATGTATTTTCCTTTTCGTTTCTTTCGTGAATTTTCCTGATTGTAAACTATCATAAATTCTTTTCAGTTTCTGAATCAATCCTATCTCTTTTTCATTATCCTCCATTTTTAAAAGATATAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTTATTTTGTTTTCTTTGTATGTGATGGACTTGTTACAGAGTGGTTATCGGGAATCATGCAAATAGTAAAAAACa                                                                                                                                                                                                                                              <  2:331091/251‑1 (MQ=255)
         aTGTATTTTCCTGTCCGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGTTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTGTGTATGTGATGGACTTGTGACAGTGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTAttt                                                                                                                                                                                                                                      <  2:38618/250‑1 (MQ=255)
                  ccTGTCCGTTTCTTGCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGCTTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTATACCATtat                                                                                                                                                                                                                            >  1:284742/1‑251 (MQ=255)
                       ccGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCTAAGAGGAAAAAACATATTGTTTAAACCATTATAgcgc                                                                                                                                                                                                                       >  2:168487/1‑251 (MQ=255)
                                    gATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCtt                                                                                                                                                                                                          <  1:211595/251‑1 (MQ=255)
                                          tCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAAt                                                                                                                                                                                                     <  1:168487/250‑1 (MQ=255)
                                                 ggTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATAAATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACa                                                                                                                                                                                                >  1:113598/1‑248 (MQ=255)
                                                    caacgaTCATAAATTCTTTTCAGTTTCTGAATCAAGCCTATCTCTTTTTCATTATCATCCATTTTGAAAAGATAGAGGGCTTTATTGTTAATTAACGTGTTTTTATGATGTTTGCTTTAATTTTCTTATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCATGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACaagaaa                                                                                                                                                                                           <  2:321502/245‑1 (MQ=255)
                                                            aTAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGc                                                                                                                                                                                   >  1:263039/1‑250 (MQ=255)
                                                                                                    tCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTGTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTAGTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATc                                                                                                                                           <  1:132857/250‑1 (MQ=255)
                                                                                                                   ttgcaaaGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATg                                                                                                                            >  1:330989/5‑250 (MQ=255)
                                                                                                                                                 aaCGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGAATTGTGACAGAGTGGTTATCGTTAATCAGGCACAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGt                                                                                             >  1:39323/1‑251 (MQ=255)
                                                                                                                                                         tttATGATGTTTGATATAATTATCTGATTTTGTTTTCTTTGTCTGTGATGGAGTTGTGAAAGAGAGTTTATCGTTAATCAGGCAAAGAGGGAAAAACATATTATTTAAACAATTAAAGCGCATAAAATATCCTTATTATTAATAAGATAGAAAAGGCATATTCTATGCCTTATCATTTTATACAACTGTAAAAAATCGCTAATTTTTTAATGAGTTTATTTGTTTAATATTATTGTAAAAGGTAGTGCAtt                                                                                     >  2:276768/1‑251 (MQ=255)
                                                                                                                                                           tATGATGTTTGCTGTCATTATCTGATTTTGTATTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTg                                                                                   >  1:299150/1‑251 (MQ=255)
                                                                                                                                                                 gTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGCGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTTTGGGAAAAGGTGTTGCATTTGGgag                                                                               >  1:235588/1‑249 (MQ=255)
                                                                                                                                                                                                                 tggcatagTGGTTATCGTTCATCAGACAACGAGGAAAAAACATATTATTTCAACCATTATAGCGCATAAAATATACTTATAATTAACAAGAAAGAAAAGGCATGTTCTCTGCATTATCATTTTCTACGACGGTAAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTATAATTATGGGAAAAGGAGATTCCTTTGTGAGAGGTAGAGTATTCCCCGGTCTGCCGACCGGGGAAGGGTAGAGATTAAttc                             >  2:75742/7‑249 (MQ=255)
                                                                                                                                                                                                                                           cAAAGAGTAAAAAACATATTATTTAAACCATTATATCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTTGCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGcc   <  2:330989/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCg  <  1:28275/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
AAAATATTTATGTATTTTCCTGTCCGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCG  >  NC_000913/4282800‑4283286

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