Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,237 C→T R25R (CGC→CGT)  dnaK → chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,2370CT96.0% 66.0 / ‑3.9 25R25R (CGC→CGT) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (11/13);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

CATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGG  >  NC_000913/12004‑12468
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
cATCTCCACCTTGATGACGGGGGTTACGACCCCATTTAGTAGTCAAACGCCGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGGTGACCGAATATATAGTGGAGCCGTTAAGATGGGAAAAATAATTGGAATCGCCCTGAATACTACCATCTCTTGTGTAGCGATTATGGATGGCACCCCTCCTCGTGTGCAGGAGAccgccg                                                                                                                                                                                                                          >  2:124460/1‑250 (MQ=255)
 aTCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGCGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGa                                                                                                                                                                                                                         >  1:234037/1‑250 (MQ=255)
           tgatgaCGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGc                                                                                                                                                                                                               <  1:138299/250‑1 (MQ=255)
            gctgaCGTGGGTTACGACCACATTTAGTAGACAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCACTATTACAGACTCACAACCACATGATGACCGAAAATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATAGATGGCACCCCTCCTCGTGTACGGGAGAACGCCAAAGGCGATCGca                                                                                                                                                                                                              >  1:59317/3‑250 (MQ=255)
                    gggttACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCtt                                                                                                                                                                                                     <  1:337658/248‑1 (MQ=255)
                                     agtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGCAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATAcc                                                                                                                                                                                    >  1:128448/1‑251 (MQ=255)
                                                       gtgagtCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAActctg                                                                                                                                                                  <  1:340687/251‑2 (MQ=255)
                                                                     aaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGc                                                                                                                                                    <  1:124460/251‑1 (MQ=255)
                                                                         gAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGAATCACAACCACATGATGAACGAATATAAAGTGGAGACGTTTAGATGGGTAAAATACTTGGTATCGACCTGGGTACTCCCAACTCTTGTGAAGAGATTAAGGGTGGAACCAATCCTCGTGTGCTGGAGCACGACGCAGGCGATCGCACCAAGCCTTCTATCATTGCATATACCAAAGATAGTGAAACTCTGGTTGGTCatctggctca                                                                                                                                                 >  1:357684/1‑241 (MQ=255)
                                                                              tGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                            <  1:179346/250‑1 (MQ=255)
                                                                              tGGGCAGTTGAAACCAGACGTTTCGCCCCTAGTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGATAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                            <  1:226531/250‑1 (MQ=255)
                                                                              tGGGCAGTTGAAACCAGACGTTTCACCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTc                                                                                                                                            <  1:317912/250‑1 (MQ=255)
                                                                                ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAgg                                                                                                                                         <  1:215022/251‑1 (MQ=255)
                                                                                            cAGACGTTTCGCCCCTATTACAGACTCACAACCACCTGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGAAGGGACGAAcc                                                                                                                             >  1:302641/1‑251 (MQ=255)
                                                                                            cAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGTCGTTTGGCTGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTTTGCTGGAGAACGCCGACGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAACCTCTGGTTGTTCAGCCGACTAGACGTCCGGCCGTGACGAAcc                                                                                                                             >  1:359941/1‑251 (MQ=255)
                                                                                                                                  atgaCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGc                                                                                                                                                                                                                     >  2:74557/1‑125 (MQ=255)
                                                                                                                                  atgaCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCAATCCTCGTGTGCTGGAGAACGCCGAAGGc                                                                                                                                                                                                                     <  1:74557/125‑1 (MQ=255)
                                                                                                                                      ccGAATAAATAGTGGAGACGTGTAGATGGGTAAAATAATTGGTATCGTCCTTGGTACTACCAACTCTTGTGGAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGGCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAACCGCCTGATTGgccgcc                                                                                  <  2:132076/251‑1 (MQ=255)
                                                                                                                                                   ggAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACg                                                                      <  1:179098/250‑1 (MQ=255)
                                                                                                                                                    gagaCGTTTAGATGGGCAAAATAAATGGTATCGCCCTGGGTACTACCACCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGATGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGACTATACCCAGGATGGTGAAACTCTAGTTGGTCAGACGGCTAAACGTCCGGCAGTGACGAACCAG‑AAAACAACCCTGTTTGCGATTCAACGCCTGATTGGACACCGCTTCCAGGACga                                                                     >  2:219083/1‑250 (MQ=255)
                                                                                                                                                                 gggTAAAATAATTGGTATAGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                       >  1:67164/1‑251 (MQ=255)
                                                                                                                                                                            ttGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTATACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGCTTGGCCGCCGATTCCAGGACGGATAAGTActggt                                                         >  2:5340/1‑234 (MQ=255)
                                                                                                                                                                            ttGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCg                                                         <  1:5340/238‑1 (MQ=255)
                                                                                                                                                                                                   caCTCGTGTGTAGCGATTATGGATTGCAACACTTCTCGTGTGCTGGAGAACGTCGTAGGCGATCGAATCTTGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCTGCTAAACGTCAGGCAGTGACGAACCCGCAAATCA‑CCCTGTTTGCTATTAAATGCCTGATTGGCAGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCCTCATGCCGTTCAAAATTATtgctgc                     <  1:277433/250‑1 (MQ=255)
                                                                                                                                                                                                                       ggatggCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATgg  <  1:48756/250‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        
CATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGG  >  NC_000913/12004‑12468

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