Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,184 G→A 11.1% T308T (ACG→ACA)  PP_3662 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,1840GA11.1% 91.4 / 4.1 36T308T (ACG→ACA) PP_3662hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base A (2/2);  total (18/18)
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.73e-01

GCACGGCGGAGGAGTTCCTGTACCTGCTGGGCATCCTCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGCCTTCCACTTCAACTGGCTGCTGAA  >  NC_002947/4160045‑4160327
                                                                                                                                           |                                                                                                                                               
gCACGGCGGAGGAGTTCCTGTACCTGCTGGGCATCCTCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTg                                                                                                                                          >  1:275639/1‑147 (MQ=255)
     gcggAGGAGTTCCTGTACCTGCTGGGCATCCTCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGc                                                                                                                                   <  1:140616/149‑1 (MQ=255)
                 cTGTACCTGCTGGGCATCCTCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTa                                                                                                                       >  5:114229/1‑149 (MQ=255)
                             ggCATCCTCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAAtcatca                                                                                                            <  2:275639/148‑1 (MQ=255)
                              gCATCCTCATGCACCCGGACAAGCACGACGTGCCGTTCCCGGTGTTGCTCACAGGCCCTCGCAGCGCCGAGCCTTATCTGCAAGAGTTGCAGGATTTCTTCTGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCttt                                                                                                          <  6:114229/149‑1 (MQ=255)
                                    tCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACACTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGAc                                                                                                    >  6:331767/1‑149 (MQ=255)
                                    tCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGATCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACACTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGAc                                                                                                    >  6:327270/1‑149 (MQ=255)
                                            ccGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCcggcgg                                                                                             >  5:14416/1‑148 (MQ=255)
                                             cGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGaa                                                                                           <  3:22340/149‑1 (MQ=255)
                                                    ccACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGccc                                                                                    >  3:281033/1‑149 (MQ=255)
                                                            tGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATg                                                                            >  5:120223/1‑149 (MQ=255)
                                                                 ttCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCg                                                                        >  4:278652/1‑148 (MQ=255)
                                                                    ccGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGcac                                                                               >  6:25044/1‑138 (MQ=255)
                                                                    ccGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGcac                                                                               <  5:25044/138‑1 (MQ=255)
                                                                    ccGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAgg                                                                    <  3:278652/149‑1 (MQ=255)
                                                                     cGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATg                                                                            <  8:35746/140‑1 (MQ=255)
                                                                     cGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATg                                                                            >  7:35746/1‑140 (MQ=255)
                                                                     cGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGc                                                                   >  8:256042/1‑149 (MQ=255)
                                                                        tggtgCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTc                                                                <  4:281033/149‑1 (MQ=255)
                                                                           tgCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGcaca                                                                              <  8:239144/132‑1 (MQ=255)
                                                                           tgCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGcaca                                                                              >  7:239144/1‑132 (MQ=255)
                                                                                aCCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGt                                                        >  1:130337/1‑149 (MQ=255)
                                                                                aCCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACCCTTTCCTCGGCCACAAGCTGGGTGAACCGGACCATCACCTGTACGAAATCATCATTGATGACCAGGCGTAACTCCCACTGCACATGGACGAAGGCCTCAAACCgg                                                         >  4:94557/1‑148 (MQ=255)
                                                                                    gCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAg                                                                                          >  5:259010/1‑111 (MQ=255)
                                                                                    gCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAg                                                                                          <  6:259010/111‑1 (MQ=255)
                                                                                      ccGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAg                                                    >  3:204484/1‑147 (MQ=255)
                                                                                                         aTCTGCAACAGTTGCACGCTTTCGTCGGCGCCACACTGTGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCTGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCaa                                <  5:327270/148‑1 (MQ=255)
                                                                                                         aTCTGCAACAGTTGCACGCTTTCGTCGGCGCCACACTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCaa                                <  5:331767/148‑1 (MQ=255)
                                                                                                          tCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAAc                               <  4:204484/148‑1 (MQ=255)
                                                                                                            tGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGAt                            >  2:227967/1‑149 (MQ=255)
                                                                                                              cAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGt                                                <  3:62804/127‑1 (MQ=255)
                                                                                                               aaCAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGcc                         >  2:28489/1‑149 (MQ=255)
                                                                                                                  aGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGCCTTc                      <  7:256042/149‑1 (MQ=255)
                                                                                                                      gCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGCCTTCCAc                   <  5:192708/148‑1 (MQ=255)
                                                                                                                          gCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGCCTTCCACTTCaa              <  6:120223/149‑1 (MQ=255)
                                                                                                                                      gcCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGCCTTCCACTTCAACTGGCTGCTGaa  <  1:227967/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GCACGGCGGAGGAGTTCCTGTACCTGCTGGGCATCCTCATGCACCCGGACAACCACGACCTGCCGTTCCCGGTGGTGCTCACCGGCCCGCGCAGCGCCGAGCCTTATCTGCAACAGTTGCACGCTTTCGTCGGCGCCACGCTGGGTGAAGCGGCCCATCGCCTGTACGAAATCATCATTGATGACCCGGCGGAAGTCGCCCGGCACATGGTCGAAGGCCTCAAAGCGGTCAAGCAGTTCCGCCGCGAGCGCAACGATGCCTTCCACTTCAACTGGCTGCTGAA  >  NC_002947/4160045‑4160327

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