Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,114 G→T 13.0% T423N (ACT→AAT)  trmF ← tRNA uridine 5‑carboxymethylaminomethyl modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,1140GT13.0% 130.6 / 9.7 46T423N (ACT→AAT) trmFtRNA uridine 5‑carboxymethylaminomethyl modification protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base T (3/3);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTG  >  NC_002947/2967‑3260
                                                                                                                                                    |                                                                                                                                                  
cagcgcTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAg                                                                                                                                                    >  2:195393/3‑148 (MQ=255)
        ggCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCaa                                                                                                                                             <  1:195393/147‑1 (MQ=255)
         gcgcAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGt                                                                                                                                          <  5:68546/149‑1 (MQ=255)
              gAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCa                                                                                                                                     <  6:25697/149‑1 (MQ=255)
                aTGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAgg                                                                                                                                   <  1:293529/149‑1 (MQ=255)
                 tGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACACTGCGGTACGGCTCCTGAGTACCCAAGGTAATCAgg                                                                                                                                   <  1:54789/149‑1 (MQ=255)
                 tGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGATTACCTTGGGTACTCAgga                                                                                                                                  <  1:139370/149‑2 (MQ=255)
                 tGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGgt                                                                                                                                  >  5:124666/1‑149 (MQ=255)
                          gcgcTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTccacca                                                                                                                         <  3:128534/149‑1 (MQ=255)
                                gTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGc                                                                                                                   >  5:43684/1‑149 (MQ=255)
                                  cATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCa                                                                                                                                     <  8:90943/129‑1 (MQ=255)
                                  cATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCa                                                                                                                                     >  7:90943/1‑129 (MQ=255)
                                   aTCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGTTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTcca                                                                                                                            <  6:2904/137‑1 (MQ=255)
                                   aTCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTcca                                                                                                                            >  5:2904/1‑137 (MQ=255)
                                      gATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGt                                                                                                             <  5:333585/149‑1 (MQ=255)
                                           ggCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTAt                                                                                                           <  6:313045/146‑1 (MQ=255)
                                           ggCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTAt                                                                                                           >  5:313045/1‑146 (MQ=255)
                                           ggCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGcc                                                                                                        <  8:115920/149‑1 (MQ=255)
                                             cccAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCtc                                                                                                      <  3:297399/149‑1 (MQ=255)
                                                aGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCtcgtc                                                                                                   >  3:131455/1‑149 (MQ=255)
                                                     aCGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggc                                                                                              >  8:148824/1‑149 (MQ=255)
                                                       ggCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgc                                                                                            >  3:89140/1‑149 (MQ=255)
                                                        gCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                           <  3:73480/149‑1 (MQ=255)
                                                          ccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                           <  4:165550/147‑1 (MQ=255)
                                                          ccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                           >  3:165550/1‑147 (MQ=255)
                                                                 ggTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGc                                                                                  <  6:43684/149‑1 (MQ=255)
                                                                       gcgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCAc                                                                            >  6:216190/1‑149 (MQ=255)
                                                                         gcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACgg                                                                          <  6:124666/149‑1 (MQ=255)
                                                                                            gcgcAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCtcgtcg                                                                                                  >  4:95995/1‑106 (MQ=255)
                                                                                             cacAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTgcggc                                                      >  5:313633/3‑149 (MQ=255)
                                                                                                 gTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTa                                                  <  4:131455/149‑1 (MQ=255)
                                                                                                      aGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTccc                                             >  1:24443/1‑149 (MQ=255)
                                                                                                       gCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCg                                            >  8:48240/1‑149 (MQ=255)
                                                                                                                 ggCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAgg                                   <  5:216190/148‑1 (MQ=255)
                                                                                                                 ggCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGc                                  >  6:292801/1‑149 (MQ=255)
                                                                                                                                cA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGct                    >  8:303628/1‑148 (MQ=255)
                                                                                                                                 a‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGaa                                                             >  6:315105/1‑106 (MQ=255)
                                                                                                                                 a‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGaa                                                             <  5:315105/106‑1 (MQ=255)
                                                                                                                                   tGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGcttct                 <  7:148824/149‑1 (MQ=255)
                                                                                                                                     cGGTACGGCTCCTGAGTACCCAATGTAATCAGGTCGTCGACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTGACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCttcttc               <  4:89140/149‑1 (MQ=255)
                                                                                                                                      ggTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCg              >  2:80294/1‑149 (MQ=255)
                                                                                                                                             gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                             <  8:89708/125‑1 (MQ=255)
                                                                                                                                             gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                             >  7:84559/3‑127 (MQ=255)
                                                                                                                                             gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                             <  8:84559/125‑1 (MQ=255)
                                                                                                                                             gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                             >  7:89708/3‑127 (MQ=255)
                                                                                                                                             gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCgg       >  1:221708/3‑149 (MQ=255)
                                                                                                                                                  gAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGgtggtg  >  3:221757/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTG  >  NC_002947/2967‑3260

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑