Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 97,931 G→T 25.0% V30V (GTC→GTA PP_0091 ← putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294797,9310GT25.0% 34.7 / 3.2 16V30V (GTC→GTAPP_0091putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base T (2/2);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GTAGATTTCATAATCGAAGCGCACCACCAGTTCATGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCGATCAGCATCATGGGCTGGCGCATGGCAAAACTCCCACTTGCGATTTACAAGT  >  NC_002947/97796‑98049
                                                                                                                                       |                                                                                                                      
gTAGATTTCATAATCGAAGCGCACCACCAGTTCATGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCg                                                                                                                        <  1:401631/136‑1 (MQ=255)
gTAGATTTCATAATCGAAGCGCACCACCAGTTCATGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCg                                                                                                                        <  1:253774/136‑1 (MQ=255)
     tttCATAATCGAAGCGCACCACCAGTTCATGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCa                                                                                                                   <  2:26573/136‑1 (MQ=255)
                       ccaccaGTTCATGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGAc                                                                                                 <  1:180575/136‑1 (MQ=255)
                            aGTTCAAGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGACGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGCCCTGCTAGGGGTCGACGggcg                                                                                            >  1:13663/1‑136 (MQ=255)
                                     ttGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCTACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGc                                                                                   >  1:285219/1‑136 (MQ=255)
                                     ttGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCTACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGc                                                                                   >  1:454806/1‑136 (MQ=255)
                                              gtCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAg                                                                          <  1:217613/136‑1 (MQ=255)
                                                   cTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGccccc                                                                     <  1:203290/136‑1 (MQ=255)
                                                          tACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTAc                                                              >  1:479448/1‑136 (MQ=255)
                                                                ccATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAg                                                                          <  1:546357/118‑1 (MQ=255)
                                                                ccATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAg                                                                          >  2:546357/1‑118 (MQ=255)
                                                                 cATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGTTCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGc                                                       <  1:29118/136‑1 (MQ=255)
                                                                   tCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCGa                                                     <  1:100442/136‑1 (MQ=255)
                                                                                                       aTCACCCGACCGGGGGTCATGGTGTACAGATCTACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCGATCAGCATCATGGGCTGGCGCATGGCAAAACTCCCAc                 <  2:285219/136‑1 (MQ=255)
                                                                                                       aTCACCCGACCGGGGGTCATGGTGTACAGATCTACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCGATCAGCATCATGGGCTGGCGCATGGCAAAACTCCCAc                 <  2:454806/136‑1 (MQ=255)
                                                                                                       aTCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCGATCAGCATCATGGGCTGGCGCATGGCAAAACTCCCAc                 <  1:412344/136‑1 (MQ=255)
                                                                                                                      gTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCTATCAGCATCATGGGCTGGCGCATGGCAAAACTCCCACTTGCGATTTACAAGt  <  1:237068/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                      
GTAGATTTCATAATCGAAGCGCACCACCAGTTCATGCTTGCCCGGTGTCACCTGAAAGTACCGGCCATCGTCCAGGCGCTTGCCGTCGAGGCGGTCGGCCATGATCACCCGACCGGGGGTCATGGTGTACAGATCGACCCAGGCCTGCTTGGGGTCGACGGGCGGCAAGGGGCTGGCGCAAGCCCCCAGTGTACTGAGGGCGATCAGCATCATGGGCTGGCGCATGGCAAAACTCCCACTTGCGATTTACAAGT  >  NC_002947/97796‑98049

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: