Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 883,127 C→G G83G (GGC→GGG rcdA → DNA‑binding transcriptional regulator RcdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273883,1270CG90.0% 56.8 / 0.6 20G83G (GGC→GGGrcdADNA‑binding transcriptional regulator RcdA
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base G (7/11);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGCGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGCTGATGTACCAGCTCTACGCGC  >  NZ_CP009273/883042‑883212
                                                                                     |                                                                                     
tgttACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCAt                                                                                   >  1:210909/1‑90 (MQ=255)
         aGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTa                                                                          <  2:51422/90‑1 (MQ=255)
                   cagTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGatat                                                                <  2:293703/90‑1 (MQ=255)
                        ttACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGGGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCt                                                           >  1:702328/1‑90 (MQ=255)
                               gaTCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTc                                                    >  2:656334/1‑90 (MQ=255)
                               gaTCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTc                                                    >  1:494909/1‑90 (MQ=255)
                                         cGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCa                                          >  2:425517/1‑90 (MQ=255)
                                           gCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAAc                                        <  2:168717/90‑1 (MQ=255)
                                             aaTATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGCGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGc                                      <  2:182523/90‑1 (MQ=255)
                                               tatCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCg                                    <  2:210909/90‑1 (MQ=255)
                                                  cAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACCGGTTGCAACGCCGGAt                                 <  2:95942/90‑1 (MQ=255)
                                                       aTTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACAt                            <  2:39105/90‑1 (MQ=255)
                                                       aTTTTTTAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACAt                            <  2:432244/90‑1 (MQ=255)
                                                        ttttttAGCGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATg                           >  1:500050/1‑90 (MQ=255)
                                                           tttAGCGATGTTAGTGATGCTCCGGGCGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAg                        <  1:154516/90‑1 (MQ=255)
                                                            ttAGCGATTTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGc                       <  1:83422/90‑1 (MQ=255)
                                                                cGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGCTGAt                   <  1:221708/90‑1 (MQ=255)
                                                                cGATGTTAGTGATGCTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGCTGAt                   <  2:554508/90‑1 (MQ=255)
                                                                              cTCCGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGCTGATGTACCAGCTCTAcg     >  2:421708/1‑90 (MQ=255)
                                                                                 cGGGGGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGCTGATGTACCAGCTCTACgcgc  <  1:143901/90‑1 (MQ=255)
                                                                                     |                                                                                     
TGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTCCGGGCGCATGCCAGGCTATCACCGATATGATCTACAGCTCACAGGTTGCAACGCCGGATAACATGGAGCTGATGTACCAGCTCTACGCGC  >  NZ_CP009273/883042‑883212

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