Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,306,692 A→T I115N (ATC→AAC)  phnM ← alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,306,6920AT88.9% 51.7 / 0.3 18I115N (ATC→AAC) phnMalpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (6/10);  total (7/11)
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.07e-01

GCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT  >  NZ_CP009273/4306618‑4306777
                                                                          |                                                                                     
gCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGAt                                                                        <  2:253053/90‑1 (MQ=255)
             cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAg                                                                          <  2:78067/75‑1 (MQ=255)
             cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAg                                                                          >  1:78067/1‑75 (MQ=255)
             cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATc                                                           >  1:329657/1‑90 (MQ=255)
                  gCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATcgccgc                                                      <  2:1297/90‑1 (MQ=255)
                                     gTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAAt                                   >  1:121496/1‑90 (MQ=255)
                                               cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa                         <  2:511199/90‑1 (MQ=255)
                                               cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa                         >  1:380492/1‑90 (MQ=255)
                                                gTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCAt                        <  2:426497/90‑1 (MQ=255)
                                                  ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                      <  1:539762/90‑1 (MQ=255)
                                                  ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                      <  1:416513/90‑1 (MQ=255)
                                                  ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc                      <  2:421329/90‑1 (MQ=255)
                                                       cGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTAc                 <  2:121496/90‑1 (MQ=255)
                                                         tctcTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACgg               >  2:401423/1‑90 (MQ=255)
                                                          ctctTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACggt              <  2:380492/90‑1 (MQ=255)
                                                                 aTGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg                           >  1:342507/1‑70 (MQ=255)
                                                                 aTGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg                           <  2:342507/70‑1 (MQ=255)
                                                                    gCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCg    >  1:146855/1‑90 (MQ=255)
                                                                      gttgttCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCt  <  2:33760/90‑1 (MQ=255)
                                                                      gttgttCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCt  >  1:338241/1‑90 (MQ=255)
                                                                      gttgttCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCt  >  2:7959/1‑90 (MQ=255)
                                                                          |                                                                                     
GCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT  >  NZ_CP009273/4306618‑4306777

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