Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,384 T→C G74G (GGT→GGC dnaK → chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,3840TC95.0% 47.1 / ‑3.9 21G74G (GGT→GGCdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (13/6);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

CCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAC  >  NC_000913/12138‑12603
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           
ccGAATAAATAGTGGAGACGTGTAGATGGGTAAAATAATTGGTATCGTCCTTGGTACTACCAACTCTTGTGGAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGGCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAACCGCCTGATTGgccgcc                                                                                                                                                                                                                         <  2:132076/251‑1 (MQ=255)
             ggAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACg                                                                                                                                                                                                             <  1:179098/250‑1 (MQ=255)
              gagaCGTTTAGATGGGCAAAATAAATGGTATCGCCCTGGGTACTACCACCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGATGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGACTATACCCAGGATGGTGAAACTCTAGTTGGTCAGACGGCTAAACGTCCGGCAGTGACGAACCAGAAAACAACCCTGTTTGCGATTCAACGCCTGATTGGACACCGCTTCCAGGACga                                                                                                                                                                                                            >  2:219083/1‑250 (MQ=255)
                           gggTAAAATAATTGGTATAGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                                                                                                                                                              >  1:67164/1‑251 (MQ=255)
                                      ttGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTATACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGCTTGGCCGCCGATTCCAGGACGGATAAGTActggt                                                                                                                                                                                                >  2:5340/1‑234 (MQ=255)
                                      ttGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCg                                                                                                                                                                                                <  1:5340/238‑1 (MQ=255)
                                                             caCTCGTGTGTAGCGATTATGGATTGCAACACTTCTCGTGTGCTGGAGAACGTCGTAGGCGATCGAATCTTGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCTGCTAAACGTCAGGCAGTGACGAACCCGCAAATCACCCTGTTTGCTATTAAATGCCTGATTGGCAGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCCTCATGCCGTTCAAAATTATtgctgc                                                                                                                                                            <  1:277433/250‑1 (MQ=255)
                                                                                 ggatggCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATgg                                                                                                                                         <  1:48756/250‑1 (MQ=255)
                                                                                                   tgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGAt                                                                                                                                                        >  1:381721/2‑217 (MQ=255)
                                                                                                   tgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGAt                                                                                                                                                        <  2:381721/216‑1 (MQ=255)
                                                                                                                 ccGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAACATTATTGCTGCTGATAACGGCGACGCATGGGCCGACGTTAAAGTCCAGAAAATGGCAccgccg                                                                                                        >  2:134485/1‑251 (MQ=255)
                                                                                                                             gCACCACGCCTTCTATAATTGCCTATACCCAGGATGGTGAAACTCTGATTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGAAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTTCAGCGTGATGTTTCCATCATGCCGTTCAAACTTATTGCTGCTGAGAACGGAGACGCATGGTTCGCACTTAAAGGCCAGACAATGGCACCGCGGCAAATCTctgct                                                                                            >  2:118430/1‑251 (MQ=255)
                                                                                                                                   cGCCTTCTATCATTGCATATACCCAGGATGGTGTAACTCTGGTTGGTCAGCAGGCTAAACGTCAGGCAGTGACGAACCCGAAAAACACCATGTTTGCGATTCACCGCATGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATAATGCCGTTCAAAATTATTGCTGCTGCTAACGGCGACGCATGGGTCGACGTTAACGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAg                                                                                        >  1:375766/1‑249 (MQ=255)
                                                                                                                                                           aGGATGGTGAAACCTGGGTTGGTCAGCCGGCAAAACGTCAGGCATGGACGACCCCGCAAAACCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCACCTTGCCGTACAAAATTATGGCTGCTGATAACGGCGACGCAGGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCt                                                              <  1:168419/251‑1 (MQ=255)
                                                                                                                                                              aTGGAGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATATCTGCTGAAGAGCTGAAAAAAATGAAGAAAACCGCTGa                                                            >  1:151618/1‑250 (MQ=255)
                                                                                                                                                                 gTGAAACTCTGGTTGGTCAGCCGGCTAAACGTAAGGAAGTGACGAACCCGCAACACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAAAAGAAAAGCCTGACGTTTCCATCATGCAGTTAACAATTATTGCTGCTGATAACGGAGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCTGCAGATCTCTGATGAAGTGCGGAAAAAAATGAAGAAAACCGCTGa                                                            >  1:368451/1‑247 (MQ=255)
                                                                                                                                                                         ctggTTGGTCAGCCGGCTAAACGACAGGCAGTGACGAACCCGCCAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGACAATGGCACCGCCGCAGATCTCTTCTGAAGTGCTGAAAAAAATGAATAAACCCGCTGAAGATTGCCTggg                                                >  2:206106/4‑251 (MQ=255)
                                                                                                                                                                                                       aGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATATCTGCTGAAGTGCTGaaaaaac                                                                            >  2:292988/1‑193 (MQ=255)
                                                                                                                                                                                                       aGTGACGAACCCGCAAAACACCCTGTTTGCGATAAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGaaaaaaa                                                                            <  1:292988/194‑1 (MQ=255)
                                                                                                                                                                                                            cGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCg              >  1:294849/1‑251 (MQ=255)
                                                                                                                                                                                                                        acacCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAc  >  1:146665/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ttAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGaa                                                                                         >  1:167510/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                           
CCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAC  >  NC_000913/12138‑12603

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