Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,911,489 C→T 8.7% R267C (CGC→TGC)  PP_3451 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,911,4890CT 8.7% 127.6 / 2.4 46R267C (CGC→TGC) PP_3451hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (2/2);  total (20/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 = 3.94e-01

TCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCCGCGAAAGGTCCAAGCC  >  NC_002947/3911341‑3911624
                                                                                                                                                    |                                                                                                                                       
tCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTgc                                                                                                                                         <  6:277710/149‑1 (MQ=255)
tCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTgc                                                                                                                                         >  6:142014/1‑149 (MQ=255)
   gggTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTgcgc                                                                                                                                       >  6:370828/1‑148 (MQ=255)
       tACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCg                                                                                                                                  >  5:130802/1‑149 (MQ=255)
               ggCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTggg                                                                                                                          >  2:124448/1‑149 (MQ=255)
                gCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGc                                                                                                                         >  1:164318/1‑149 (MQ=255)
                  gTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGccc                                                                                                                       >  6:11329/1‑149 (MQ=255)
                    tgccgtgccgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCgg                                                                                                                     <  3:79769/149‑1 (MQ=255)
                      ccgtgccgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCggg                                                                                                                    <  5:142014/148‑1 (MQ=255)
                      ccgtgccgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCggg                                                                                                                    <  3:181583/148‑1 (MQ=255)
                       cgtgccgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGt                                                                                                                  <  3:345530/149‑1 (MQ=255)
                        gtgccgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGCGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTa                                                                                                                 <  3:222784/149‑1 (MQ=255)
                            cgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCATGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCtt                                                                                                             <  4:203682/149‑1 (MQ=255)
                            cgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCtt                                                                                                             <  2:343529/149‑1 (MQ=255)
                             gggtgggtGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTg                                                                                                            <  5:108154/149‑1 (MQ=255)
                              ggtgggtgTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGc                                                                                                           <  7:100539/149‑1 (MQ=255)
                                 gggtgTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCca                                                                                                         >  1:271526/1‑148 (MQ=255)
                                            gcgAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAAtt                                                                                             <  7:199213/149‑1 (MQ=255)
                                              gAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTAc                                                                                           <  5:370828/149‑1 (MQ=255)
                                              gAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGT‑CCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCGGGGGTATCTTGCCACAGCAATGAATTAc                                                                                           >  6:52256/1‑148 (MQ=255)
                                               aGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTAcc                                                                                          >  2:198308/1‑149 (MQ=255)
                                                  cGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAgg                                                                                       <  2:164318/149‑1 (MQ=255)
                                                      ccGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGggcc                                                                                    <  5:11329/148‑1 (MQ=255)
                                                        gCCAGGTGAGCGCGCCTGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGggccggc                                                                                 <  1:124448/149‑1 (MQ=255)
                                                         ccAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGggccggcc                                                                                >  5:54809/1‑149 (MQ=255)
                                                                    cgcCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGAt                                                                     <  6:130802/149‑1 (MQ=255)
                                                                           gCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCca                                                              >  2:67592/1‑149 (MQ=255)
                                                                                      gggCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTcccccc                                                   >  6:320646/1‑148 (MQ=255)
                                                                                       ggCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGTGCTACCGGcc                                                                                                                               >  1:39948/1‑72 (MQ=255)
                                                                                       ggCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGTGCTACCGGcc                                                                                                                               >  7:299076/1‑72 (MQ=255)
                                                                                       ggCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGTGCTACCGGcc                                                                                                                               <  8:299076/72‑1 (MQ=255)
                                                                                       ggCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGTGCTACCGGcc                                                                                                                               <  2:39948/72‑1 (MQ=255)
                                                                                             gCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTa                                            <  2:160311/149‑1 (MQ=255)
                                                                                                 gccgAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGt                                        >  3:275974/1‑149 (MQ=255)
                                                                                                     aGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCt                                    <  6:54809/149‑1 (MQ=255)
                                                                                                         cAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAg                                >  2:333506/1‑149 (MQ=255)
                                                                                                         cAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCAAACCTTACGGTCCCTGTAg                                >  8:222811/1‑149 (MQ=255)
                                                                                                                ggtAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCgg                          <  2:271526/148‑1 (MQ=255)
                                                                                                                         tGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCcgc                <  4:275974/149‑1 (MQ=255)
                                                                                                                           ccATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTAAACCcgctc              <  1:198308/149‑3 (MQ=255)
                                                                                                                            cATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCCGCGa              >  4:32671/1‑148 (MQ=255)
                                                                                                                              tGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCCGCGAAAg           <  1:333506/149‑1 (MQ=255)
                                                                                                                                 cGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGCGGTATCTTGCCACAGCATTGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGCAGCGGGTTTACCCGCGAAAGGTc        >  2:301728/1‑149 (MQ=255)
                                                                                                                                  gAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCca                                                   <  4:163467/105‑1 (MQ=255)
                                                                                                                                  gAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCca                                                   >  3:163467/1‑105 (MQ=255)
                                                                                                                                  gAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTAccc                  <  5:248075/138‑1 (MQ=255)
                                                                                                                                  gAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTAccc                  >  6:248075/1‑138 (MQ=255)
                                                                                                                                    ggCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCCGCGAAAGGTCCaa     <  1:67592/149‑1 (MQ=255)
                                                                                                                                       aGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCCGCGAAAGGTCCAAGcc  >  7:21544/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
TCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCACAGGTCCCCCACACCGTACGGTCCCTGTAGGAGCGGGTTTACCCGCGAAAGGTCCAAGCC  >  NC_002947/3911341‑3911624

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