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% 38.2 / ‑1.2 19F11F (TTC→TTTyjcFpentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (9/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

ATTTATGTATTTTCCTGTCCGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACG  >  NC_000913/4282805‑4283277
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
atTTATGTATTTTCCTGTCCGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACAtatt                                                                                                                                                                                                                                 >  1:172072/1‑250 (MQ=255)
                  ccGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAgcgc                                                                                                                                                                                                              <  1:321444/251‑1 (MQ=255)
                                                                                     tctCACTTTTTCATGATCATCCATTTTGAAAAGAGAGCGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAACACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAATAAAGAAAAGACATATTCTCTGCATTATAATTTTCTACGACTg                                                                                                                                           >  2:148484/3‑251 (MQ=255)
                                                                                               tCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATta                                                                                                                                                                                                                    >  1:219010/1‑168 (MQ=255)
                                                                                               tCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATta                                                                                                                                                                                                                    <  2:219010/168‑1 (MQ=255)
                                                                                                   gATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAACGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCgatc                                                                                                                              >  2:138276/1‑247 (MQ=255)
                                                                                                                              cTTTATTGTCAATTAACGTGCTTTTATGATGTTTGCTGTACTTATCTGATTTTGTTTTCTTTGTATGTTATGGACTTGTGACAGAGTGGTTTTAGTTAATCAGGCAAAGAGGTAAAGACATATTATTTAAACAATTATAGCGCATAAAATATCCTTATAATTAACACGAAAGAAAAGGCATATTCTCAGCATTATCATTTTATACGACTGACAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAat                                                                                                   >  2:297703/1‑250 (MQ=255)
                                                                                                                                     gTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATggt                                                                                           >  1:91023/1‑250 (MQ=255)
                                                                                                                                     gTCAATTAACGTGTTTTGATGATGGTTTCTGTTATTATCTGTGTTTGTTTTCTTGGTATGGGATGGACTTGTGACATAGTGGTTATCGTTAATCAGGCAAAGAGTAAATAACATATTATTTAAACCATTATAGCTCATAAAATATCCTTATAATTAACAATAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGATTGTCAAAAATCGCTCATTTTTTAATTAGGTTATTTGTTTAATATTATgg                                                                                            <  2:285203/250‑1 (MQ=255)
                                                                                                                                          ttAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGaaa                                                                                        >  2:178326/1‑249 (MQ=255)
                                                                                                                                            aaCGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGt                                                                                    <  1:12963/251‑1 (MQ=255)
                                                                                                                                               gtgtTTTTATGATGTTTGCTGTAATTATCTGCTTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGAGAAAAGGTGAt                                                                                 >  1:268141/1‑251 (MQ=255)
                                                                                                                                                 gtTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTTTTTAATATTATGGGAAAAGGTGATg                                                                                <  1:112524/250‑1 (MQ=255)
                                                                                                                                                                     aaTTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGAAAGCGTGGTTATCGTTATTCAGGCAAAGAGGAAAAAACATATTCTTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACCAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGAATGTCAAAAATCGCTCATTTTTTATTGAGTTTATTTGTTTCATATTATGGGAAAAGGTGATGCTTTTGGGAGAGGAAGAGttt                                                           >  2:121458/1‑249 (MQ=255)
                                                                                                                                                                                        cTTTGTCTGTGAAGGAGATGTGACATCGTGGTTAGCGTTAATCAGACAAAGCGGAAAAAAAATAAAATTTAAACCATTATAGCGCATAGAATATACTTAGAATAAACAAGTAAGAAACGGCATATTCTCAGCATTATCATTTTCAACGACTGTCAAACATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGAGATGCATTTGGTAGATGTAGAGTATTCCCCGGCCAgactaaagg                                        >  2:480243/1‑245 (MQ=255)
                                                                                                                                                                                           tGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCgggg                                      <  2:243062/250‑1 (MQ=255)
                                                                                                                                                                                                    tGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTa                                                                                                             <  2:105374/170‑1 (MQ=255)
                                                                                                                                                                                                    tGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGAAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTa                                                                                                             >  1:105374/1‑170 (MQ=255)
                                                                                                                                                                                                                               aaTCAGGCAAAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACg  <  1:178326/250‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ATTTATGTATTTTCCTGTCCGTTTCTTTCGTGATTTTTCCTGATGGTAAACTATCATAAATTCTTTTCATTTTCTGAATCAAGCCTATCACTTTTTCATGATCATCCATTTTGAAAAGAGAGAGGGCTTTATTGTCAATTAACGTGTTTTTATGATGTTTGCTGTAATTATCTGATTTTGTTTTCTTTGTATGTGATGGACTTGTGACAGAGTGGTTATCGTTAATCAGGCAAAGAGGGAAAAACATATTATTTAAACCATTATAGCGCATAAAATATCCTTATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACG  >  NC_000913/4282805‑4283277

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