Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,241,031 T→G 7.4% L182R (CTG→CGG)  PP_4618 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,241,0310TG 7.4% 152.2 / 3.7 54L182R (CTG→CGG) PP_4618hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base G (2/2);  total (27/27)
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.74e-01

CCGACCACCTGCTCGATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCGCGCCTCTCTCTACCATCATCTGG  >  NC_002947/5240887‑5241168
                                                                                                                                                |                                                                                                                                         
ccGACCACCTGCTCGATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCAc                                                                                                                                      <  2:891757/150‑1 (MQ=255)
 cGACCACCTGCTCGATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACg                                                                                                                                     >  1:1140281/1‑150 (MQ=255)
  gACCACCTGCTCGATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGc                                                                                                                                    >  2:805668/1‑150 (MQ=255)
       ccTGCTCGATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCTCGCTGCAGATGCTGCAACGCCTGTCACGCAACCt                                                                                                                               >  2:1048338/1‑150 (MQ=255)
              gATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCaa                                                                                                                        >  1:50716/1‑150 (MQ=255)
                tCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCATCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAAc                                                                                                                                 >  1:330177/1‑139 (MQ=255)
                tCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCATCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAAc                                                                                                                                 <  2:330177/139‑1 (MQ=255)
                  ccTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGc                                                                                                                    <  1:1048338/150‑1 (MQ=255)
                   cTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGca                                                                                                                   <  2:190851/150‑1 (MQ=255)
                   cTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGca                                                                                                                   >  2:1146830/1‑150 (MQ=255)
                   cTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGca                                                                                                                   <  1:1181336/150‑1 (MQ=255)
                             cggcGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAAc                                                                                                         <  1:474343/150‑1 (MQ=255)
                                cgcCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGAc                                                                                                      <  2:942232/150‑1 (MQ=255)
                                 gcCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAACGCCAGACCGCCGCAGATGCATCGACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACa                                                                                                     >  1:445035/1‑150 (MQ=255)
                                        tCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGTCGCTGACAGCCGAGCACCTGTGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAaccac                                                                                                >  1:1052067/1‑148 (MQ=255)
                                        tCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGTCGCTGACAGCCGAGCACCTGTGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAaccac                                                                                                <  2:1052067/148‑1 (MQ=255)
                                           cAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAAc                                                                                                                                 <  1:136829/112‑1 (MQ=255)
                                           cAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAAc                                                                                                                                 >  2:136829/1‑112 (MQ=255)
                                           cAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTCACGACAACCACCTggc                                                                                           >  1:864505/1‑150 (MQ=255)
                                            aGCGACTGCCGCTATTTGCTTGCTCTCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTggcg                                                                                          <  2:374498/150‑1 (MQ=255)
                                            aGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTggcg                                                                                          <  2:1140281/150‑1 (MQ=255)
                                            aGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTg                                                                                             >  2:283513/1‑147 (MQ=255)
                                            aGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTg                                                                                             <  1:283513/147‑1 (MQ=255)
                                                aCTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCAt                                                                                      <  2:1156293/150‑1 (MQ=255)
                                                    ccGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAAc                                                                                                                                 <  1:1042953/103‑1 (MQ=255)
                                                    ccGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAAc                                                                                                                                 >  2:1042953/1‑103 (MQ=255)
                                                    ccGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGAt                                                                                  <  2:327141/150‑1 (MQ=255)
                                                     cgcgAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTaa                                                                                                          >  1:329811/1‑125 (MQ=255)
                                                     cgcgAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTaa                                                                                                          <  2:329811/125‑1 (MQ=255)
                                                       cgAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGa                                                                               <  1:1146830/150‑1 (MQ=255)
                                                               tGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGt                                                                                                            <  1:1163866/113‑1 (MQ=255)
                                                               tGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGt                                                                                                            >  2:1163866/1‑113 (MQ=255)
                                                                ggCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCg                                                                      <  2:242684/150‑1 (MQ=255)
                                                                       cAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATg                                                               >  2:686944/1‑150 (MQ=255)
                                                                            cGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCATCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGtt                                                          >  1:1001651/1‑150 (MQ=255)
                                                                             ggCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCATCGCCTGTCACGCAACCTGCGCCAACTGCACGCCATTAACGACATCCTCCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTg                                                         <  1:750760/150‑1 (MQ=255)
                                                                             ggCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTg                                                         >  1:128854/1‑150 (MQ=255)
                                                                                ccGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGc                                                      >  2:571355/1‑150 (MQ=255)
                                                                                    aCCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCagag                                                                        <  2:970247/128‑1 (MQ=255)
                                                                                    aCCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCagag                                                                        >  1:970247/1‑128 (MQ=255)
                                                                                           gcTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGt                                           >  2:853608/1‑150 (MQ=255)
                                                                                                 aGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCaga                                                                         >  2:573688/1‑114 (MQ=255)
                                                                                                 aGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCaga                                                                         <  1:573688/114‑1 (MQ=255)
                                                                                                    cGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTa                                  <  1:571355/150‑1 (MQ=255)
                                                                                                           cTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACt                           >  2:930949/1‑150 (MQ=255)
                                                                                                                        gCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCGCGCctctctct              <  2:445035/150‑1 (MQ=255)
                                                                                                                                   aTGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCGCGCCTCTCTCTACCATCATCTg   >  2:1053841/1‑150 (MQ=255)
                                                                                                                                    tgctgcAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCGCGCCTCTCGCTACCATCATCTgg  >  1:248659/1‑150 (MQ=255)
                                                                                                                                        gcAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCgcgc                      >  2:484746/1‑126 (MQ=255)
                                                                                                                                        gcAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCgcgc                      <  1:484746/126‑1 (MQ=255)
                                                                                                                                            cGCCGGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGtt                                                          >  1:244532/1‑86 (MQ=255)
                                                                                                                                            cGCCGGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGtt                                                          <  2:244532/86‑1 (MQ=255)
                                                                                                                                            cGCCGGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGtt                                                          >  1:244261/1‑86 (MQ=255)
                                                                                                                                            cGCCGGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGtt                                                          <  2:244261/86‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                         
CCGACCACCTGCTCGATCCCTGCGGCGGCCGGCGCCTGTATCCCAGCGACTGCCGCGAATTGCTGGCGCGCCAGTTCGGCCCGCACCTGGCGCTGACAGCCGAGCACCTGCGCAGCGCCAGCCCGCTGCAGATGCTGCAACGCCTGTCACGCAACCTGCGCCAACTGCACATCAGTAACGACAACCACCTGGCGGCATTGATCGATTCAGAGCGGATCATGCAGTTGGGCCCGGTGCAGGTGAGCGACTACATGACTCGCGCCTCTCTCTACCATCATCTGG  >  NC_002947/5240887‑5241168

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