Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,972,776 C→T 9.3% T22M (ACG→ATG)  PP_5725 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,972,7760CT 9.3% 119.0 / 3.0 43T22M (ACG→ATG) PP_5725hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base T (2/2);  total (22/21)
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.27e-01

CGGCAGCGCGGGCGGCAACGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACAAGTCATGAAGAATTTTTAA  >  NC_002947/5972639‑5972913
                                                                                                                                          |                                                                                                                                         
cGGCAGCGCGGGCGGCAACGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGcc                                                                                                                                        >  5:3041/1‑141 (MQ=255)
cGGCAGCGCGGGCGGCAACGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGcc                                                                                                                                        <  6:3041/141‑1 (MQ=255)
       gcgggcggCAACGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCa                                                                                                                         <  5:234640/149‑1 (MQ=255)
        cgggcggCAACGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAg                                                                                                                        <  3:76646/149‑1 (MQ=255)
                aaCGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCAGAGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGc                                                                                                                  >  1:253872/1‑148 (MQ=255)
                aaCGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCGTGACGGCCAATCAGCGCGGc                                                                                                                  >  5:289873/1‑147 (MQ=255)
                    tCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGttt                                                                                                            <  6:289873/149‑1 (MQ=255)
                        gggCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGAtg                                                                                                        >  6:17308/1‑149 (MQ=255)
                           cGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATgtgt                                                                                                     <  5:17308/149‑1 (MQ=255)
                            ggTGGGTACCGTATCAACAACAAGAAGTATGCGTGAGTTGCGACCA‑TGTCCGAGGCACGCATCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGtt                                                                                                    <  4:22792/149‑1 (MQ=255)
                                gggACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCa                                                                                                                                      >  6:144457/1‑111 (MQ=255)
                                gggACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCa                                                                                                                                      <  5:144457/111‑1 (MQ=255)
                                         tCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATc                                                                                       >  7:295836/1‑149 (MQ=255)
                                           aacaacaaGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATgcgc                                                                                            >  6:220565/1‑142 (MQ=255)
                                           aacaacaaGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATgcgc                                                                                            <  5:220565/142‑1 (MQ=255)
                                              aacaaGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCaa                                                                                  >  4:316165/1‑149 (MQ=255)
                                                  agaagTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAgt                                                                              >  2:89297/1‑149 (MQ=255)
                                                                   gCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAggtggt                                                               <  6:34279/147‑1 (MQ=255)
                                                                   gCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAggtggt                                                               >  5:34279/1‑147 (MQ=255)
                                                                   gCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGgtggtg                                                              >  3:168269/1‑148 (MQ=255)
                                                                      aCCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACt                                                          <  4:168269/149‑1 (MQ=255)
                                                                       ccA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATTACGGACAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTg                                                         >  2:188215/1‑149 (MQ=255)
                                                                       ccA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATgcgc                                                                                            >  5:155266/1‑114 (MQ=255)
                                                                       ccA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATgcgc                                                                                            <  6:155266/114‑1 (MQ=255)
                                                                             tCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTa                                                     <  1:173549/148‑1 (MQ=255)
                                                                               cGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGtt                                                                                                    <  7:71644/99‑1 (MQ=255)
                                                                               cGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGtt                                                                                                    >  8:71644/1‑99 (MQ=255)
                                                                                     aCGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGt                                             >  3:103083/1‑148 (MQ=255)
                                                                                          gCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCa                                       <  4:103083/149‑1 (MQ=255)
                                                                                          gCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCa                                       <  1:89297/149‑1 (MQ=255)
                                                                                            ttCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGa                                                                                                          >  2:181887/1‑80 (MQ=255)
                                                                                            ttCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGa                                                                                                          <  1:181887/80‑1 (MQ=255)
                                                                                              cATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTGGGTGATCCAAcc                                    <  6:192740/148‑1 (MQ=255)
                                                                                                ttAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACg                                 >  2:34602/1‑149 (MQ=255)
                                                                                                             gCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACaa                    <  3:316165/149‑1 (MQ=255)
                                                                                                                       gATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACAAGTCATgaag           >  4:209442/1‑148 (MQ=255)
                                                                                                                             cTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACAAGTCATGAAGAAtttt     >  3:8832/1‑148 (MQ=255)
                                                                                                                                   gCCGATATGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGgtggtg                                                              <  5:276452/85‑1 (MQ=255)
                                                                                                                                   gCCGATATGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGgtggtg                                                              >  6:276452/1‑85 (MQ=255)
                                                                                                                                   gCCGATATGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGgtggtg                                                              <  7:220669/85‑1 (MQ=255)
                                                                                                                                   gCCGATATGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGgtggtg                                                              >  8:220669/1‑85 (MQ=255)
                                                                                                                                     cGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACAAGTCATGAAGAATTTTTaa  >  2:56200/1‑143 (MQ=255)
                                                                                                                                     cGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACAAGTCATGAAGAATTTTTaa  <  1:56200/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CGGCAGCGCGGGCGGCAACGTCACGGGCGGTTGGGACCGTATCAACAACAAGAAGTATGCGTGAGGTGCGACCA‑TGTCCGAGGCACGCAGCTTCATTAATTCGCACACGCAGTCATATGATTCGCTGAAGGCCGATACGCCCATGACGGCCAATCAGCGCGGCAAGTTTGATGTGTTGAATGCGCATATCGTCAATAGTGTGGTGGTAGGTGGTGAACTGGTTATTGTAGGTGATCCAACCACGCCCTCCTGTACAAGTCATGAAGAATTTTTAA  >  NC_002947/5972639‑5972913

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