Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,160,127 A→G P340P (CCT→CCC)  ygiQ ← radical SAM superfamily protein YgiQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,160,1270AG95.8% 74.2 / ‑4.6 24P340P (CCT→CCC) ygiQradical SAM superfamily protein YgiQ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (10/13);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACT  >  NC_000913/3159990‑3160250
                                                                                                                                         |                                                                                                                           
cagcCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCggg                                                                                                                           <  2:21251/140‑1 (MQ=255)
  gcCACGCATAATGTTGACCGAAAAACGGAGCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCg                                                                                                                         <  1:252337/140‑1 (MQ=255)
    cACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGcgggcggg                                                                                                                       <  2:394469/140‑1 (MQ=255)
    cACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGAATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGcgggcggg                                                                                                                       <  2:566887/140‑1 (MQ=255)
      cGCATAATGTTGACCTAAAAACTGAGCATATCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGtt                                                                                                                     <  1:343983/140‑1 (MQ=255)
          tAATGTTGACCGAAACACGGATTATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATc                                                                                                                 <  1:331572/140‑1 (MQ=255)
           aaTGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGCATAGCGGGCGGGATGATcc                                                                                                                >  1:141710/1‑140 (MQ=255)
                  aCCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATa                                                                                                         >  1:83195/1‑140 (MQ=255)
                       aaaaCGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGt                                                                                                    <  2:657728/140‑1 (MQ=255)
                                               aaTACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATgcgc                                                                            >  1:556090/1‑140 (MQ=255)
                                                   cGGGCATTGCCATAGTCCGGATATGGCACTCGCTTGTCTGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCGTCAATGCGCGGGc                                                                        <  1:351013/140‑1 (MQ=255)
                                                      gCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGACGCCGTGTATTTGAATCAATGCGCGGGcaca                                                                     >  2:211623/1‑140 (MQ=255)
                                                         ttGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGcc                                                                  <  2:457750/140‑1 (MQ=255)
                                                         ttGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGCAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAgac                                                                  >  2:682615/1‑138 (MQ=255)
                                                         ttGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGC‑AAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCt                                                                 >  1:93220/1‑140 (MQ=255)
                                                          tGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCt                                                                 >  1:629192/1‑140 (MQ=255)
                                                           gCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTg                                                                >  2:712773/1‑140 (MQ=255)
                                                                 ggCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttgg                                                          <  2:551519/140‑1 (MQ=255)
                                                                      gATGTGGCACGCGCTCGTATGGTAGCGCACAAACGCTATCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCg                                                     <  1:211623/140‑1 (MQ=255)
                                                                                      gTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                     <  2:629192/140‑1 (MQ=255)
                                                                                           gCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                <  2:646865/140‑1 (MQ=255)
                                                                                             aGCGC‑AAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                             <  2:93220/140‑1 (MQ=255)
                                                                                               cgcAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                            >  1:156864/1‑140 (MQ=255)
                                                                                                gcAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                           >  1:181768/1‑140 (MQ=255)
                                                                                                   aaaaCGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                        <  2:156864/140‑1 (MQ=255)
                                                                                                                  tctTCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGcc         <  2:556090/140‑1 (MQ=255)
                                                                                                                     tCGGTAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCTCTTTATCGCCCtt      <  2:518553/140‑1 (MQ=255)
                                                                                                                         tAGAAAGCGGAATAGCGGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATTCGCTAAGCATGGGCGTCCAGCACTTTATCGCCCTTCACt  >  2:289906/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                           
CAGCCACGCATAATGTTGACCGAAAAACGGATCATTTCGTAAGCCGGAATACGGGCATTGCCATAGGCCGGATGTGGCACGCGCTTGTATGGCAGCGCAAAAACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACT  >  NC_000913/3159990‑3160250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑