Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,422,785 C→G 11.3% T453T (ACC→ACG)  PP_1246 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,422,7850CG11.3% 153.4 / 10.9 53T453T (ACC→ACG) PP_1246hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base G (2/4);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

ATCACCCGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGACCTCGGCCAGCGTAGT  >  NC_002947/1422645‑1422934
                                                                                                                                            |                                                                                                                                                     
aTCACCCGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACGGGTGTGcac                                                                                                                                              <  1:599616/150‑3 (MQ=255)
  cACCCGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGa                                                                                                                                            <  1:791165/150‑1 (MQ=255)
   aCCCGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAg                                                                                                                                           >  1:451058/1‑150 (MQ=255)
    cccGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTTGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGc                                                                                                                                          <  1:434243/150‑1 (MQ=255)
     ccGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCg                                                                                                                                         >  2:17885/1‑150 (MQ=255)
            cGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCTCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACcag                                                                                                                                  <  2:606232/150‑1 (MQ=255)
            cGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACcag                                                                                                                                  <  2:56425/150‑1 (MQ=255)
             gCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACcagc                                                                                                                                 >  1:945755/1‑150 (MQ=255)
                 cGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGtt                                                                                                                             <  2:839569/150‑1 (MQ=255)
                   agcagaaCCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGc                                                                                                                           >  1:38462/6‑150 (MQ=255)
                     aTCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAg                                                                                                                         >  2:471005/1‑150 (MQ=255)
                          ccGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACgg                                                                                                                    >  2:956270/1‑150 (MQ=255)
                          ccGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACgg                                                                                                                    >  2:183002/1‑150 (MQ=255)
                          ccGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACgg                                                                                                                    <  2:470964/150‑1 (MQ=255)
                             ccAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACGGGTGTGcacg                                                                                                                                             >  1:970122/1‑119 (MQ=255)
                             ccAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACGGGTGTGcacg                                                                                                                                             <  2:970122/122‑4 (MQ=255)
                                aGGCCTGTGCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTAGCCCGTGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAg                                                                                                              <  2:1118228/150‑1 (MQ=255)
                                    cTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGAc                                                                                                          >  2:140997/1‑150 (MQ=255)
                                      gTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACtt                                                                                                        >  2:783356/1‑150 (MQ=255)
                                       ttcttcGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACGTc                                                                                                       >  2:113146/1‑150 (MQ=255)
                                                gtgCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGgcgc                                                                                              <  2:1109039/150‑1 (MQ=255)
                                                     cgACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAAc                                                                                         <  1:140997/150‑1 (MQ=255)
                                                        cGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGa                                                                                      <  2:763646/150‑1 (MQ=255)
                                                         ggAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAg                                                                                     <  1:956270/150‑1 (MQ=255)
                                                          gAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAgg                                                                                    >  2:775603/1‑150 (MQ=255)
                                                               aaCTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTg                                                                               <  2:33271/150‑1 (MQ=255)
                                                               aaCTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTg                                                                               >  1:108203/1‑150 (MQ=255)
                                                               aaCTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTg                                                                               <  2:108203/150‑1 (MQ=255)
                                                               aaCTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTg                                                                               <  1:783356/150‑1 (MQ=255)
                                                                  tGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACGGGTGTGcacgc                                                                                                                                            <  2:723625/86‑5 (MQ=255)
                                                                  tGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACGGGTGTGcacgc                                                                                                                                            >  1:723625/1‑82 (MQ=255)
                                                                                 tcacagcagcGTGTGCATGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATc                                                             >  1:38557/4‑150 (MQ=255)
                                                                                    cagcagcGTGTGCCTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCATACCCGTGTGCTGGCGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAg                                                          <  2:603472/150‑1 (MQ=255)
                                                                                     agcagcGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCAt                                                                                                                 >  1:1101201/1‑94 (MQ=255)
                                                                                     agcagcGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCAt                                                                                                                 <  2:1101201/94‑1 (MQ=255)
                                                                                                 gTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGcgcttcgctt                                             <  1:85983/150‑1 (MQ=255)
                                                                                                  tGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCgcttcgcttc                                            <  1:977460/150‑1 (MQ=255)
                                                                                                   gAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGa                                                                                      <  2:2443/107‑1 (MQ=255)
                                                                                                   gAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGa                                                                                      >  1:2443/1‑107 (MQ=255)
                                                                                                       agcGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCAGGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCAc                                       <  1:688007/150‑1 (MQ=255)
                                                                                                        gcGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGAGCGCGCAGATGGCCGGTATCAAGGTGCGCTTAGCTTCGCTACt                                      >  1:63021/1‑150 (MQ=255)
                                                                                                              gcgcCAACAGGGGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGagca                                <  2:66979/150‑1 (MQ=255)
                                                                                                                  cAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCaa                            >  2:348279/1‑150 (MQ=255)
                                                                                                                  cAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGATTCGCTTCGCCACTGGAGCAGCaa                            >  2:864449/1‑150 (MQ=255)
                                                                                                                    aCAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCaaaa                          >  2:600016/1‑150 (MQ=255)
                                                                                                                      aGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGCTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAAcc                        >  1:139515/1‑150 (MQ=255)
                                                                                                                          gCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGc                    <  2:620868/150‑1 (MQ=255)
                                                                                                                             ggCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATCAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGAc                 >  2:98386/1‑150 (MQ=255)
                                                                                                                               cccGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCt                      >  2:94113/1‑143 (MQ=255)
                                                                                                                               cccGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCt                      <  1:94113/143‑1 (MQ=255)
                                                                                                                                tctccagCACACGGGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGACCTc              <  1:241968/144‑1 (MQ=255)
                                                                                                                                ccGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGACCTc              <  1:1105886/150‑1 (MQ=255)
                                                                                                                                          acCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGACCTCGGCCAGCGTa    >  2:676497/1‑150 (MQ=255)
                                                                                                                                            ccGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGCGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGACCTCGGCCAGCGTAGt  >  2:749493/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
ATCACCCGCAACCGCATCGGCATCAACCGCCAAGGCCTGTTCTTCGATGTGCGCGACGGAGCCAACTGGGCGGTACAGCGCGAGCAGCAGCGTGTGCGTGAGCAGCGTCAGCGCCAACAGGAGCTGGCCCGCGTGCACACCCGTGTGCTGGAGCGTTACCAGCAGTTGCAGGACGGCATGAGCGACTTCAAGGGGCGCGAAACCGAGGCGTTGGCGCAGATGGCCGGTATCAAGGTGCGCTTCGCTTCGCCACTGGAGCAGCAAAACCCTGCGACCTCGGCCAGCGTAGT  >  NC_002947/1422645‑1422934

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