Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,241,026 C→T 26.4% G239G (GGG→GGA PP_3718 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,241,0260CT26.4% 25.9 / 8.1 19G239G (GGG→GGAPP_3718aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (3/2);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGA  >  NC_002947/4240889‑4241175
                                                                                                                                         |                                                                                                                                                     
ccAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                             <  2:558911/148‑2 (MQ=255)
ccAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                             >  1:558911/1‑147 (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)
                                                                                                                                         cccgcccgAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGATCTCGGCCACCAGGAAGTCGAGGCCATCTGCCAAGTCCATTTCTGCTGGCGCGCGGTACGGGTTGGGGTTGga  >  1:598700/1‑150 (MQ=255)
                                                                                                                                         cccgcccgAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTg                           >  2:211705/1‑125 (MQ=255)
                                                                                                                                         cccgcccgAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTg                           <  1:211705/125‑1 (MQ=255)
                                                                                                                                         cccgcccgAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgg         <  2:361950/143‑1 (MQ=255)
                                                                                                                                         cccgcccgAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttgg         >  1:361950/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
CCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGA  >  NC_002947/4240889‑4241175

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: