Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,878,241 C→T 8.7% L252L (CTG→CTA)  PP_5152 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,878,2410CT 8.7% 128.2 / 2.7 46L252L (CTG→CTA) PP_5152hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (2/2);  total (23/23)
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.92e-01

ATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCATCGGTGTTC  >  NC_002947/5878094‑5878385
                                                                                                                                                   |                                                                                                                                                
aTAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGATTGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGtccagctcc                                                                                                                                                  >  7:70195/1‑148 (MQ=255)
   gAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCTATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCcagca                                                                                                                                              <  7:162956/149‑1 (MQ=255)
   gAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCcagca                                                                                                                                              <  2:116322/149‑1 (MQ=255)
          gcCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGtt                                                                                                                                       <  6:167329/149‑1 (MQ=255)
          gcCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGtt                                                                                                                                       >  2:306946/1‑149 (MQ=255)
           ccGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCTAGCAGCATGTTg                                                                                                                                      >  3:295902/1‑149 (MQ=255)
           ccGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCTAGCAGCATGTTg                                                                                                                                      >  7:229964/1‑149 (MQ=255)
           ccGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTg                                                                                                                                      <  1:144804/149‑1 (MQ=255)
               cTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGt                                                                                                                                  <  6:298222/149‑1 (MQ=255)
                   cTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGc                                                                                                                              <  5:196498/149‑1 (MQ=255)
                    tGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTc                                                                                                                                 <  2:221807/145‑1 (MQ=255)
                    tGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTc                                                                                                                                 >  1:221807/1‑145 (MQ=255)
                       ccATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGa                                                                                                                           >  7:239274/1‑148 (MQ=255)
                                 gTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGc                                                                                                                 >  3:80881/1‑148 (MQ=255)
                                  ttGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCaa                                                                                                               <  5:119144/149‑1 (MQ=255)
                                     ccccTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTAAATCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGATGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAgag                                                                                                            >  6:33647/1‑147 (MQ=255)
                                       ccTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGACCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGc                                                                                                          >  5:276140/1‑149 (MQ=255)
                                          gtcgtcCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTGATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGACTCGGCCGATACAACCAGCAAGTTGCCAg                                                                                                       >  4:40828/1‑149 (MQ=255)
                                           tcgtcCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCa                                                                                                        <  2:1635/147‑1 (MQ=255)
                                             gtcCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTAc                                                                                                    >  3:4269/1‑149 (MQ=255)
                                             gtcCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTAc                                                                                                    >  2:191173/1‑149 (MQ=255)
                                                 aTGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGtccagctcca                                                                                                                                                 <  6:259298/100‑1 (MQ=255)
                                                 aTGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGtccagctcca                                                                                                                                                 >  5:259298/1‑100 (MQ=255)
                                                 aTGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCgg                                                                                                <  4:4269/149‑1 (MQ=255)
                                                            ttCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACAcc                                                                                       >  1:295514/1‑147 (MQ=255)
                                                                  ccGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCAGGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAgcg                                                                                  >  1:91616/1‑146 (MQ=255)
                                                                   cgcgGAATGGTGCCTTCAAGCACGTTTAAGCCGCCCAGCAGCGTCATGAATATGACTTGCTGGCCCTGATCTTCCAGCTCCAGCAGCATGTTGGATTAGGCCGATATCACCAGCAAGTGGACAGTACGCGGGTCGACACCAAGCGCAGa                                                                              >  8:114334/1‑149 (MQ=255)
                                                                             tGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTgcg                                                                     >  3:386042/1‑148 (MQ=255)
                                                                                    aaGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGa                                                             <  6:335532/149‑1 (MQ=255)
                                                                                    aaGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGa                                                             <  8:320718/149‑1 (MQ=255)
                                                                                       cccGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGCCAGGTTGCGCAAGTCGAGtt                                                          >  8:264877/1‑149 (MQ=255)
                                                                                        ccGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTc                                                         <  1:191173/149‑1 (MQ=255)
                                                                                           ttGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATc                                                      <  2:295514/149‑1 (MQ=255)
                                                                                           ttGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATc                                                      <  4:386042/149‑1 (MQ=255)
                                                                                                gccgccCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATTGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAg                                                 >  5:59164/1‑149 (MQ=255)
                                                                                                gccgccCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAg                                                 >  3:105446/1‑149 (MQ=255)
                                                                                                        gcagcGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGgc                                         >  2:282122/1‑149 (MQ=255)
                                                                                                           gcGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGgctgc                                      <  6:276140/149‑1 (MQ=255)
                                                                                                             gTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCTAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCtt                                    <  8:229964/149‑1 (MQ=255)
                                                                                                             gTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCTAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCtt                                    <  4:295902/149‑1 (MQ=255)
                                                                                                              tCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACgg                                          <  3:301111/142‑1 (MQ=255)
                                                                                                              tCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACgg                                          >  4:301111/1‑142 (MQ=255)
                                                                                                              tCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTtc                                   <  3:40828/149‑1 (MQ=255)
                                                                                                                aTGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGtt                                                          >  7:275252/1‑124 (MQ=255)
                                                                                                                aTGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGtt                                                          <  8:275252/124‑1 (MQ=255)
                                                                                                                      cTGACTTGCTTGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCtt                           <  8:239274/149‑1 (MQ=255)
                                                                                                                                   ccTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCgg              <  1:306946/149‑1 (MQ=255)
                                                                                                                                    cTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGCCAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGc             >  8:224036/1‑149 (MQ=255)
                                                                                                                                      gCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCa           >  6:343114/1‑149 (MQ=255)
                                                                                                                                       cTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCAt          <  2:91616/149‑1 (MQ=255)
                                                                                                                                         cGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCATCg        <  7:264877/149‑1 (MQ=255)
                                                                                                                                             cagctccagcAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCATCGgtgt    >  7:293729/1‑149 (MQ=255)
                                                                                                                                                 tccagcAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCATCGGTGTTc  >  5:160617/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATGGTGCCTTCAAGCCCGTTGAAGCCGCCCAGCAGCGTCATGAAGCTGACTTGCTGGCCCTGCTCGTCCAGCTCCAGCAGCATGTTGGAGTCGGCCGATAGCACCAGCAAGTGGCCAGTACGCGGGTCGACACCCAGCGCAGACAGGTTGCGCAAGTCGAGTTCATCGCTAGGCAACGGCTGCTTCTCACCCTTGAGCGGGCTGCGGCCATCGGTGTTC  >  NC_002947/5878094‑5878385

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