Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,437,066 T→C D608D (GAT→GAC mnmC → bifunctional tRNA (5‑methylaminomethyl‑2‑thiouridine)(34)‑methyltransferase MnmD/FAD‑dependent 5‑carboxymethylaminomethyl‑2‑thiouridine(34) oxidoreductase MnmC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,437,0660TC92.1% 110.4 / 0.0 38D608D (GAT→GACmnmCbifunctional tRNA (5‑methylaminomethyl‑2‑thiouridine)(34)‑methyltransferase MnmD/FAD‑dependent 5‑carboxymethylaminomethyl‑2‑thiouridine(34) oxidoreductase MnmC
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (20/15);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

CGCTGCGGTGTGCGTTGTGCCACCCGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAA  >  NZ_CP009273/2436929‑2437197
                                                                                                                                         |                                                                                                                                   
cgcTGCGGTGTGCGTTGTGCCACCCGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGAt                                                                                                                                     >  1:13542/1‑138 (MQ=255)
cgcTGCGGTGTGCGTTGTGCCACCCGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGAt                                                                                                                                     <  2:13542/138‑1 (MQ=255)
                 tgCCACCCGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCt                                                                                                                   >  2:746450/1‑139 (MQ=255)
                      cccGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTTTGTTTGAGGGTTTAgg                                                                                                              >  2:8432/1‑139 (MQ=255)
                      cccGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAgg                                                                                                              >  1:448005/1‑139 (MQ=255)
                           gATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGGtctc                                                                                                         >  1:863349/1‑139 (MQ=255)
                             tcatcTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTcgc                                                                                                       <  1:8432/139‑1 (MQ=255)
                                       aTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGc                                                                                                                       >  1:326240/1‑113 (MQ=255)
                                       aTGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGc                                                                                                                       <  2:326240/113‑1 (MQ=255)
                                         ggTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCt                                                                                           >  1:273572/1‑139 (MQ=255)
                                             ggCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGcc                                                                                        <  2:748361/138‑1 (MQ=255)
                                              gCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGcccc                                                                                      >  2:833600/1‑139 (MQ=255)
                                                 aTGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCt                                                                                   <  1:110189/139‑1 (MQ=255)
                                                 aTGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCt                                                                                   <  1:833600/139‑1 (MQ=255)
                                                  tGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCtg                                                                                  <  2:863349/139‑1 (MQ=255)
                                                                     acacTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTggcgg                                                               <  2:273572/139‑1 (MQ=255)
                                                                       acTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGgcggcg                                                             <  1:468739/139‑1 (MQ=255)
                                                                            tGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTtct                                                                                                          <  2:413665/89‑1 (MQ=255)
                                                                            tGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTtct                                                                                                          >  1:413665/1‑89 (MQ=255)
                                                                            tGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGc                                                                                    >  1:1713/1‑111 (MQ=255)
                                                                            tGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGc                                                                                    <  2:1713/111‑1 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGAt                                                                                                                                     >  1:587490/1‑55 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGAt                                                                                                                                     <  2:587490/55‑1 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCt                                                                                   >  1:216281/1‑105 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCt                                                                                   <  2:216281/105‑1 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGcgac                                                 >  1:112227/1‑139 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGcgac                                                 >  2:652679/1‑139 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGcgac                                                 >  2:263957/1‑139 (MQ=255)
                                                                                   tGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGGTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGcgac                                                 >  1:766601/1‑139 (MQ=255)
                                                                                           tGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTcgc                                                                                                       <  1:849825/77‑1 (MQ=255)
                                                                                           tGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTcgc                                                                                                       >  2:849825/1‑77 (MQ=255)
                                                                                                  gCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatg                                  >  1:731499/1‑139 (MQ=255)
                                                                                                    aGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatgga                                <  1:393704/139‑1 (MQ=255)
                                                                                                             aGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGtt                                                                                             <  2:479661/69‑1 (MQ=255)
                                                                                                             aGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGtt                                                                                             >  1:479661/1‑69 (MQ=255)
                                                                                                             aGGCGGTAAGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTAc                       >  1:35577/1‑139 (MQ=255)
                                                                                                                    aaGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatgga                                <  2:653473/123‑1 (MQ=255)
                                                                                                                    aaGCGCGCCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCgatgga                                >  1:653473/1‑123 (MQ=255)
                                                                                                                         cgcCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGt            >  2:230051/1‑138 (MQ=255)
                                                                                                                         cgcCGGTTTTTGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGt            >  1:473281/1‑138 (MQ=255)
                                                                                                                                   tGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGaa  <  2:35577/138‑1 (MQ=255)
                                                                                                                                   tGACGACCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGaa  <  2:473281/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
CGCTGCGGTGTGCGTTGTGCCACCCGCGATCATCTGCCAATGGTAGGCAATGTTCCCGATTATGAGGCAACACTCGTGGAATATGCGTCGTTGGCGGAGCAGAAAGATGAGGCGGTAAGCGCGCCGGTTTTTGACGATCTCTTTATGTTTGCGGCTTTAGGTTCTCGCGGTTTGTGTTCTGCCCCGCTGTGTGCCGAGATTCTGGCGGCGCAGATGAGCGACGAACCGATTCCGATGGATGCCAGTACGCTGGCGGCGTTAAACCCGAA  >  NZ_CP009273/2436929‑2437197

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