Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,988,780 T→C 7.8% K188E (AAA→GAA)  PP_2613 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,988,7800TC 7.8% 139.7 / 2.9 51K188E (AAA→GAA) PP_2613hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base C (2/2);  total (27/24)
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.48e-01

GGGGCGGCATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTGGCAACGCTGTGGTCAGG  >  NC_002947/2988633‑2988927
                                                                                                                                                    |                                                                                                                                                   
ggggCGGCATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTc                                                                                                                                                    >  4:3073/1‑149 (MQ=255)
  ggcggcATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCTAGTTCAAGCTGCTGGAGGGTCATAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCa                                                                                                                                                  <  5:29106/149‑1 (MQ=255)
     ggcATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTc                                                                                                                                               <  3:113579/149‑1 (MQ=255)
      gcATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCa                                                                                                                                              <  7:42716/149‑1 (MQ=255)
       cATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTc                                                                                                                                               <  4:188087/147‑1 (MQ=255)
       cATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCaa                                                                                                                                             >  3:188087/1‑148 (MQ=255)
               aGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGc                                                                                                                                     <  3:3073/149‑1 (MQ=255)
                atccTCATTTGCCAATGAGAGTACAAAGTGCTCGTTCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTTGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGTCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGc                                                                                                                                     <  4:61202/146‑1 (MQ=255)
                   cTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGcccc                                                                                                                                  >  7:26383/1‑148 (MQ=255)
                   cTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGcccc                                                                                                                                  >  3:315542/1‑148 (MQ=255)
                   cTCATGTGACAATGAGAGGACAAAGTGCTCGTGCATGAATGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGTTCAAAGGAAGTTCTGATTTTGATCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCTGc                                                                                                                                     >  3:61202/1‑145 (MQ=255)
                         tgCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGa                                                                                                                           <  7:286370/149‑1 (MQ=255)
                         tgCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGa                                                                                                                           >  8:181218/1‑149 (MQ=255)
                         tgCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGGTGa                                                                                                                           <  3:27387/149‑1 (MQ=255)
                          gCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGcc                                                                                                                                    <  7:16082/139‑1 (MQ=255)
                          gCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGcc                                                                                                                                    >  8:16082/1‑139 (MQ=255)
                           ccAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGAc                                                                                                                          <  4:315542/148‑1 (MQ=255)
                                  gagGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATg                                                                                                                  >  8:327177/1‑149 (MQ=255)
                                  gagGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATg                                                                                                                  <  4:224318/149‑1 (MQ=255)
                                          aGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCGCCCAGTTGACTCGCCATGCCAAcc                                                                                                            >  6:14261/1‑147 (MQ=255)
                                                                  gCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTg                                                                                  >  3:326312/1‑149 (MQ=255)
                                                                  gCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTg                                                                                  <  4:326312/149‑1 (MQ=255)
                                                                   ccGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTCCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGac                                                                                                       >  3:86323/1‑127 (MQ=255)
                                                                   ccGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTCCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGac                                                                                                       <  4:86323/127‑1 (MQ=255)
                                                                   ccGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTCCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGac                                                                                                       <  2:202919/127‑1 (MQ=255)
                                                                   ccGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTCCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGac                                                                                                       >  1:202919/1‑127 (MQ=255)
                                                                     gATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCTTTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCgaga                                                                                                        <  7:12074/124‑1 (MQ=255)
                                                                     gATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCgaga                                                                                                        >  8:12074/1‑124 (MQ=255)
                                                                       tCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTg                                                                                  >  6:346484/1‑144 (MQ=255)
                                                                       tCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTg                                                                                  <  5:346484/144‑1 (MQ=255)
                                                                                        gAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCAt                                                            >  2:173400/1‑149 (MQ=255)
                                                                                        gAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCAt                                                            >  3:276192/1‑149 (MQ=255)
                                                                                           ggTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCa                                                         <  7:327177/149‑1 (MQ=255)
                                                                                             tCAAAGGCAGTTCTGATTTTGAGCGGAGCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAgc                                                                >  1:180183/1‑141 (MQ=255)
                                                                                             tCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAgc                                                                <  2:180183/140‑1 (MQ=255)
                                                                                                    cAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTg                                                >  4:146243/1‑149 (MQ=255)
                                                                                                       ttCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCTTTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAgg                                                   <  1:284859/143‑1 (MQ=255)
                                                                                                       ttCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAgg                                                   >  2:284859/1‑143 (MQ=255)
                                                                                                             ttttGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGt                                       >  5:52593/1‑149 (MQ=255)
                                                                                                             ttttGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGt                                       >  3:77788/1‑149 (MQ=255)
                                                                                                                 gAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCg                                    <  6:52593/148‑1 (MQ=255)
                                                                                                                   gCGG‑GCAGCAATCGGCTTCTTTTCGTGTGTTTTCCATTCATGAGCCGGCACCCAGTTGACACGCCATGCCAACCGAGACACCGGGCCAAGCCAGCCCTTGCGATGCTCGTCGACCAGCGCATTAAGCAAGGCTCCAAGATGGTCCGgtg                                 >  2:166196/1‑149 (MQ=255)
                                                                                                                        gcagcaATCGGCTTGTTTTCGTGTGTTTTCCTTTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGa                              >  6:262946/1‑148 (MQ=255)
                                                                                                                             attCGTCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGt                        <  1:166196/147‑1 (MQ=255)
                                                                                                                             aaTCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGt                        >  2:221783/1‑149 (MQ=255)
                                                                                                                                  gCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTg                   <  8:26383/149‑1 (MQ=255)
                                                                                                                                   cTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTgg                  <  4:276192/149‑1 (MQ=255)
                                                                                                                                    ttGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTAGCCATGCCAACCGAGACACCGGGCCAAGCAAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTgg                  >  5:116051/1‑148 (MQ=255)
                                                                                                                                       ttttCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTGGCAAc              >  4:306102/1‑149 (MQ=255)
                                                                                                                                            gtgtgtTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGt                    <  7:259724/138‑1 (MQ=255)
                                                                                                                                            gtgtgtTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGt                    >  8:259724/1‑138 (MQ=255)
                                                                                                                                                 ttttCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTGGCAACGCTGTGGt      >  5:163345/1‑147 (MQ=255)
                                                                                                                                                   ttCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCa                                                                                         <  2:175057/62‑1 (MQ=255)
                                                                                                                                                   ttCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCa                                                                                         >  1:175057/1‑62 (MQ=255)
                                                                                                                                                    tCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTGGCAACGCTGTGGTCAgg  <  3:306102/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGGCGGCATGGCCAAGACCTCATGTGCCAATGAGAGGACAAAGTGCTCGTGCATGCCTGCTTGCAGCCGATCAAGTTCAAGCTGCTGGAGGGTCAAAGGCAGTTCTGATTTTGAGCGG‑GCAGCAATCGGCTTGTTTTCGTGTGTTTTCCATTCATGAGCCGGCCCCCAGTTGACTCGCCATGCCAACCGAGACACCGGGCCAAGCCAGGCCTTGCGATGGTCGTCGTCCAGCGCATTCAGCCAGGCTGCAAGATGGTCCGGTGTGACATGGTGAGTGGCAACGCTGTGGTCAGG  >  NC_002947/2988633‑2988927

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