Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,276,105 C→T 8.0% S40L (TCG→TTG)  PP_2874 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,276,1050CT 8.0% 141.0 / 2.6 50S40L (TCG→TTG) PP_2874hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

TTCCAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGAT  >  NC_002947/3275967‑3276243
                                                                                                                                          |                                                                                                                                          
ttCTAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTTgcgc                                                                                                                                        <  4:213371/143‑1 (MQ=255)
ttCTAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTTgcgc                                                                                                                                        <  4:209140/143‑1 (MQ=255)
ttCAAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTTgcgc                                                                                                                                        >  3:209140/1‑143 (MQ=255)
ttCAAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTTgcgc                                                                                                                                        >  3:213371/1‑143 (MQ=255)
  ccAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTggggg                                                                                                                                 <  8:288523/148‑1 (MQ=255)
  ccAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGc                                                                                                                                >  1:16027/1‑149 (MQ=255)
    aGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGa                                                                                                                              <  5:318137/149‑1 (MQ=255)
        aCGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACg                                                                                                                            >  2:231218/1‑147 (MQ=255)
             ggaTTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCGGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGt                                                                                                                     <  6:98715/149‑1 (MQ=255)
             ggaTTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGt                                                                                                                     <  6:244795/149‑1 (MQ=255)
                ttGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACACAGCTGTTTc                                                                                                                  >  3:7021/1‑149 (MQ=255)
                    aTGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTgg                                                                                                               <  8:18112/148‑1 (MQ=255)
                       cTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGtt                                                                                                           >  2:208490/1‑149 (MQ=255)
                        tCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTg                                                                                                          >  3:11835/1‑149 (MQ=255)
                          aTGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTa                                                                                                        >  1:2172/1‑149 (MQ=255)
                                 tCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACATCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCggg                                                                                                 <  5:286527/149‑1 (MQ=255)
                                    gggAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTggggg                                                                                                                                 <  1:133293/114‑1 (MQ=255)
                                    gggAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTggggg                                                                                                                                 >  2:133293/1‑114 (MQ=255)
                                         gaaACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATgggcg                                                                                         >  6:339580/1‑149 (MQ=255)
                                          aaaCCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATgggcgg                                                                                        <  1:208490/149‑1 (MQ=255)
                                            aCCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATgggcgggc                                                                                      >  7:152247/1‑149 (MQ=255)
                                                    gCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGt                                                                                                         <  2:330069/122‑1 (MQ=255)
                                                    gCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGt                                                                                                         >  1:330069/1‑122 (MQ=255)
                                                      aGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGaa                                                                            <  4:11835/149‑1 (MQ=255)
                                                      aGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGaa                                                                            <  1:231218/149‑1 (MQ=255)
                                                        tGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAaga                                                                          <  8:152247/149‑1 (MQ=255)
                                                        tGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAaga                                                                          <  3:228484/149‑1 (MQ=255)
                                                              gACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTa                                                                    <  5:339580/149‑1 (MQ=255)
                                                                   acaGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGTCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCcg                                                               <  5:180548/149‑1 (MQ=255)
                                                                   acaGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCcg                                                               >  1:175958/1‑149 (MQ=255)
                                                                             aGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGCCGTGCAGGTGCGCTTGTCGTGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCggg                                                     >  1:239133/1‑149 (MQ=255)
                                                                                     ccACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCa                                             >  3:50211/1‑149 (MQ=255)
                                                                                       acgacgAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCAt                                            >  1:294792/1‑148 (MQ=255)
                                                                                         gacgaAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGAc                                         <  4:50211/149‑1 (MQ=255)
                                                                                         gacgaAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGAc                                         >  5:229618/1‑149 (MQ=255)
                                                                                            gaAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCg                                      >  5:220829/1‑149 (MQ=255)
                                                                                                         aCGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTgg                         <  2:2172/149‑1 (MQ=255)
                                                                                                              accacACTGGTCGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGgtgt                    >  4:163133/1‑149 (MQ=255)
                                                                                                              accacACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCga                                                  <  7:151560/119‑1 (MQ=255)
                                                                                                              accacACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCga                                                  >  8:151560/1‑119 (MQ=255)
                                                                                                              accacACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTggg                        <  1:157327/145‑1 (MQ=255)
                                                                                                              accacACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTggg                        >  2:157327/1‑145 (MQ=255)
                                                                                                               ccacACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTc                   <  2:175958/149‑1 (MQ=255)
                                                                                                                cacaCTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTc                   >  2:87640/1‑148 (MQ=255)
                                                                                                                             cAGGTGCGCTTGTCGTGGCTGGGGGCTACGCAGCTGTTTCTGTGGTTGTAAGCGGGGGATGGGGGGGGGTTGGTGGTAGAGTATTACCCCGCGCTGCCGGGCGAGTCCATGACCCGTTGGCTGGTGTGGTGTGTGGACCGCGCCCggg      <  2:239133/148‑1 (MQ=255)
                                                                                                                              aGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGgt                             >  5:86465/1‑124 (MQ=255)
                                                                                                                              aGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGgt                             <  6:86465/124‑1 (MQ=255)
                                                                                                                              aGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCg    <  2:202237/149‑1 (MQ=255)
                                                                                                                                   cgcTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGAt  >  2:119491/1‑146 (MQ=255)
                                                                                                                                   cgcTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGAt  <  1:119491/146‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TTCCAGCAACGGAGGATTGAATGCTCATGCGAGTCAGGGAAGAAACCTATTGGCAGTGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGAT  >  NC_002947/3275967‑3276243

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