Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,241,025 C→T 31.5% V240M (GTG→ATG)  PP_3718 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,241,0250CT31.5% 20.9 / 10.2 19V240M (GTG→ATG) PP_3718aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (4/2);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTG  >  NC_002947/4240889‑4241167
                                                                                                                                        |                                                                                                                                              
ccAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                     >  1:558911/1‑147 (MQ=255)
ccAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                     <  2:558911/148‑2 (MQ=255)
      aGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtc                                                                                                                             <  2:417794/150‑1 (MQ=255)
                    gaagaaGGTGCCCAGGCGCCAGAACGAGGACACCACTTCGTCGGCGACGAAGAGGATGTCGAAGATATAGGACAGCTGCAACATGCGCTGAAAGTAGCCCTCGGGTGGAATGATCATCCCGCCCGACCGCATGATTGGCTCGGTgatgaa                                                                                                               >  2:66494/1‑150 (MQ=255)
                                                 tCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                  >  1:158622/1‑150 (MQ=255)
                                                           gtcgGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTc                                                                        <  2:190955/150‑1 (MQ=255)
                                                                            tGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGccacca                                                       <  1:63713/150‑1 (MQ=255)
                                                                                     tCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGa                                              >  1:480240/1‑150 (MQ=255)
                                                                                          cacaGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAg                                                                                                                                      >  1:299372/1‑57 (MQ=18)
                                                                                          cacaGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAg                                                                                                                                      <  2:299372/57‑1 (MQ=18)
                                                                                                    acaTGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGt                               >  1:302895/1‑150 (MQ=255)
                                                                                                              gAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgtcgtc                     <  2:158622/150‑1 (MQ=255)
                                                                                                                     cccTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGgcgcg              >  2:74197/1‑150 (MQ=255)
                                                                                                                      gcTCGGGCGGGGTGATCATTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGgcgcgc             >  2:309883/2‑150 (MQ=255)
                                                                                                                       cTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGcgcgcg            >  1:263984/1‑150 (MQ=255)
                                                                                                                        tCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCgg           <  1:74197/150‑1 (MQ=255)
                                                                                                                                aaTGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGACCTCGGc                                                            <  2:614884/93‑1 (MQ=255)
                                                                                                                                aaTGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGc                                                            >  1:614884/1‑93 (MQ=255)
                                                                                                                                 aTGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttg  <  1:498739/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                              
CCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTG  >  NC_002947/4240889‑4241167

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: