Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,375,486 C→A 16.3% A32S (GCG→TCG)  PP_5637 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,375,4860CA16.3% 119.9 / 12.8 49A32S (GCG→TCG) PP_5637hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base A (5/3);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

ACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCACCAATACTGTATA  >  NC_002947/4375339‑4375634
                                                                                                                                                   |                                                                                                                                                     
accacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGa                                                                                                                                                     >  1:620157/1‑150 (MQ=255)
 ccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTCGGGTTAAACACCCAGGCCTGGGTCGCGAc                                                                                                                                                    >  1:557963/1‑150 (MQ=255)
             atatAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGGTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAgg                                                                                                                                        <  2:858771/150‑1 (MQ=255)
               ataGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGTCTGGGTCGAGTCGCGACCCaggcctg                                                                                                                                      <  2:44452/150‑6 (MQ=255)
                  gCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGc                                                                                                                                   <  2:1061463/150‑1 (MQ=255)
                         ttagACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATg                                                                                                                            <  1:908656/147‑1 (MQ=255)
                                      aCCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTgg                                                                                                               >  1:551280/1‑150 (MQ=255)
                                        cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGAc                                                                                                             <  1:546451/150‑1 (MQ=255)
                                              gatgttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgc                                                                                                       <  1:823455/150‑1 (MQ=255)
                                                 gttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgcg                                                                                                    >  1:808730/1‑150 (MQ=255)
                                                          gAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTg                                                                                           >  1:1028656/1‑150 (MQ=255)
                                                                  ccATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTcc                                                                                   >  2:631026/1‑150 (MQ=255)
                                                                     tGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTa                                                                                <  2:828854/150‑1 (MQ=255)
                                                                       ccGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAgg                                                                              <  2:551280/150‑1 (MQ=255)
                                                                        cGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGCGGGGCTGGACTTTAAGTAggg                                                                             >  2:780603/1‑149 (MQ=255)
                                                                            gACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTc                                                                         >  1:1185310/1‑150 (MQ=255)
                                                                              ccTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTg                                                                       >  1:9881/1‑150 (MQ=255)
                                                                               cTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGa                                                                      >  2:631178/1‑150 (MQ=255)
                                                                                     aTCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCaa                                                                >  2:113201/1‑150 (MQ=255)
                                                                                      tCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAg                                                               <  2:808730/150‑1 (MQ=255)
                                                                                        gTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGc                                                             <  1:629153/150‑1 (MQ=255)
                                                                                         tGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGcc                                                            >  1:96037/1‑150 (MQ=255)
                                                                                           cAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCAc                                                          <  2:620157/150‑1 (MQ=255)
                                                                                             gCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGc                                                        >  1:70017/1‑150 (MQ=255)
                                                                                                gggggATAATGTCGGGCCCAGGCGTAGGGTTAAACAGCCAGGCCTGGGTCGCGGCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCAt                                                     >  2:718725/1‑150 (MQ=255)
                                                                                                 ggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGa                                                                                         >  2:318705/1‑113 (MQ=255)
                                                                                                 ggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGa                                                                                         <  1:318705/113‑1 (MQ=255)
                                                                                                    gATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCATGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTc                                                 <  1:780603/150‑1 (MQ=255)
                                                                                                           tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                   <  2:1267059/77‑1 (MQ=255)
                                                                                                           tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                   >  1:1267059/1‑77 (MQ=255)
                                                                                                            cGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTc                                         <  1:1183779/150‑1 (MQ=255)
                                                                                                             gggCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTACGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCg                                        >  2:1216829/1‑150 (MQ=255)
                                                                                                               gCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGaa                                      >  1:626142/1‑150 (MQ=255)
                                                                                                                cccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAAc                                     <  2:754169/150‑1 (MQ=255)
                                                                                                                   aGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                                  <  2:70017/150‑1 (MQ=255)
                                                                                                                            gTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTa                         <  1:113201/150‑1 (MQ=255)
                                                                                                                                 acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                    >  2:456934/1‑150 (MQ=255)
                                                                                                                                 acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                    <  2:1028656/150‑1 (MQ=255)
                                                                                                                                 acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                    <  1:631026/150‑1 (MQ=255)
                                                                                                                                   acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑gg                 >  2:219108/1‑150 (MQ=255)
                                                                                                                                     ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCa               >  2:1246590/1‑150 (MQ=255)
                                                                                                                                      agcctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGcc                                                       >  2:155157/6‑110 (MQ=255)
                                                                                                                                      agcctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGcc                                                       <  1:155157/105‑1 (MQ=255)
                                                                                                                                      agcctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCAc              >  1:210887/6‑150 (MQ=255)
                                                                                                                                         ctcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCa                                                           >  1:97934/3‑103 (MQ=255)
                                                                                                                                         ctcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCa                                                           <  2:97934/101‑1 (MQ=255)
                                                                                                                                         ctcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCACCaa           >  2:1285606/3‑150 (MQ=255)
                                                                                                                                         ctcTGGGTCGAGTCGCGACCCAGAGGCTAGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAGGGAAGGAC‑TGCGGGCACCaa           >  2:614568/3‑150 (MQ=255)
                                                                                                                                               gTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCACCAATACTGt     >  2:1033303/1‑150 (MQ=255)
                                                                                                                                                 cgcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCTGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCACCAATACTGtat   >  2:149110/1‑150 (MQ=255)
                                                                                                                                                  gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAg                                                               <  1:1011359/90‑1 (MQ=255)
                                                                                                                                                  gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAg                                                               >  2:1011359/1‑90 (MQ=255)
                                                                                                                                                  gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCACCAATACTGtata  <  1:117251/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGC‑GGCACCAATACTGTATA  >  NC_002947/4375339‑4375634

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