Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,478,830 G→C 17.4% I33M (ATC→ATG PP_4816 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,478,8300GC17.4% 54.0 / 6.6 23I33M (ATC→ATGPP_4816hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/11);  new base C (2/2);  total (10/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.29e-01

TTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGG  >  NC_002947/5478690‑5478965
                                                                                                                                            |                                                                                                                                       
ttGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCg                                                                                                                                >  1:118476/1‑150 (MQ=255)
                    gACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg                                                                                                            <  1:443630/150‑1 (MQ=255)
                              ggcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc                                                                                                  <  1:279018/150‑1 (MQ=255)
                              ggcACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATc                                                                                                  <  1:102152/150‑1 (MQ=255)
                                 aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctgtt                                                                                                                              >  1:602673/1‑115 (MQ=255)
                                 aCCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATCATCGGcctgtt                                                                                                                              <  2:602673/119‑5 (MQ=255)
                                   cagcagcagCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg                                                                                                            >  1:128684/1‑135 (MQ=255)
                                   cagcagcagCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAg                                                                                                            <  2:128684/135‑1 (MQ=255)
                                                                         ccTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCAc                                                       <  2:515254/150‑1 (MQ=255)
                                                                                 aGCGGCTTGCCGACCTAACCGCGGGCGATGAGGAAGAGCGGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGc                                               >  2:441858/1‑150 (MQ=255)
                                                                                     gCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGcgc                                                                                                                                 <  2:352431/64‑3 (MQ=255)
                                                                                     gCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGACGATGATCGGcgc                                                                                                                                 >  1:352431/1‑62 (MQ=255)
                                                                                           cGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa                                           >  1:496671/1‑144 (MQ=255)
                                                                                           cGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATaa                                           <  2:496671/144‑1 (MQ=255)
                                                                                                  aCCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCAc                                                       <  2:410314/125‑1 (MQ=255)
                                                                                                  aCCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCAc                                                       >  1:410314/1‑125 (MQ=255)
                                                                                                   ccGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAACCGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCggg                             <  1:339106/150‑1 (MQ=255)
                                                                                                                    gagCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATaa            >  1:49028/1‑150 (MQ=255)
                                                                                                                        agcagcaACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCAc        <  2:118476/150‑1 (MQ=255)
                                                                                                                              aacaacAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCgg  <  1:111549/150‑1 (MQ=255)
                                                                                                                              aacaacAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCAAGCTCgg  >  1:588932/1‑150 (MQ=255)
                                                                                                                                  ccgggCCGATCATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCAt                          >  1:75007/4‑122 (MQ=255)
                                                                                                                                  ccgggCCGATCATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCAt                          <  2:75007/119‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
TTGGCATATACCATCACCCCGACTGCCGGCGGCACCAGCAGCAGCGGGAGATAGCGCAGCAAGCTGCTTGCGGCCTCGTTCAGCGGCTTGCCGACCTCACCGCGGGCCATGAGGAAGAGCAGCAGCAACAACAGGCCGATGATCGGCCCGGGTAAAAACGGCACGAACAGGTGATTGATCGCCGTCCCCAGCAGCTGGAACAGTACCAGCCAGGTCAAACCACGCAACAGCATAAGCATCTCCCTCGGGCATGGCCGCCATTATAAGCACGCTCGG  >  NC_002947/5478690‑5478965

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