Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,706,273 T→C S7P (TCC→CCC)  srmB → ATP‑dependent RNA helicase SrmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,706,2730TC95.0% 56.4 / ‑6.2 20S7P (TCC→CCC) srmBATP‑dependent RNA helicase SrmB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (13/6);  minor base A (0/1);  total (13/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 = 1.00e+00

TACAAACAGATGAAGATTGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTTCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCAGGCTGC  >  NZ_CP009273/2706193‑2706352
                                                                                |                                                                               
tACAAACAGATGAAGATTGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTg                                                                        >  1:276484/1‑90 (MQ=255)
    aaCAGATGAAGATTGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACt                                                                    >  2:64077/1‑90 (MQ=255)
       aGATGAAGATTGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTcga                                                                 >  2:229375/1‑90 (MQ=255)
               aTTGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGcc                                                         <  2:276484/90‑1 (MQ=255)
                ttGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCct                                                        >  2:121575/1‑90 (MQ=255)
                         tAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGc                                               <  2:431960/90‑1 (MQ=255)
                         tAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGc                                               <  2:368525/90‑1 (MQ=255)
                         tAACGTCTATAATCAGCGCCCCAAACAGAGGTAGAACATGACTGTACCGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGc                                               <  2:341910/90‑1 (MQ=255)
                           aCGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGccc                                             <  1:266401/90‑1 (MQ=255)
                              tCTATAATCAGCGCCCCACACAAAGGTACACCATGCCGGCACCGGCTTTTACCGAAATTGAACCCGACGAAAGCCTGCTGGACGCCCTcc                                          <  2:218495/90‑1 (MQ=255)
                                      cAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTTGACGAAAGCCTGCTGGAAGCCCTCCAGGATaaa                                  >  2:113076/1‑90 (MQ=255)
                                      cAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATaaa                                  >  1:176218/1‑90 (MQ=255)
                                       aGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAg                                 >  2:145289/1‑90 (MQ=255)
                                                        tAGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGAcc                <  1:139359/90‑1 (MQ=255)
                                                         aGAACATGACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCg               >  2:53309/1‑90 (MQ=255)
                                                                gACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCa        >  2:324001/1‑90 (MQ=255)
                                                                gACTGTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCAGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCa        >  2:367002/1‑90 (MQ=255)
                                                                    gTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCAGGCt    >  1:11842/1‑90 (MQ=255)
                                                                    gTAACGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCAGGCt    >  2:48044/1‑90 (MQ=255)
                                                                      aaCGACTTTTCCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCAGGCTgc  >  1:391895/1‑90 (MQ=255)
                                                                                |                                                                               
TACAAACAGATGAAGATTGCAGCCGTAACGTCTATAATCAGCGCCCCACACAGAGGTAGAACATGACTGTAACGACTTTTTCCGAACTTGAACTCGACGAAAGCCTGCTGGAAGCCCTCCAGGATAAAGGTTTCACTCGCCCGACCGCCATTCAGGCTGC  >  NZ_CP009273/2706193‑2706352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: