Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 431,365 C→T 12.5% V318I (GTT→ATT)  PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947431,3650CT12.5% 100.8 / 6.0 40V318I (GTT→ATT) PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (15/20);  new base T (3/2);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

GGTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCGCTCGGTCA  >  NC_002947/431218‑431514
                                                                                                                                                   |                                                                                                                                                     
ggTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCa                                                                                                                                                     <  1:306042/150‑1 (MQ=255)
    gCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCTTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGt                                                                                                                                                 >  2:21136/1‑150 (MQ=255)
                     ccagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCa                                                                                                                                            <  2:327799/138‑1 (MQ=255)
                     ccagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCa                                                                                                                                            >  1:327799/1‑138 (MQ=255)
                      cagTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATg                                                                                                                               <  2:61512/150‑1 (MQ=255)
                           cGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTa                                                                                                                          <  1:210845/150‑1 (MQ=255)
                                    ggCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTgcgcg                                                                                                                 <  2:602980/150‑1 (MQ=255)
                                    ggCACACGGTATGGGCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCCTGGGGTATGGTgcgcg                                                                                                                 >  2:177454/1‑150 (MQ=255)
                                              aTGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAAc                                                                                                       <  1:126470/150‑1 (MQ=255)
                                                  tCGTTGCGCCACGTGATGAGCTGGGTGATCTGTGGGGACGATGAGCCATGTGCCAGGATGCGGTTGACGAGGTGGCTGAGTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGGACCAGa                                                                                                   <  2:68198/150‑1 (MQ=255)
                                                                 aTGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAgg                                                                                    >  2:259365/1‑150 (MQ=255)
                                                                 aTGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAgg                                                                                    <  1:259365/150‑1 (MQ=255)
                                                                   gATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGga                                                                                  <  1:304159/150‑1 (MQ=255)
                                                                    aTCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGgag                                                                                 <  1:21136/150‑1 (MQ=255)
                                                                       tGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAAc                                                                              <  1:489874/150‑1 (MQ=255)
                                                                         ggTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCa                                                                                                     <  1:138395/125‑1 (MQ=255)
                                                                         ggTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCa                                                                                                     >  2:138395/1‑125 (MQ=255)
                                                                                       acgaTGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGaaa                                                              <  2:353239/150‑1 (MQ=255)
                                                                                         gaTGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTc                                                                         >  1:237223/1‑137 (MQ=255)
                                                                                         gaTGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTc                                                                         <  2:237223/137‑1 (MQ=255)
                                                                                               ccATGGGCCAGCACGAGCAACACAAGCTGGCTGATTTCGCCGGCCAGCGAAACCAGCATATCCAGGCTTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGGTGCAGGGAGAAGAGGTCGCGCTCGGAAACCAcgccg                                                      >  2:6752/1‑150 (MQ=255)
                                                                                                    ggCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAg                                                 <  2:481808/150‑1 (MQ=255)
                                                                                                    ggCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAg                                                 <  2:397334/150‑1 (MQ=255)
                                                                                                     gCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAgg                                                >  2:534866/1‑150 (MQ=255)
                                                                                                     gCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAgg                                                >  2:295346/1‑150 (MQ=255)
                                                                                                            tGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGt                                         >  1:454632/1‑150 (MQ=255)
                                                                                                                ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCAt                                                                                                                                >  2:254295/1‑59 (MQ=255)
                                                                                                                ttCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCAt                                                                                                                                <  1:254295/59‑1 (MQ=255)
                                                                                                                   tACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTccac                                  >  2:151403/1‑150 (MQ=255)
                                                                                                                        atacgctgGTTTCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGc                             >  2:181494/5‑150 (MQ=255)
                                                                                                                          tggctgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAg                           >  2:60006/1‑150 (MQ=255)
                                                                                                                           ggctgATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGa                          <  1:534866/150‑1 (MQ=255)
                                                                                                                              tggtttCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGa                                                                                                   >  2:340965/4‑74 (MQ=37)
                                                                                                                              tggtttCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGa                                                                                                   <  1:340965/71‑1 (MQ=37)
                                                                                                                                gtttCGCTGCGCGGCGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATgg                     >  2:263389/2‑150 (MQ=255)
                                                                                                                                aTTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATgg                     <  2:454632/150‑1 (MQ=255)
                                                                                                                                aTTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATgg                     <  1:295346/150‑1 (MQ=255)
                                                                                                                                          gcAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGAAGAAGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGcc           >  1:315230/1‑150 (MQ=255)
                                                                                                                                              cGAAATCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCGCTc       <  1:181494/150‑1 (MQ=255)
                                                                                                                                                 aaCCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCGCTCGGt    >  2:170617/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGAAGGGAGAACAGGTCGCCACTGGATACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCGCTCGGTca  >  2:19298/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTCGCCACGTTCGGCCAGGGCCAGTTCGATGACCCGGCACACGGTATGGTCGTTGAGCAAGGTAATGATCTGGGTGATCTGCGTCGACGATGCGCCATGGGCCAGCATGCGTTCTACAAGCTGGCTGATTTCGCCGCGCAGCGAAACCAGCGTATCCAGGCGTGGCGCATGGCGTATGGTGCGCGCCAGGTGAACCAGATCGACCCGTTGCAGGGAGAACAGGTCGCGCTCGGAAACCACGCCGCACAGGCGTCCGTTTTCCACCAGGCAGACATGGGCAATATGCCGCTCGGTCA  >  NC_002947/431218‑431514

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: