Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,291,043 A→G 8.5% Q95R (CAG→CGG)  PP_2893 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,291,0430AG 8.5% 132.4 / 2.2 48Q95R (CAG→CGG) PP_2893hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/18);  new base G (2/2);  total (27/21)
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.84e-01

GACTTCCAGGGGCGCCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCGATCTGCTTG  >  NC_002947/3290898‑3291192
                                                                                                                                                 |                                                                                                                                                     
gACTTCCAGGGGCGCCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTgcaggca                                                                                                                                                   <  1:375245/150‑1 (MQ=255)
   ttCCAGGGGCGCCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGGAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTacaggcaggc                                                                                                                                                >  2:241524/1‑148 (MQ=255)
     ccAGGGGCGCCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCg                                                                                                                                              >  2:329647/1‑150 (MQ=255)
            cgcCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATAGACCAGTTGCAGGCAGACCGCATTGAc                                                                                                                                       <  1:22896/150‑1 (MQ=255)
             gcCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACa                                                                                                                                      >  1:393121/1‑150 (MQ=255)
              ccTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACAt                                                                                                                                     <  1:645061/150‑1 (MQ=255)
               ccggCGGGAAGTTGGCAGGCCAGCCGGCTAACCTCGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCCGGCAGCCCGCCTTGACATg                                                                                                                                    <  1:124359/148‑1 (MQ=255)
                       cAGTGGGCAGGGCAGCCGGCTAACCTGGGCCGACGGGGGCAGCACGGTAGAGATGCATGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGCCCAGTTGCAGGCAGACCGCATTGACATGGCCTTTgc                                                                                                                            <  2:142570/150‑1 (MQ=255)
                                 gCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCAcc                                                                                                                  <  1:329647/150‑1 (MQ=255)
                                  ccAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCATTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCg                                                                                                                 >  2:321289/1‑150 (MQ=255)
                                           tAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCg                                                                                                        >  1:153165/1‑150 (MQ=255)
                                              ccTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCgaga                                                                                                     <  2:393121/150‑1 (MQ=255)
                                              ccTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCgaga                                                                                                     <  1:820762/150‑1 (MQ=255)
                                               cTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGagag                                                                                                    >  2:4070/1‑150 (MQ=255)
                                                    ccGAAGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGccc                                                                                               >  2:678077/1‑150 (MQ=255)
                                                              cagcaCGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAAtt                                                                                     >  1:147354/1‑150 (MQ=255)
                                                                                ggCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGcc                                                                   >  1:421060/1‑150 (MQ=255)
                                                                                 gcgcCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCGGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCg                                                                                             >  1:826200/1‑123 (MQ=255)
                                                                                 gcgcCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCGGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCg                                                                                             <  2:826200/123‑1 (MQ=255)
                                                                                   gcCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCt                                                                >  1:714774/1‑150 (MQ=255)
                                                                                    ccGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTg                                                               <  1:638329/150‑1 (MQ=255)
                                                                                     cGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTgg                                                              >  1:440320/1‑150 (MQ=255)
                                                                                           tGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCg                                                        >  2:4616/1‑150 (MQ=255)
                                                                                           tGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCg                                                        >  1:99635/1‑150 (MQ=255)
                                                                                            ggCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGa                                                       >  2:640548/1‑150 (MQ=255)
                                                                                             gCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGaa                                                      <  2:592011/150‑1 (MQ=255)
                                                                                                tgttCGCCTGCCTGCTTGCTGCGCGCGACCTTATGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTg                                                   >  2:533733/4‑150 (MQ=255)
                                                                                                tGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCa                                                           <  2:75907/142‑1 (MQ=255)
                                                                                                tGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCa                                                           >  1:75907/1‑142 (MQ=255)
                                                                                                tGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTg                                                   >  1:174325/1‑150 (MQ=255)
                                                                                                tGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTg                                                   <  2:153165/150‑1 (MQ=255)
                                                                                                        tgcctgcTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGt                                           >  1:171005/1‑150 (MQ=255)
                                                                                                          cctgcTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCa                                         >  1:314752/1‑150 (MQ=255)
                                                                                                           ctgcTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACGTGGCCTTTGCGCCAAGCACCGATCCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCaa                                        <  2:97297/150‑1 (MQ=255)
                                                                                                              cttgctGCGCGCGACCTTGTGCATTGACCAGTTGCGGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGACCATCGa                                                       >  1:627875/1‑132 (MQ=255)
                                                                                                              cttgctGCGCGCGACCTTGTGCATTGACCAGTTGCGGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGACCATCGa                                                       <  2:627875/132‑1 (MQ=255)
                                                                                                                  ctgcGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCCCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCggg                                 <  1:813563/150‑1 (MQ=255)
                                                                                                                                gtgCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGa                   <  1:121039/150‑1 (MQ=255)
                                                                                                                                  gCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATg                 <  2:171005/150‑1 (MQ=255)
                                                                                                                                   cATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGAt                  <  2:6002/148‑1 (MQ=255)
                                                                                                                                   cATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGAt                  >  1:6002/1‑148 (MQ=255)
                                                                                                                                   cATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATgg                <  2:174325/150‑1 (MQ=255)
                                                                                                                                     ttGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCa              >  2:345033/1‑150 (MQ=255)
                                                                                                                                      tGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGTAg             >  2:220022/1‑150 (MQ=255)
                                                                                                                                       gACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGTCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGc            <  1:4070/150‑1 (MQ=255)
                                                                                                                                       gACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGc            >  1:796385/1‑150 (MQ=255)
                                                                                                                                        aCCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCg           >  1:293350/1‑150 (MQ=255)
                                                                                                                                         ccAGTTGCAGGCAGACAGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGc                                                                                   >  2:121446/1‑77 (MQ=255)
                                                                                                                                         ccAGTTGCAGGCAGACAGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGc                                                                                   <  1:121446/77‑1 (MQ=255)
                                                                                                                                          cAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCGAt         >  1:397996/1‑150 (MQ=255)
                                                                                                                                                 aggcagACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGGCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCGATCTGCTTg  <  2:421060/150‑1 (MQ=255)
                                                                                                                                                 aggcagACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCGATCTGCTTg  >  1:660210/1‑150 (MQ=255)
                                                                                                                                                 aggcagACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCGATCTGCTTg  >  1:610806/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACTTCCAGGGGCGCCTGGCGGGCAGTTGGCAGGCCAGCCGGCTAACCTGGGCCGACGGCGGCAGCACGGTAGAGATGCAGGCGCCGTTGCTGGCCTGGTCGCCTGCCTGCTTGCTGCGCGCGACCTTGTGCATTGACCAGTTGCAGGCAGACCGCATTGACATGGCCTTTGCGCCAAGCACCGAACCCACCGAGAGCGCCCCGCTGCAATTGCCGACTCTGCGCTTGCCGCTGGCCATCGAACTGGGCGAGGTCAAGGTCGGGCAGCTGCGCCTGGATGGCAGCGATCTGCTTG  >  NC_002947/3290898‑3291192

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