Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,613,378 T→C M1A (GTG→GCG) † marR → multiple antibiotic resistance transcriptional regulator MarR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,613,3780TC97.3% 135.3 / ‑3.5 37M1A (GTG→GCG) †marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

TCAATTCATTCATTTGACTTATACTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGTGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAGAAGAAAGATCGCC  >  NZ_CP009273/1613301‑1613458
                                                                             |                                                                                
tCAATTCATTCATTTGACTTATACTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAg                                                                      <  1:326710/90‑1 (MQ=255)
tCAATTCATTCATTTGACTTATACTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAg                                                                      <  1:456015/90‑1 (MQ=255)
             ttGCCTTATACTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAAt                                                         <  2:526849/90‑1 (MQ=255)
             ttGACTTATACTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAAt                                                         >  2:199826/1‑90 (MQ=255)
                     taCTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGTGAAAAGTACCAGCGATCTGTTCAATGAAattat                                                 <  1:231618/90‑1 (MQ=255)
                     taCTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAattat                                                 <  2:428447/90‑1 (MQ=255)
                     taCTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAattat                                                 <  2:1054869/90‑1 (MQ=255)
                         tGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCa                                             >  2:672968/1‑90 (MQ=255)
                         tGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGAACCAGCGATCTGTTCAATGAAATTATTCCa                                             >  2:1112327/1‑90 (MQ=255)
                              gggCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTggg                                        <  2:94398/90‑1 (MQ=255)
                                    tattatCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCtt                                  <  1:199826/90‑1 (MQ=255)
                                      ttatCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTaa                                <  1:1187289/90‑1 (MQ=255)
                                        atCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATc                              >  1:1109736/1‑90 (MQ=255)
                                        atCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATc                              >  2:522469/1‑90 (MQ=255)
                                        atCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATc                              >  2:235498/1‑90 (MQ=255)
                                           cccTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCat                           <  1:382188/90‑1 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGaaa                                                                              >  2:782348/1‑37 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGaaa                                                                              <  1:782348/37‑1 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCatat                         >  1:955263/1‑90 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCatat                         >  2:609080/1‑90 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCat                           >  1:349672/1‑88 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCat                           >  1:385886/1‑88 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCat                           <  2:349672/88‑1 (MQ=255)
                                             cTGCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCat                           <  2:385886/88‑1 (MQ=255)
                                               gCAACTAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATgg                       >  1:3548/1‑90 (MQ=255)
                                                    tAATTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAAt                  >  1:179968/1‑90 (MQ=255)
                                                     aaTTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATc                 >  1:992209/1‑90 (MQ=255)
                                                     aaTTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATc                 <  2:695470/90‑1 (MQ=255)
                                                      aTTACTTGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCa                >  2:185164/1‑90 (MQ=255)
                                                            tGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTaa                                <  2:1086390/68‑1 (MQ=255)
                                                            tGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTaa                                >  1:1086390/1‑68 (MQ=255)
                                                            tGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAgaagaa          >  2:234189/1‑90 (MQ=255)
                                                            tGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAgaagaa          >  1:649555/1‑90 (MQ=255)
                                                            tGCCAGGGCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAgaagaa          >  2:118577/1‑90 (MQ=255)
                                                                 gggCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAGAAGAAAGATc     >  1:928014/1‑90 (MQ=255)
                                                                   gCAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAGAAGAAAGATCGc   <  2:1197149/90‑1 (MQ=255)
                                                                    cAACTAATGCGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAGAAGAAAGATCGcc  <  2:928014/90‑1 (MQ=255)
                                                                             |                                                                                
TCAATTCATTCATTTGACTTATACTTGCCTGGGCAATATTATCCCCTGCAACTAATTACTTGCCAGGGCAACTAATGTGAAAAGTACCAGCGATCTGTTCAATGAAATTATTCCATTGGGTCGCTTAATCCATATGGTTAATCAGAAGAAAGATCGCC  >  NZ_CP009273/1613301‑1613458

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