Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,778 T→C S3S (AGT→AGC mprA → DNA‑binding transcriptional repressor MprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,7780TC76.2% 35.0 / 10.1 21S3S (AGT→AGCmprADNA‑binding transcriptional repressor MprA
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base C (4/12);  total (5/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAG  >  NC_000913/2810701‑2810844
                                                                             |                                                                  
aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTc                                                                  <  2:469308/80‑1 (MQ=255)
   cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGttt                                                              <  1:469309/81‑1 (MQ=255)
   cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGttt                                                              <  2:469310/81‑1 (MQ=255)
         aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGccc                                                        <  2:469311/81‑1 (MQ=255)
            tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCAt                                                      <  1:469312/80‑1 (MQ=255)
                 tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAAc                                                 >  2:469344/1‑80 (MQ=255)
                         aCTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTa                                         <  2:469313/80‑1 (MQ=255)
                             tCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTaaaa                                      <  2:469314/79‑1 (MQ=255)
                               gTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                  <  2:469315/81‑1 (MQ=255)
                                 tACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                >  2:469347/1‑81 (MQ=255)
                                           gCTGAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                      <  1:469316/81‑1 (MQ=255)
                                              gAAATTAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                    >  2:469343/1‑80 (MQ=255)
                                               aaaTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                    <  2:469317/79‑1 (MQ=255)
                                                  ttAATGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa               <  1:469318/81‑1 (MQ=255)
                                                      tGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt           <  1:469319/81‑1 (MQ=255)
                                                      tGAGGTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt           >  2:469345/1‑81 (MQ=255)
                                                         ggTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt        >  2:469346/1‑81 (MQ=255)
                                                          gTCATACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt       <  1:469320/81‑1 (MQ=255)
                                                             aTACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt     <  1:469321/80‑1 (MQ=255)
                                                              tACCCAAATGGATAGCTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa   <  1:469322/81‑1 (MQ=255)
                                                                cccAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg  <  2:469323/80‑1 (MQ=255)
                                                                             |                                                                  
AAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAG  >  NC_000913/2810701‑2810844

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