Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,097,447 T→C 7.8% S14G (AGC→GGC)  PP_0954 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,097,4470TC 7.8% 139.5 / 2.1 51S14G (AGC→GGC) PP_0954hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base C (2/2);  total (25/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 = 9.20e-01

ATCATAGAACCCTGAGTGATGATCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGATTCTAACGGC  >  NC_002947/1097299‑1097589
                                                                                                                                                    |                                                                                                                                               
atcatAGAACCCTGAGTGATGATCACGTCGCCGGAGTAAGTGGCCACGCATTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCt                                                                                                                                                 >  5:162354/1‑149 (MQ=255)
   atAGAACCCTGAGTGATGATCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAg                                                                                                                                              <  2:230192/149‑1 (MQ=255)
       aaCCCTGAGTGATGATCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTACGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCt                                                                                                                                           >  7:258977/1‑148 (MQ=255)
             gAGTGATGATCACTGCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGa                                                                                                                                         <  6:247850/144‑1 (MQ=255)
             gAGTGATGATCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGa                                                                                                                                         >  5:247850/1‑144 (MQ=255)
                  atgatCACGTCGCCGGTGTAAGTGGCCACTCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAgggg                                                                                                                               <  6:162354/149‑1 (MQ=255)
                  atgatCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAggg                                                                                                                                <  2:267489/148‑1 (MQ=255)
                  atgatCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAggg                                                                                                                                <  3:169653/148‑1 (MQ=255)
                     atCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAg                                                                                                                            >  5:96616/1‑149 (MQ=255)
                          gTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGttt                                                                                                                       >  4:59939/1‑149 (MQ=255)
                              ccGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGCGGGCGCTCGCGCTTCCCAGTGCTGCGCCCAGGCTGAGCAAAAGGGGGAGGGTTTTAAc                                                                                                                   >  6:266661/1‑149 (MQ=255)
                              ccGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCCCAGGCTGAGCAAAAGGGGGAGGGTTTTAAc                                                                                                                   >  8:327470/1‑149 (MQ=255)
                                ggTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGa                                                                                                                 <  6:96616/149‑1 (MQ=255)
                                  tgtAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGc                                                                                                               <  2:66746/149‑1 (MQ=255)
                                               gCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGGTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGtttt                                                                                                                      <  8:5478/129‑1 (MQ=255)
                                               gCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGtttt                                                                                                                      >  7:5478/1‑129 (MQ=255)
                                                           aTCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCCCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTag‑a                                                                                      <  7:327470/148‑1 (MQ=255)
                                                           aTCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCCCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTag‑a                                                                                      <  5:266661/148‑1 (MQ=255)
                                                             ccAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTag‑agag                                                                                   >  1:221651/1‑149 (MQ=255)
                                                               aGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGc                                                                                  >  7:243756/1‑148 (MQ=255)
                                                                  tGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGt                                                                              <  8:243756/149‑1 (MQ=255)
                                                                           tCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGc                                                                     >  6:229909/1‑149 (MQ=255)
                                                                            cGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGcc                                                                    >  7:18332/1‑149 (MQ=255)
                                                                              gCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATa                                                                                                         <  4:81012/111‑1 (MQ=255)
                                                                                       cGGATAGGCTGGTCACGGTCGTGCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCaaa                                                          <  8:258977/148‑1 (MQ=255)
                                                                                        ggATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATa                                                        >  2:97337/1‑149 (MQ=255)
                                                                                             ggCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCtt                                                   >  7:183128/1‑149 (MQ=255)
                                                                                                tGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATAAGCTTTCa                                                >  3:43023/1‑149 (MQ=255)
                                                                                                     aCGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCt                                           >  8:44534/1‑149 (MQ=255)
                                                                                                      cGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTAGCCTTCTTTCAAATACGCTTTCATGCCtt                                          >  2:270529/1‑149 (MQ=255)
                                                                                                       ggTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCAATCTTGCCTTCTTTCAAATACGCTTTCATGCCttt                                         >  3:305688/1‑149 (MQ=255)
                                                                                                         tCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGc                                       <  5:103933/149‑1 (MQ=255)
                                                                                                         tCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGc                                       >  6:103933/1‑149 (MQ=255)
                                                                                                         tCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGc                                       <  2:221651/149‑1 (MQ=255)
                                                                                                           gTTCGGCAGGGCGAGGGCGCTCTCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGa                                                                                                                             <  8:170255/62‑1 (MQ=255)
                                                                                                           gTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGccc                                     <  4:33040/149‑1 (MQ=255)
                                                                                                               ggCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTa                                  <  4:43023/148‑1 (MQ=255)
                                                                                                                  agggcgagggcgCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAGCAGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAgtg                               <  1:270529/149‑1 (MQ=255)
                                                                                                                   gggcgagggcgCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAgtggt                             <  8:18332/149‑1 (MQ=255)
                                                                                                                       gagggcgCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTg                                        <  4:272976/134‑1 (MQ=255)
                                                                                                                       gagggcgCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTg                                        >  3:272976/1‑134 (MQ=255)
                                                                                                                          ggcgCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCg                      >  6:293868/1‑149 (MQ=255)
                                                                                                                          ggcgCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCg                      >  8:167886/1‑149 (MQ=255)
                                                                                                                            cgcTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGgc                    <  5:229909/149‑1 (MQ=255)
                                                                                                                                    cTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGa            <  1:97337/149‑1 (MQ=255)
                                                                                                                                       cccAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGAtt          <  7:44534/148‑1 (MQ=255)
                                                                                                                                       cccAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGATTc         >  8:59143/1‑149 (MQ=255)
                                                                                                                                              ctgcGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGcc                 <  4:273607/134‑1 (MQ=255)
                                                                                                                                              ctgcGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGcc                 >  3:273607/1‑134 (MQ=255)
                                                                                                                                              ctgcGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGATTCTAACGGc  >  6:193302/1‑149 (MQ=255)
                                                                                                                                                 cgcTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGATTCTAAc     >  6:169719/1‑143 (MQ=255)
                                                                                                                                                 cgcTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGATTCTAAc     <  5:169719/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATCATAGAACCCTGAGTGATGATCACGTCGCCGGTGTAAGTGGCCACGCCTTGCTTGTCATCCAGGTGAGCGTTGTCGGCCTGGATGCGGATAGGCTGGTCACGGTCGTTCGGCAGGGCGAGGGCGCTCGCGCTTCCCAGTGCTGCGCTCAGGCTGAGCAAAAGGGGGAGGGTTTTAACGAGCCTCATACTGTCCTCTTACGTTAG‑AGAGCAAGTCCATCTTGCCTTCTTTCAAATACGCTTTCATGCCTTTGCCCGTAGTGGTGCCTCCGGCGCCGTCGATTCTAACGGC  >  NC_002947/1097299‑1097589

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