Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,389,946 C→G 11.3% A23G (GCC→GGC)  PP_1212 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,389,9460CG11.3% 147.8 / 10.0 53A23G (GCC→GGC) PP_1212hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base G (3/3);  total (24/29)
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.77e-01

GCAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGC  >  NC_002947/1389801‑1390085
                                                                                                                                                 |                                                                                                                                           
gcAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcggcgc                                                                                                                                          >  6:277054/1‑149 (MQ=255)
          ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                <  5:277054/149‑1 (MQ=255)
          ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                >  6:107419/1‑149 (MQ=255)
          ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATc                                                                                                                                >  5:275360/1‑149 (MQ=255)
             aCCATATATAATCCCTCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTTCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAt                                                                                                                             <  2:95345/149‑1 (MQ=255)
             aCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCg                                                                                                                               <  5:51164/147‑1 (MQ=255)
             aCCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGa                                                                                                                              <  7:102527/148‑1 (MQ=255)
                    ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGAt                                                                                                                                 >  1:153264/1‑138 (MQ=255)
                    ataATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGAt                                                                                                                                 <  2:153264/138‑1 (MQ=255)
                        tCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTg                                                                                                                  >  3:278830/1‑149 (MQ=255)
                          ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                 >  8:209687/1‑148 (MQ=255)
                                  tgctgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGcc                                                                                                         <  6:286266/148‑1 (MQ=255)
                                     tgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCa                                                                                                     >  2:327374/1‑149 (MQ=255)
                                           ccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGaaa                                                                                               <  7:318756/149‑1 (MQ=255)
                                                agccagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgct                                                                                          <  6:275360/149‑1 (MQ=255)
                                                       gTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTc                                                                                    <  7:209687/148‑1 (MQ=255)
                                                           tGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                               <  4:278830/149‑1 (MQ=255)
                                                           tGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                               <  1:327374/149‑1 (MQ=255)
                                                             gTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                               >  6:128647/1‑147 (MQ=255)
                                                              ttGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                            >  5:60678/1‑149 (MQ=255)
                                                               tGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                           >  4:285140/1‑149 (MQ=255)
                                                                 aCGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                            <  6:276730/146‑1 (MQ=255)
                                                                 aCGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAg                                                                            >  5:276730/1‑146 (MQ=255)
                                                                    agagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                      <  5:107419/149‑1 (MQ=255)
                                                                           cccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAg                                                               <  6:67116/149‑1 (MQ=255)
                                                                            ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                            <  4:181591/87‑1 (MQ=255)
                                                                            ccATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGAtt                                                                                                                            >  3:181591/1‑87 (MQ=255)
                                                                                  cccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGATGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAAc                                                        >  8:22969/1‑149 (MQ=255)
                                                                                   ccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                               >  5:142132/1‑93 (MQ=255)
                                                                                   ccTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCa                                                                                                               <  6:142132/93‑1 (MQ=255)
                                                                                        aTGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTg                                                   <  7:121508/148‑1 (MQ=255)
                                                                                        aTGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGAAGAAGATCAGCGCCGAGACGCAGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAACCTGCAGTCGGTGCCAGCCTTTCGCCTCAGCGGTAACGGTTgg                                                  >  8:143003/1‑149 (MQ=255)
                                                                                         tGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTgg                                                  <  3:8636/148‑1 (MQ=255)
                                                                                            cTATCAATGTGCGTCCTGCGAGAACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGACAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                               >  2:231585/1‑148 (MQ=255)
                                                                                                 aaTGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                         >  5:178843/1‑149 (MQ=255)
                                                                                                  aTGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCAGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                         >  1:265885/1‑148 (MQ=255)
                                                                                                     tgCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGAc                                      <  7:22969/148‑1 (MQ=255)
                                                                                                       cGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                   <  6:60678/149‑1 (MQ=255)
                                                                                                         tCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCaa                                 <  6:178843/149‑1 (MQ=255)
                                                                                                         tCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCaa                                 <  3:285140/149‑1 (MQ=255)
                                                                                                                    ccGCATGGAAGTGCTGCAGAAGATCAGAGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGTGCg                       >  5:339270/1‑148 (MQ=255)
                                                                                                                           gAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGa                <  7:143003/148‑1 (MQ=255)
                                                                                                                              gTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCt            <  2:265885/149‑1 (MQ=255)
                                                                                                                               tgctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTg           <  6:280965/149‑1 (MQ=255)
                                                                                                                                    cAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGACAGGCTTTCGCCTGAGCTGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCa        <  1:231585/147‑1 (MQ=255)
                                                                                                                                      gaagaTCAGCGCCGCGCGGCTGATCGATTGCCCGGCCTGTCAGGCGCGGGCACTGAAAAAGCTGCTGTCGGTGCCAGGGTTTCGCGTGAGCGGTAGCGGTGGGTACGAATCCGAGTTCTAGGCGGGGGGGGAAAAGAATCTGGGAggcg    >  2:299950/1‑149 (MQ=255)
                                                                                                                                      gaagaTCAGAGCCGCGCCGCTGATCGATTGCCCGGCCAGTCACGCGCCGGCACTGAAAAAGCAGCTGTCGGTGACAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAACGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAgccg    >  2:284625/1‑149 (MQ=255)
                                                                                                                                         gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                          >  7:300685/1‑124 (MQ=255)
                                                                                                                                         gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCggg                          <  8:300685/124‑1 (MQ=255)
                                                                                                                                         gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc  >  8:262900/1‑148 (MQ=255)
                                                                                                                                         gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                               >  4:171005/1‑71 (MQ=255)
                                                                                                                                         gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGc                                                                               <  3:171005/71‑1 (MQ=255)
                                                                                                                                           tCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                  >  6:235674/1‑130 (MQ=255)
                                                                                                                                           tCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                  <  5:235674/130‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GCAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGC  >  NC_002947/1389801‑1390085

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