Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,244,522 T→C 12.9% V23A (GTG→GCG)  dusC → tRNA‑dihydrouridine synthase DusC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,244,5220TC12.9% 73.8 / 3.6 32V23A (GTG→GCG) dusCtRNA‑dihydrouridine synthase DusC
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base C (2/2);  total (15/17)
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.87e-01

GTCCGGCGCTGATCGGCGCCTGAGTGGCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTGCAATTGCTGGG  >  NC_002947/2244381‑2244666
                                                                                                                                             |                                                                                                                                                 
ctcaggCGCCGATCAGCGCCTGAGTGGCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCAGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGAGGG‑Cgg                                                                                                                                            >  4:262669/5‑148 (MQ=255)
     gcgcTGATCGGCGCCTGAGTGGCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTagc                                                                                                                                        >  2:53191/1‑145 (MQ=255)
        cTGATCGGCGCCTGAGTGGCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTg                                                                                                                                   <  2:263864/149‑1 (MQ=255)
           atgggCGCCTGAGTGGCAGGCATTTTTTGTTAAAGTGCCGGCTTTTTCGGGTGATAGCCAGACATGCAAATTGCCCTGGCCCCCTTGGAGTGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑GGGTATCGACTGGt                                                                                                                                 <  8:251089/145‑1 (MQ=255)
           aTCGGCGCCTGAGTGGCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGt                                                                                                                                 <  6:151116/148‑1 (MQ=255)
                          gCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCa                                                                                                                  <  8:173847/148‑1 (MQ=255)
                                   ttttGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCgtgt                                                                                                         <  7:223298/148‑1 (MQ=255)
                                          aaaCTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGc                                                                                                 >  3:137254/1‑149 (MQ=255)
                                             cTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGAATGGTGAGTCACAGAGTTCATTCGCGTGTGCCACCGCCTg                                                                                              >  8:116807/1‑149 (MQ=255)
                                                    cTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGGCGACAACATCCTGCGCGACGTATTGCCTCGTGTGGGACGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCg                                                                                        <  1:249757/149‑1 (MQ=255)
                                                             cTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGt                                                                              <  4:137254/149‑1 (MQ=255)
                                                                agaCCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTc                                                                            <  6:100837/148‑1 (MQ=255)
                                                                    cccgcaaTGCAAATTGCCCTGGCCCCCATGGAGCGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACaa                                                                       <  8:116416/143‑1 (MQ=255)
                                                                       gACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTcgc                                                                                                             <  5:150554/108‑1 (MQ=255)
                                                                       gACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTcgc                                                                                                             >  6:150554/1‑108 (MQ=255)
                                                                                    cccTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTg                                                       >  2:289602/1‑149 (MQ=255)
                                                                                    cccTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTg                                                       >  3:256969/1‑149 (MQ=255)
                                                                                         gCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTc                                                   <  3:313190/148‑1 (MQ=255)
                                                                                           ccccATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCgtgt                                                                                                         >  8:240390/1‑92 (MQ=255)
                                                                                           ccccATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCgtgt                                                                                                         <  7:240390/92‑1 (MQ=255)
                                                                                           ccccATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAgg                                                <  6:120092/149‑1 (MQ=255)
                                                                                           ccccATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAgg                                                >  6:310153/1‑149 (MQ=255)
                                                                                            cccATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCTACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGgcgc                                                               <  8:188774/133‑1 (MQ=255)
                                                                                            cccATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCTACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGgcgc                                                               >  7:188774/1‑133 (MQ=255)
                                                                                               aTGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGgcgc                                            >  2:157522/1‑149 (MQ=255)
                                                                                                    ggggCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTc                                                                                <  3:50148/108‑1 (MQ=255)
                                                                                                    ggggCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTAGCCGGCGTCCTc                                                                                >  4:50148/1‑108 (MQ=255)
                                                                                                                            cgACGTATTGACTCGTGCGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATgc                  >  8:36361/1‑146 (MQ=255)
                                                                                                                            cgACGTATTGACTCGTGCGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATgc                  >  8:41380/1‑146 (MQ=255)
                                                                                                                            cgACGTATTGACTCGTGCGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATgc                  <  7:41380/146‑1 (MQ=255)
                                                                                                                            cgACGTATTGACTCGTGCGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATgc                  <  7:36361/146‑1 (MQ=255)
                                                                                                                                     gACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTGCAATTGc      >  6:202840/1‑149 (MQ=255)
                                                                                                                                          gtgtGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTg             >  8:183849/1‑137 (MQ=255)
                                                                                                                                          gtgtGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTg             <  7:183849/137‑1 (MQ=255)
                                                                                                                                          gtgtGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTGCAATTGCTgg   <  5:202840/147‑1 (MQ=255)
                                                                                                                                          gtgtGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGGGTGCGACCGCCGGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTGCAATTGCTggg  <  7:116807/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GTCCGGCGCTGATCGGCGCCTGAGTGGCCGGCATTTTTTGTTAAACTGCCGGCTTTTTCCGCTGAGACCCCGACATGCAAATTGCCCTGGCCCCCATGGAGGGGCTGGTCGACAACATCCTGCGCGACGTATTGACTCGTGTGGG‑CGGTATCGACTGGTGCGTCACCGAGTTCATTCGCGTGTGCGACCGCCTGTTGCCGGCGTCCTCGTTCGACAAGCTGGCGCCTGAACTGCGTCAGGGCGCGCGCACGGCTGCCGGTGTGCCCATGCGCGTGCAATTGCTGGG  >  NC_002947/2244381‑2244666

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