Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,120 A→T 12.2% L421* (TTG→TAG) ‡ 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,1200AT12.2% 130.0 / 8.4 49L421* (TTG→TAG) ‡trmFtRNA uridine 5‑carboxymethylaminomethyl modification protein
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base T (3/3);  total (25/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 = 8.24e-01

GGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTT  >  NC_002947/2975‑3265
                                                                                                                                                  |                                                                                                                                                 
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                                                                                                >  3:165550/1‑147 (MQ=255)
                                                  ccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                                <  4:165550/147‑1 (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                                  >  7:89708/3‑127 (MQ=255)
                                                                                                                                     gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                                  <  8:84559/125‑1 (MQ=255)
                                                                                                                                     gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCgg            >  1:221708/3‑149 (MQ=255)
                                                                                                                                          gAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGgtggtg       >  3:221757/1‑149 (MQ=255)
                                                                                                                                               cccAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt  >  3:371011/1‑149 (MQ=255)
                                                                                                                                               cccAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGt   <  1:80294/148‑1 (MQ=255)
                                                                                                                                                ccAAGGTAATCTGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt  >  1:286216/1‑148 (MQ=255)
                                                                                                                                                ccAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt  >  4:254315/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTT  >  NC_002947/2975‑3265

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: ‑