Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,148,473 C→A 9.4% S137S (TCC→TCA) ‡ msuE → FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,148,4730CA 9.4% 121.0 / 2.3 40S137S (TCC→TCA) ‡msuEFMN reductase
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base A (2/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGG  >  NC_002947/3148337‑3148621
                                                                                                                                        |                                                                                                                                                    
tGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCCCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTg                                                                                                                                          <  1:283499/149‑1 (MQ=255)
     aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTc                                                                                                                                       >  4:173378/1‑147 (MQ=255)
     aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgg                                                                                                                                              <  4:93798/140‑1 (MQ=255)
     aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgg                                                                                                                                              >  3:93798/1‑140 (MQ=255)
      tCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAg                                                                                                                                    <  7:331697/149‑1 (MQ=255)
                   cGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                                                                                        <  5:44419/148‑1 (MQ=255)
                        gCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCAc                                                                                                                                 <  8:96328/134‑1 (MQ=255)
                        gCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCAc                                                                                                                                 >  7:96328/1‑134 (MQ=255)
                         ccTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                                                                                        <  2:305458/142‑1 (MQ=255)
                         ccTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                                                                                        >  1:305458/1‑142 (MQ=255)
                          cTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGa                                                                                                                >  7:199726/1‑149 (MQ=255)
                           tGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                                                                               >  1:186292/1‑149 (MQ=255)
                             gCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                                                                               <  2:186292/147‑1 (MQ=255)
                              ccGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCAc                                                                                                            >  2:263920/1‑149 (MQ=255)
                                ggTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCa                                                                                                          <  3:35367/149‑1 (MQ=255)
                                     gCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgggc                                                                                                                                            <  7:234427/110‑3 (MQ=38)
                                     gCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgggc                                                                                                                                            >  8:234427/1‑108 (MQ=38)
                                           gcgcGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAgg                                                                                               >  7:113616/1‑149 (MQ=255)
                                                          ggCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtggt                                                                                  <  8:113616/147‑1 (MQ=255)
                                                          ggCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtggt                                                                                  <  1:263920/147‑1 (MQ=255)
                                                               aCCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgcc                                                                           >  7:322653/1‑149 (MQ=255)
                                                                 cTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgccgc                                                                         >  2:153520/1‑149 (MQ=255)
                                                                  tGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgccgc                                                                         >  6:346612/1‑148 (MQ=255)
                                                                      aCCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCg                                                                    <  1:153520/149‑1 (MQ=255)
                                                                                            cAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCtt                                              >  1:295238/1‑149 (MQ=255)
                                                                                              aGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTtga                                            <  2:295238/149‑1 (MQ=255)
                                                                                                 cGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTtgatga                                         <  3:173378/149‑1 (MQ=255)
                                                                                                       gCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCggtg                                                                 >  8:65245/1‑119 (MQ=255)
                                                                                                       gCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCggtg                                                                 <  7:65245/119‑1 (MQ=255)
                                                                                                       gCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                                    >  3:249536/1‑148 (MQ=255)
                                                                                                        cTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                                    <  8:199726/147‑1 (MQ=255)
                                                                                                         tACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCAc                                 >  5:272225/1‑149 (MQ=255)
                                                                                                              gTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTagcgagc                                                        <  8:132704/121‑1 (MQ=255)
                                                                                                              gTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGtt                            <  4:249536/149‑1 (MQ=255)
                                                                                                                 ggCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgt                                                                                     >  2:24179/1‑89 (MQ=255)
                                                                                                                 ggCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTgc                          <  5:346612/148‑1 (MQ=255)
                                                                                                                       ccGGTGGATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCt                   <  3:180872/149‑1 (MQ=255)
                                                                                                                            gtATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                                    <  2:18470/127‑1 (MQ=255)
                                                                                                                            gtATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                                    >  1:18470/1‑127 (MQ=255)
                                                                                                                                  cgcgTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTttct        >  8:28448/1‑149 (MQ=255)
                                                                                                                                        cTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAgg  <  7:28448/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
TGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGG  >  NC_002947/3148337‑3148621

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