Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,315,520 G→T R7S (CGC→AGC)  phnO ← aminoalkylphosphonate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,5200GT94.4% 54.4 / ‑2.6 18R7S (CGC→AGC) phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (6/11);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGC  >  NC_000913/4315384‑4315648
                                                                                                                                        |                                                                                                                                
gcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAg                                                                                                                               <  1:186143/140‑1 (MQ=255)
    gATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTca                                                                                                                           <  1:453532/140‑1 (MQ=255)
                         cgcgCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTtctc                                                                                                      <  2:388949/140‑1 (MQ=255)
                           cgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                                                                                    >  1:25623/1‑140 (MQ=255)
                           cgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                                                                                    >  2:150053/1‑140 (MQ=255)
                                  tGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATg                                                                                             >  1:416189/1‑140 (MQ=255)
                                    gCGTTAAAACCCACGCGAAACGCGTGGCGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGa                                                                                           <  1:296480/140‑1 (MQ=255)
                                         aaaaCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGc                                                                                      >  2:154244/1‑140 (MQ=255)
                                                         gcgTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcga                                                                      <  2:262662/140‑1 (MQ=255)
                                                                  tCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGTCCGACTGgc                                                             <  2:308863/140‑1 (MQ=255)
                                                                         ccGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCt                                                      <  2:503763/140‑1 (MQ=255)
                                                                         ccGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCt                                                      <  2:416189/140‑1 (MQ=255)
                                                                                               cAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATgg                                <  2:139486/140‑1 (MQ=255)
                                                                                                  cgcgTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCaa                             <  2:160308/140‑1 (MQ=255)
                                                                                                  cgcgTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCaa                             <  1:150053/140‑1 (MQ=255)
                                                                                                       aaaCCGCGTCGGTGTCGTACTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGGTGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCtg                        >  1:182339/1‑140 (MQ=255)
                                                                                                                         aCTGCGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGggcggc      <  1:154244/140‑1 (MQ=255)
                                                                                                                             cGTGGCCGGGCTAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGgcgcgc  >  2:152273/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                
GCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTGGCGCAGGCTGCCGTCATTGTTGAGCGTATGGCAATCCTGTGGAGTATAGCGGGCGGCGCGC  >  NC_000913/4315384‑4315648

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