Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,652,524 G→A 8.3% G1415G (GGC→GGT)  PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,652,5240GA 8.3% 136.3 / 2.7 48G1415G (GGC→GGT) PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/27);  new base A (2/2);  total (19/29)
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.10e-01

GCCCAGATCCACCGGTGAGTTGGACGCATCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCGGCCTTTTTCTC  >  NC_002947/1652377‑1652670
                                                                                                                                                   |                                                                                                                                                  
gCCCAGATCCCCCGGTGAGTTGGACGCCTCCAACGACAGGTCTCGGCGATAGGCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGc                                                                                                                                                   <  4:236919/149‑1 (MQ=255)
       tCCACCGGTGAGTTGGACGCATCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGAc                                                                                                                                            <  3:16699/149‑1 (MQ=255)
                tcgTTGGATGCGTCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTc                                                                                                                                   >  1:90893/3‑149 (MQ=255)
                   ttGGACGCATCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCggt                                                                                                                                <  5:136484/149‑3 (MQ=255)
                   ttGGACGCATCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGAt                                                                                                                                >  7:197599/1‑149 (MQ=255)
                   ttGGACGCATCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGAt                                                                                                                                <  6:130455/149‑1 (MQ=255)
                           aTCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCgc                                                                                                                        <  1:329452/149‑1 (MQ=255)
                                aCGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCtt                                                                                                                   <  1:231111/149‑1 (MQ=255)
                                aCGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCtt                                                                                                                   <  2:90893/149‑1 (MQ=255)
                                         ctcGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTg                                                                                                          <  3:79657/149‑1 (MQ=255)
                                            ggCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTg                                                                                                          <  1:249931/146‑1 (MQ=255)
                                            ggCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGGCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGt                                                                                                        <  5:38215/148‑1 (MQ=255)
                                                    tCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCa                                                                                                                                                 <  8:223665/99‑1 (MQ=255)
                                                      aTTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTACCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCca                                                                                             >  2:282885/1‑149 (MQ=255)
                                                      aTTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTACCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCca                                                                                             >  2:287662/1‑149 (MQ=255)
                                                             cTTAGCGTTTCCTGGGAGACCGAGTTGGCGTCCAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGAGCCCCCACCAGAAGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATAGTc                                                                                      <  2:109203/149‑1 (MQ=255)
                                                              ttAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCa                                                                                     >  8:279311/1‑149 (MQ=255)
                                                              ttAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCa                                                                                     <  8:27071/149‑1 (MQ=255)
                                                                 gCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGc                                                                                  >  1:229362/1‑149 (MQ=255)
                                                                      tCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCt                                                                                 >  4:167682/1‑145 (MQ=255)
                                                                      tCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCt                                                                                 <  3:167682/145‑1 (MQ=255)
                                                                       ccTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCa                                                                            >  8:346819/1‑149 (MQ=255)
                                                                          gggTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGc                                                                         >  2:114585/1‑149 (MQ=255)
                                                                                 cGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTc                                                                                      <  2:260156/129‑1 (MQ=255)
                                                                                       ggCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGa                                                                       <  7:122597/138‑1 (MQ=255)
                                                                                       ggCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGa                                                                       >  8:122597/1‑138 (MQ=255)
                                                                                       ggCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAgg                                                             <  7:346819/148‑1 (MQ=255)
                                                                                       ggCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTACCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAgg                                                             <  1:287662/148‑1 (MQ=255)
                                                                                       ggCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTACCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAgg                                                             <  1:282885/148‑1 (MQ=255)
                                                                                          gTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGa                                                         <  2:229362/149‑1 (MQ=255)
                                                                                                   cACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGac                                                                      <  3:334255/127‑1 (MQ=255)
                                                                                                   cACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGac                                                                      >  4:334255/1‑127 (MQ=255)
                                                                                                    aCGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGc                                                 <  5:49303/147‑1 (MQ=255)
                                                                                                              tGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTgg                                      <  8:197599/148‑1 (MQ=255)
                                                                                                                   gACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGc                                <  1:114585/149‑1 (MQ=255)
                                                                                                                        tAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCtt                             <  1:195337/147‑1 (MQ=255)
                                                                                                                        tAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGc                           >  1:169817/1‑149 (MQ=255)
                                                                                                                               gtccccACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGc                    >  4:12791/3‑149 (MQ=255)
                                                                                                                                 ccccACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGcc                   >  1:25724/1‑148 (MQ=255)
                                                                                                                                  cccACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTg                            <  3:247130/138‑1 (MQ=255)
                                                                                                                                  cccACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTg                            >  4:247130/1‑138 (MQ=255)
                                                                                                                                   ccaccaGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGaa                >  7:175925/1‑149 (MQ=255)
                                                                                                                                     accaGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTCTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAAcc              >  4:313957/1‑149 (MQ=255)
                                                                                                                                       caGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCgg            >  6:334841/1‑149 (MQ=255)
                                                                                                                                       caGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCgg            >  1:234659/1‑149 (MQ=255)
                                                                                                                                         gATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCGGcc          <  2:25724/149‑1 (MQ=255)
                                                                                                                                                gTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCGGCCTTTTtct   <  7:279311/149‑1 (MQ=255)
                                                                                                                                                 ttGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCGGCCTTTTtctc  <  3:12791/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCCAGATCCACCGGTGAGTTGGACGCATCCAACGACAGGTCTCGGCGATAGTCATTGCCGCTTAGCGTTTCCTGGGTGACCGAGTTGGCGTCGAGCCTCACGTCGCCTTTGGCCGACTGTAAGCGTGCCCCCACCAGATGAGCGTTGCCACCGACCTGCAGGTCGATGCCGTCGCTGCTTTTAAGCGTGGTCTGACGCGCCACATTGTCACGCTGTGCATGCGACACGTCGAGGCTGACGGTCGGCGAGAACCCTGGGTCTGCCTTGCTCAGGGCCGAACCGGCCTTTTTCTC  >  NC_002947/1652377‑1652670

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