Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,875,244 G→A 23.5% V129V (GTG→GTA PP_4284 → transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,875,2440GA23.5% 22.8 / 4.3 17V129V (GTG→GTAPP_4284transporter
Reads supporting (aligned to +/- strand):  ref base G (4/9);  new base A (2/2);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

GCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGCCGGCGTCGGATTGTTTCTCGGGCTGATCGGCCTGAAGACTGCAGGCA  >  NC_002947/4875102‑4875374
                                                                                                                                              |                                                                                                                                  
gcctgcTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCa                                                                                                                             <  1:348731/150‑1 (MQ=255)
 cctgcTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTAAAGCGTTGGCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCAt                                                                                                                            >  2:356488/1‑150 (MQ=255)
               tGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCNCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGt                                                                                                              <  1:9407/150‑1 (MQ=255)
                  aCGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGa                                                                                                           <  1:573829/150‑1 (MQ=255)
                   cGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGaa                                                                                                          >  2:413488/1‑150 (MQ=255)
                                    gCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGgctgc                                                                                         <  1:526444/150‑1 (MQ=255)
                                                    ggggCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGt                                                                         <  1:413488/150‑1 (MQ=255)
                                                            agggcttctgTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTgc                                                                 <  2:373064/146‑1 (MQ=255)
                                                                         cTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTATTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATggg                                                    >  1:320221/1‑150 (MQ=255)
                                                                         cTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTATTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATggg                                                    >  1:320204/1‑150 (MQ=255)
                                                                             cccGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGcc                                                <  1:3995/149‑1 (MQ=255)
                                                                                                                       ggggCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGaa                                                                                                <  1:422754/60‑1 (MQ=255)
                                                                                                                       ggggCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGaa                                                                                                >  2:422754/1‑60 (MQ=255)
                                                                                                                           cGGTATTCGTCTCGGGCGTATTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGCCGGCGTCGGATTGTTTCTCGGGCTGATCGGCCTGAAGACTGCAGGca  <  2:320204/150‑1 (MQ=255)
                                                                                                                           cGGTATTCGTCTCGGGCGTATTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGCCGGCGTCGGATTGTTTCTCGGGCTGATCGGCCTGAAGACTGCAGGca  <  2:320221/150‑1 (MQ=255)
                                                                                                                                      tcGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGCCGGCGTCGGATTGTTTCTCGGGCTGATCGGCCTGAAGACTGCAGGc   <  2:362714/138‑1 (MQ=255)
                                                                                                                                      tcGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGCCGGCGTCGGATTGTTTCTCGGGCTGATCGGCCTGAAGACTGCAGGc   >  1:362714/1‑138 (MQ=255)
                                                                                                                                              |                                                                                                                                  
GCCTGCTGATGGGCCTGTACGCCAACTGGCCGGTGGGCCTGGCGCCGGGCATGGGGCTCAATGCCTTCTTTACCTACACCGTGGTCGGCACCATGGGCTACAACTGGGAAACGGCGCTGGGGGCGGTATTCGTCTCGGGCGTGTTGTTCATGTTCCTGACCTTGTCGAAAGTGCGCGAATGGCTGCTCAACAGCATTCCGGTCAGCCTGCGCCATGCCATGGGGGCCGGCGTCGGATTGTTTCTCGGGCTGATCGGCCTGAAGACTGCAGGCA  >  NC_002947/4875102‑4875374

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: