Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 874,543 Δ1 bp 8.7% coding (119/840 nt) PP_0759 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947874,5430A. 8.7% 130.5 / 7.7 46coding (119/840 nt)PP_0759hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/22);  new base . (2/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

GACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGC  >  NC_002947/874403‑874690
                                                                                                                                             |                                                                                                                                                   
gACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACa                                                                                                                                            >  1:3656/1‑150 (MQ=255)
            gATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAg                                                                                                                                <  1:28433/150‑1 (MQ=255)
                 ggCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTAc                                                                                                                                        >  2:13202/1‑137 (MQ=255)
                 ggCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTAc                                                                                                                                        <  1:13202/137‑1 (MQ=255)
                 ggCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGc                                                                                                                           >  1:469304/1‑150 (MQ=255)
                     cATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCCGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTCcatcatcggcgc                                                                                                                       <  2:612795/150‑7 (MQ=255)
                      aTGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCg                                                                                                                      >  2:206059/1‑150 (MQ=255)
                       tgtgGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTTGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGa                                                                                                                     <  1:513660/150‑1 (MQ=255)
                           gATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGAGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGa                                                                                                                 >  1:272758/1‑150 (MQ=255)
                            aTAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                                <  2:656535/150‑1 (MQ=255)
                            aTAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                                <  1:58317/150‑1 (MQ=255)
                               gggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                             <  2:469304/150‑1 (MQ=255)
                               gggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAg                                                                                                             <  2:3656/150‑1 (MQ=255)
                                ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTNCTGAGTGCTNGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAgg                                                                                                            >  1:395786/1‑150 (MQ=255)
                                ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCaa                                                                                                                                              >  1:355583/1‑116 (MQ=255)
                                ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCaa                                                                                                                                              <  2:355583/116‑1 (MQ=255)
                                ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACa                                                                                                                                       <  2:218806/123‑1 (MQ=255)
                                ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACa                                                                                                                                       >  1:218806/1‑123 (MQ=255)
                                       tttCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTCCAGCATGCCGACTACCCATGAGGGCAAccg                                                                                                     <  2:503343/150‑2 (MQ=255)
                                         tCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGa                                                                                                   <  1:127377/150‑1 (MQ=255)
                                           aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTACTGTTGGCGTTCGACGCCAAc                                                                                                                                             >  2:208919/1‑106 (MQ=255)
                                           aTGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTACTGTTGGCGTTCGACGCCAAc                                                                                                                                             <  1:208919/106‑1 (MQ=255)
                                                tAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGAGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTAcgc                                                                                            <  2:272758/150‑1 (MQ=255)
                                                        ttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGt                                                                                                <  2:729414/138‑1 (MQ=255)
                                                        ttGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGt                                                                                                >  1:729414/1‑138 (MQ=255)
                                                              tcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATcagc                                                                              >  1:303377/1‑150 (MQ=255)
                                                                  ccTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAc                                                                          >  2:100374/1‑150 (MQ=255)
                                                                     gctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCCGTGGTATCAGCAAACCCt                                                                       >  2:581599/1‑150 (MQ=255)
                                                                              tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGc                                                              >  2:237170/1‑150 (MQ=255)
                                                                               cgCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCa                                                             <  2:395786/149‑1 (MQ=255)
                                                                                   tCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGa                                                         >  2:623660/1‑150 (MQ=255)
                                                                                                        cGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGt                                    <  1:40619/150‑1 (MQ=255)
                                                                                                            gAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAg                                >  2:210626/1‑150 (MQ=255)
                                                                                                                 ccaccaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAg                                                                                        >  1:112594/1‑89 (MQ=255)
                                                                                                                 ccaccaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAg                                                                                        <  2:112594/89‑1 (MQ=255)
                                                                                                                              tGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTGTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAACCTGGAGGAAGTCTCGACCAAGGGGGTTCGCGACAAGGCCAGCGCCGATGCCAAGCgtgt              >  2:611531/1‑150 (MQ=255)
                                                                                                                               gCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGtgtg             <  2:818690/150‑1 (MQ=255)
                                                                                                                                 tGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGaa                                                                <  2:160254/97‑1 (MQ=255)
                                                                                                                                 tGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGaa                                                                >  1:160254/1‑97 (MQ=255)
                                                                                                                                 tGTTGGCGTTGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCa           >  2:583472/1‑150 (MQ=255)
                                                                                                                                 tGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGCGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCACGGGTGTGCTCGACAACTCCAGCGCCGATGCCAAGCGTGTGCa           >  2:536279/1‑150 (MQ=255)
                                                                                                                                  gTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAg          <  1:100374/150‑1 (MQ=255)
                                                                                                                                   ttGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAgg         <  1:210626/150‑1 (MQ=255)
                                                                                                                                           cGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGTTGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGAAACGGGGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGc  <  1:581599/150‑1 (MQ=255)
                                                                                                                                            gAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgtgt              <  1:606690/137‑1 (MQ=255)
                                                                                                                                            gAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgtgt              >  2:606690/1‑137 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTGC  >  NC_002947/874403‑874690

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: