Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,041,183 T→A 15.9% Q17H (CAA→CAT PP_5293 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,041,1830TA15.9% 57.1 / 3.8 25Q17H (CAA→CATPP_5293hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base A (2/2);  total (12/13)
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.85e-01

GAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCG  >  NC_002947/6041041‑6041329
                                                                                                                                              |                                                                                                                                                  
gAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCgaa                                                                                                                                             >  2:336231/1‑148 (MQ=255)
gAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCgaa                                                                                                                                             >  2:304043/1‑148 (MQ=255)
             ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATggg                                                                                                                                <  2:97379/150‑1 (MQ=255)
             ctctTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCAATGCCGAAATCCGGCAATggg                                                                                                                                <  1:336231/150‑1 (MQ=255)
                  tGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCgggg                                                                                                                           <  2:508332/150‑1 (MQ=255)
                         ggCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGa                                                                                                                    <  2:209963/150‑1 (MQ=255)
                            gTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAg                                                                                                                 <  2:464702/150‑1 (MQ=255)
                                    aaGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAAt                                                                                                         >  2:354801/1‑150 (MQ=255)
                                      gggTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGTTGATCGACAGGTTCGAATCg                                                                                                       >  2:436853/1‑150 (MQ=255)
                                              cTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCAAGAGGGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGcc                                                                                               >  1:163216/1‑150 (MQ=255)
                                                        aGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACAGAGGCGATTGCCGGATTTCGGCATTTGGCGGGGTGATCGACAGGTTCGAATCACTCATGCCCCGCAactcc                                                                                     >  1:568202/1‑150 (MQ=255)
                                                                        gCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGATACAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcagtc                                                                     >  1:162831/1‑150 (MQ=255)
                                                                           ggCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTtcagtcaac                                                                  >  2:403273/1‑148 (MQ=255)
                                                                           ggCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGc                                                                  <  2:158110/150‑1 (MQ=255)
                                                                              ccAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAg                                                               >  1:518258/1‑150 (MQ=255)
                                                                                                      ccAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGa                                                                                                                        >  1:108858/1‑69 (MQ=255)
                                                                                                      ccAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGa                                                                                                                        <  2:108858/69‑1 (MQ=255)
                                                                                                       cAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtc                                      <  2:518258/150‑1 (MQ=255)
                                                                                                         ggCCCAACAGCGCGGCAAGGTGGACGCTCTTTCCCTATGGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAAtctc                                    <  2:162831/150‑1 (MQ=255)
                                                                                                            ccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgct                                                           >  2:216954/1‑124 (MQ=255)
                                                                                                            ccAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGtgct                                                           <  1:216954/124‑1 (MQ=255)
                                                                                                                 aGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGc                            <  2:266923/150‑1 (MQ=255)
                                                                                                                 aGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGc                            <  1:354801/150‑1 (MQ=255)
                                                                                                                               gCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCaaa              <  2:175229/150‑1 (MQ=255)
                                                                                                                                           ccATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGGTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCg  >  2:79419/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GAAAGTTCAGCGCCTCTTTGCCCTGGGCGTCGATAAAAGGGTCTTGCTCGCGCTTGAGGTGCCATAGCACCAGCGGGCCCAGCAGGTGCCCAAGCGGCACTACCAGGCCCAACAGCGCGGAGAGGTGGCAGAACATCGCCCATTGCCGGATTTCGGCATTGGGCGGGGTGATCGACAGGTTCGAATCGCTCATGCCCCGCACCTCCTTGGCCGGTTTCAGTCAGCCAGTGCTGCCTGCTGCAGCTCGAAAATCTCGTTCATGCCCTTCTGCGCCAAAGCGAGCATGGCG  >  NC_002947/6041041‑6041329

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: