Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,973,479 A→G 8.3% L132P (CTT→CCT)  PP_2600 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,973,4790AG 8.3% 132.1 / 2.6 48L132P (CTT→CCT) PP_2600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base G (2/2);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGCCCCC  >  NC_002947/2973335‑2973625
                                                                                                                                                    |                                                                                                                                                  
aTGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGc                                                                                                                                                >  2:126790/1‑149 (MQ=255)
 tGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCt                                                                                                                                               <  8:297160/149‑1 (MQ=255)
         aaGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTgg                                                                                                                                       <  5:136428/149‑1 (MQ=255)
         aaGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTgg                                                                                                                                       >  4:90194/1‑149 (MQ=255)
         aaGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTgg                                                                                                                                       <  4:23615/149‑1 (MQ=255)
          aGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGc                                                                                                                                      <  4:273613/149‑1 (MQ=255)
            gcggcAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGt                                                                                                                                    >  1:260729/1‑149 (MQ=255)
             cggcAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAGGCGCAGCCCTTGGCg                                                                                                                                     >  8:195193/1‑147 (MQ=255)
             cggcAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAGGCGCAGCCCTTGGCg                                                                                                                                     <  7:195193/147‑1 (MQ=255)
                   ggTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCca                                                                                                                             >  6:390305/1‑149 (MQ=255)
                       ggCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgc                                                                                                                          <  5:390305/148‑1 (MQ=255)
                        gCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgcat                                                                                                                        >  5:182551/1‑149 (MQ=255)
                        gCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCcatgcat                                                                                                                        <  2:107254/149‑1 (MQ=255)
                         cACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGCTCTAGCGCAAGCGTGGCAGCCGGCCGCAAACCTCGCGGGCATCGTAGCG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTCCGCCTCGCGTATCGCcatgcatg                                                                                                                       >  5:269752/1‑149 (MQ=255)
                            cGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGAGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTcc                                                                                                                    >  4:147072/1‑149 (MQ=255)
                             ggCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCg                                                                                                                   <  3:90194/149‑1 (MQ=255)
                              gCCATGATCGGCATGCCCATGAACGAGGGTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCgg                                                                                                                  <  6:294385/149‑1 (MQ=255)
                                              ccATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGAtgccgttg                                                                                                        >  5:369580/1‑143 (MQ=255)
                                              ccATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGAtgccgttg                                                                                                        <  6:369580/143‑1 (MQ=255)
                                                   aaCGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCaaa                                                                                             >  2:55148/1‑149 (MQ=255)
                                                      gAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTg                                                                                          >  4:227859/1‑149 (MQ=255)
                                                              gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTcc                                                                                                                    <  5:360978/115‑1 (MQ=255)
                                                              gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTcc                                                                                                                    >  6:360978/1‑115 (MQ=255)
                                                              gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCg                                                                                  >  8:115643/1‑149 (MQ=255)
                                                              gcgcATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCg                                                                                  <  2:260729/149‑1 (MQ=255)
                                                                  aTGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTa                                                                              <  1:55148/149‑1 (MQ=255)
                                                                          gcgcAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGTCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTtctc                                                                        <  5:255352/147‑1 (MQ=255)
                                                                          gcgcAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTtctc                                                                        <  7:115643/147‑1 (MQ=255)
                                                                            gcAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaa                                                                     >  4:225897/1‑148 (MQ=255)
                                                                              aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                  >  3:349355/1‑149 (MQ=255)
                                                                              aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                  <  3:227859/149‑1 (MQ=255)
                                                                              aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                  <  3:225897/149‑1 (MQ=255)
                                                                              aaGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGaaaaa                                                                  <  6:182551/149‑1 (MQ=255)
                                                                                gCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGTGCCTCCTCAGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGa                                                                      <  1:126790/147‑1 (MQ=255)
                                                                                          agGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATg                                                      >  1:330315/1‑149 (MQ=255)
                                                                                                      tCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCggg                                                  <  1:141943/141‑1 (MQ=255)
                                                                                                      tCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCggg                                                  >  2:141943/1‑141 (MQ=255)
                                                                                                            gCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCaagg                                      >  8:189564/1‑147 (MQ=255)
                                                                                                             cATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCaaggaag                                   <  3:147072/149‑1 (MQ=255)
                                                                                                              aTCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCaaggaagg                                  >  5:348749/1‑149 (MQ=255)
                                                                                                                 gTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCAt                               <  4:349355/149‑1 (MQ=255)
                                                                                                                            gcggcgATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGAcgcg                        <  2:330315/149‑1 (MQ=255)
                                                                                                                                  aTCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCAGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGccac                                                            <  6:184406/107‑1 (MQ=255)
                                                                                                                                  aTCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCAGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGccac                                                            >  5:184406/1‑107 (MQ=255)
                                                                                                                                           gatacgccAGGCGCAGCCCTTGGGGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGtt         >  5:77742/5‑149 (MQ=255)
                                                                                                                                            atacgccAGGCGCAGCCCTTGGTGTATCGCCATGCATGTACGGTATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCTAAAAACTCCACCATATGCGGGCACTCCTCAAGGAAGGCATCGACGCTCTCGGCATTGTAGttt        <  6:77742/146‑1 (MQ=255)
                                                                                                                                               cgccAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGcc     >  7:168147/1‑149 (MQ=255)
                                                                                                                                                  cAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGccccc  <  2:133190/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTTGCGCAAGCGTGGCAGCAGGGCGCCAACCTCGCGGGCATCGTAGGG‑‑‑‑CGCGGCGATCACCTGGTGGCCGCCAAGGGCTGCGCCTGGCGTATCGCCATGCATGTCCGGGATGCCGTTGCCGTCGGCAAACTGCAACGCCGTGTACTTCTCGAAAAACGCCACCATATGCGGGCACGCCTCAAGGAAGGCATCGACGCGCTCGGCATTGTAGTTTGCCCCC  >  NC_002947/2973335‑2973625

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