Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,184,158 A→C 9.8% M83L (ATG→CTG)  PP_4564 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,184,1580AC 9.8% 102.2 / 4.8 41M83L (ATG→CTG) PP_4564lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (19/18);  new base C (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGT  >  NC_002947/5184014‑5184307
                                                                                                                                                |                                                                                                                                                     
cTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTc                                                                                                                                                  >  1:201020/1‑150 (MQ=255)
  ggTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGc                                                                                                                                                <  1:506355/150‑1 (MQ=255)
      ccGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcgac                                                                                                                                            >  2:46091/1‑150 (MQ=255)
              gatGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAggg                                                                                                                                    >  2:375996/3‑149 (MQ=255)
                           cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                       >  2:436777/1‑150 (MQ=255)
                           cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                       >  1:411366/1‑150 (MQ=255)
                              gACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTg                                                                                                                    >  2:531271/1‑150 (MQ=255)
                                cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgca                                                                                                                  <  1:64057/150‑1 (MQ=255)
                                       cgcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGAGATGGTCGCCGCGGCCCAGGTCACCGACCAGGCCCTGCAGCAAAGc                                                                                                           <  1:398598/150‑1 (MQ=255)
                                                aGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAg                                                                                                  >  1:221556/1‑150 (MQ=255)
                                                 gggCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAgg                                                                                                 >  1:334614/1‑150 (MQ=255)
                                                         ctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAg                                                                                         >  1:94941/1‑150 (MQ=255)
                                                           ggATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTGCCTCGACTGGCTCGACCGGCTACGTGGCATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGATGCCCGTgttgc                                                                                       <  2:562884/150‑5 (MQ=255)
                                                                     ttGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGa                                                                             >  2:192653/1‑150 (MQ=255)
                                                                       gcgcGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                           <  2:94941/150‑1 (MQ=255)
                                                                             gCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                     <  1:476589/150‑1 (MQ=255)
                                                                             gCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                     >  1:564934/1‑150 (MQ=255)
                                                                             gCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                     <  2:197934/150‑1 (MQ=255)
                                                                                      ctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                            >  1:173818/1‑150 (MQ=255)
                                                                                            ttGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTCa                                                                                                                                  <  1:243393/74‑1 (MQ=25)
                                                                                            ttGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTCa                                                                                                                                  >  2:243393/1‑74 (MQ=25)
                                                                                               cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                             <  1:119690/92‑1 (MQ=255)
                                                                                               cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                             >  2:119690/1‑92 (MQ=255)
                                                                                                  gACACAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAAATGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAc                                                <  1:378239/150‑1 (MQ=255)
                                                                                                   aCTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAcc                                               >  2:310310/1‑150 (MQ=255)
                                                                                                     tCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCat                                                                                                                                      <  1:50362/61‑2 (MQ=37)
                                                                                                     tCAACTGCCCGGCTACCTCGACTGGCGCGACCGGCTACGTGACCTGGTCGTCGGCGACCat                                                                                                                                      >  2:50362/1‑60 (MQ=25)
                                                                                                       aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                      <  2:459177/59‑1 (MQ=255)
                                                                                                       aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                      >  1:459177/1‑59 (MQ=255)
                                                                                                       aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                           >  2:49668/1‑150 (MQ=255)
                                                                                                               ggCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTgcgc                                   <  1:531271/150‑1 (MQ=255)
                                                                                                                gCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGcgcg                                  <  1:436777/150‑1 (MQ=255)
                                                                                                                gCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGAGCATGAGGCAATTGGTCGATGGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGcgcg                                  >  1:293457/1‑150 (MQ=255)
                                                                                                                    ccTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGAGCGGCAt                              <  2:564934/150‑1 (MQ=255)
                                                                                                                      tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGa                            <  2:173818/150‑1 (MQ=255)
                                                                                                                               tCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCGCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAcc                   <  2:89570/150‑1 (MQ=255)
                                                                                                                                cGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCa                  <  1:310310/150‑1 (MQ=255)
                                                                                                                                 gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAg                 >  1:248229/1‑150 (MQ=255)
                                                                                                                                             gACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATg     <  2:334614/150‑1 (MQ=255)
                                                                                                                                             gACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATg     <  1:46091/150‑1 (MQ=255)
                                                                                                                                               cATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCg   >  1:243470/1‑150 (MQ=255)
                                                                                                                                                aTGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGt  <  2:411366/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTGGTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGT  >  NC_002947/5184014‑5184307

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