Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,713,741 A→G I40T (ATC→ACC)  anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,713,7410AG78.6% 50.7 / 11.4 28I40T (ATC→ACC) anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (10/12);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGCCCGATTTCATAGTTCATCCT  >  NZ_CP009273/1713616‑1713869
                                                                                                                             |                                                                                                                                
tAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCa                                                                                                                      >  2:698275/1‑138 (MQ=255)
        gCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCaaa                                                                                                                    <  1:94553/132‑1 (MQ=255)
        gCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCaaa                                                                                                                    >  2:94553/1‑132 (MQ=255)
            tGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCtg                                                                                                         >  2:693009/1‑139 (MQ=255)
              aCCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCtgtg                                                                                                       >  2:487030/1‑139 (MQ=255)
                                        ttGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGcc                                                                             <  1:774590/139‑1 (MQ=255)
                                             tATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacac                                                                        <  1:487030/139‑1 (MQ=255)
                                             tATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacac                                                                        <  1:693009/139‑1 (MQ=255)
                                             tATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacac                                                                        <  1:698275/139‑1 (MQ=255)
                                                  aGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGc                                                                                                               <  2:844004/95‑1 (MQ=255)
                                                  aGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGc                                                                                                               >  1:844004/1‑95 (MQ=255)
                                                       tCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCaa                                                                <  1:262729/137‑1 (MQ=255)
                                                       tCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAaca                                                              >  1:70488/1‑139 (MQ=255)
                                                       tCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAac                                                               <  1:98432/138‑1 (MQ=255)
                                                       tCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAac                                                               <  1:639540/138‑1 (MQ=255)
                                                                 gAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGc                                                    >  1:437892/1‑139 (MQ=255)
                                                                   aaGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTg                                                  <  2:119861/139‑1 (MQ=255)
                                                                        tAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGcc                                             <  2:856362/139‑1 (MQ=255)
                                                                                tgGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATaa                                     >  1:550423/1‑139 (MQ=255)
                                                                                       ggCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGGTCGGCCAACTCAAACTTGCCAACTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCaa                               >  2:765222/1‑138 (MQ=255)
                                                                                                  aGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCa                                                                                                             >  1:449575/1‑49 (MQ=255)
                                                                                                  aGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCa                                                                                                             <  2:449575/49‑1 (MQ=255)
                                                                                                  aGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGccc                    >  2:305717/1‑138 (MQ=255)
                                                                                                       aGCCTGTTTCAGAGATACCGGGATCGGCCAACTCa                                                                                                                      <  1:586922/35‑1 (MQ=255)
                                                                                                       aGCCTGTTTCAGAGATACCGGGATCGGCCAACTCa                                                                                                                      >  2:586922/1‑35 (MQ=255)
                                                                                                                   agaTACCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGCCCGATTTCATAGTTCATCCt  >  1:55519/1‑139 (MQ=255)
                                                                                                                       aCCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCt                                            <  2:351637/93‑1 (MQ=255)
                                                                                                                       aCCGGGGTCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCt                                            >  1:351637/1‑93 (MQ=255)
                                                                                                                             |                                                                                                                                
TAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGCCCGATTTCATAGTTCATCCT  >  NZ_CP009273/1713616‑1713869

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