Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,505,096 C→T 7.5% A309T (GCA→ACA)  yhfS ← FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,505,0960CT 7.5% 144.8 / 3.3 54A309T (GCA→ACA) yhfSFNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (2/2);  total (27/27)
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.52e-01

ATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCGGCACCGCACCACCGTTAAG  >  NC_000913/3504948‑3505243
                                                                                                                                                    |                                                                                                                                                   
aTACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGc                                                                                                                                                     <  2:329469/149‑1 (MQ=255)
 tcctctCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAGCGCACAATGTTCTGATGGTGGATTCGCCAGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCAGATTTCGACTCTGcc                                                                                                                                                    >  2:787921/3‑148 (MQ=255)
        aCGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAATCGGGTTAATGCGAATAGCACAATGTTCTGATTGTCGATTCGCATGGCGAAACGTTCCGGAAAGGCGATAAAATAGCGGAGGGATTTCATATTTCGTCTCTGCaccgaccg                                                                                                                                             >  2:943745/1‑149 (MQ=255)
          gCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCAGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCggg                                                                                                                                           >  1:920153/1‑149 (MQ=255)
            aaaaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGTACCCACCGGGTa                                                                                                                                         >  1:605338/1‑149 (MQ=255)
            aaaaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGTACCCACCGGGTa                                                                                                                                         >  1:554457/1‑149 (MQ=255)
              aaTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaag                                                                                                                                       <  1:808939/149‑1 (MQ=255)
                            tctcTTCACCGCTGCGATTCGGGTTAATTCGACTCGCACTAAGAGCTGATTGAGGATTCGCCTGTCGAAACGTTCCGGAAAGGCGATAAACGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCgcgc                                                                                                                         <  1:460326/149‑1 (MQ=255)
                              tctTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCtt                                                                                                                       >  2:528626/1‑149 (MQ=255)
                                ttCCCCGCTGCGATTCGTGTTAATGCGCATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGAAAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCtttt                                                                                                                     <  1:181753/149‑1 (MQ=255)
                                ttCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATATAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCGACCGGGTAAGGCAAGGCACCGCGAtttt                                                                                                                     <  2:258512/149‑1 (MQ=255)
                                ttCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCtttt                                                                                                                     >  2:542999/1‑149 (MQ=255)
                                    ccGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTgcgc                                                                                                                 <  1:242680/149‑1 (MQ=255)
                                    ccGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTgcgc                                                                                                                 <  1:854298/149‑1 (MQ=255)
                                          cGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGGCGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTc                                                                                                           <  1:16386/149‑1 (MQ=255)
                                          cGATTCGGGTTAATGCGAATCGCACAATGTACTGATTGTGGATTCGCCTGGCGAAACGTTCCGTAAAGGCGATAAAAGAGCGTCGGGATTTCATATTTCGACTCTGCACCCAGCGGTTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTc                                                                                                           <  1:422884/149‑1 (MQ=255)
                                                   ttAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTGCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGCCTCTGAACCCACCGGGTAAGCCAAGTCCACGCGCTTTTGCGCCTCTTCCAGCActct                                                                                                  >  1:603703/1‑149 (MQ=255)
                                                       tGCGATTCGCACAATGTTCTGAGTGTGGATTCGCCTGGCGAAACGTTCCGTAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCATCCACCGGGTAAGGAAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggcg                                                                                              <  2:276136/149‑1 (MQ=255)
                                                        gCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGTACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggcgg                                                                                             <  2:554457/149‑1 (MQ=255)
                                                        gCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGTACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTggcgg                                                                                             <  2:605338/149‑1 (MQ=255)
                                                               gCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCgg                                                                                      <  1:415254/149‑1 (MQ=255)
                                                                cacaATGTTCTGATTGTGTATTCGCCTGGCGACACGTTCCGGAAAGGCGATAAACGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGc                                                                                     <  2:233371/149‑1 (MQ=255)
                                                                cacaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGc                                                                                     <  2:958712/149‑1 (MQ=255)
                                                                cacaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGc                                                                                     >  1:452722/1‑149 (MQ=255)
                                                                   aaTGTTCTGATTGTGGATTCTCCAGGCGAAACGTTCCTGAAAGGAGATAAAATAGAGGCGGTATTTCATATTTCGAATCTGCACCCACCGGGTAAGGCAAGGCACCGCGATTTTTCGCCTCTTCCAGCACTCTGGCGGCAATCGGCtga                                                                                  <  2:562236/149‑1 (MQ=255)
                                                                   aaTGTTCTGATTGTGGATGCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCtga                                                                                  <  1:57157/149‑1 (MQ=255)
                                                                        tCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCAGTCTGGCGGCATTCGGCTGATGaaa                                                                             >  2:42083/1‑149 (MQ=255)
                                                                              tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTGATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCTCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAAGCGGCTGATGAGACTCGAc                                                                       >  2:207645/1‑149 (MQ=255)
                                                                                  gATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATc                                                                   >  1:917264/1‑149 (MQ=255)
                                                                                   aTTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTAATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCa                                                                  >  2:441321/1‑149 (MQ=255)
                                                                                      ggccgtgcgcAACGGTCCCGAAACGCGCTAAAAGCGCGGCGGGATTTCATATATAGACGCTGCCCCCACCCGGTAAGGCAAGGCACCGCGCGTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAaca                                                               <  1:672014/143‑1 (MQ=255)
                                                                                              gAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGAc                                                       <  2:920153/149‑1 (MQ=255)
                                                                                               aaaCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACt                                                      >  1:320631/1‑149 (MQ=255)
                                                                                               aaaCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCACTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACt                                                      >  1:884048/1‑149 (MQ=255)
                                                                                               aaaCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCAGCACCCACCGGGTAAGGCAAGGCACCGCGCGTTTGCGCCGCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACt                                                      <  1:104946/149‑1 (MQ=255)
                                                                                                   gTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTgcgc                                                  <  1:592416/149‑1 (MQ=255)
                                                                                                        ggAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAg                                             >  2:717847/1‑149 (MQ=255)
                                                                                                              gCGATAAAAGAGCGGCGTGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCa                                       <  2:376651/149‑1 (MQ=255)
                                                                                                                 aTAAAAGTGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTTGGGGCAATCGGCTGATTAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACcg                                    <  2:603703/149‑1 (MQ=255)
                                                                                                                   aaaaGAGCGGCGGGATTTCATATTTCGACTCTGCACCCCCAGGGTAAGGCAAGGCACCGCGCTTTAGCGCCTCTTCCCGCACTCTGGAGGCAATCGGCTGCTGCAACTCGACAATCAGCACCGTCGACTGCGCAATAGCAATCACcgcg                                  >  1:340096/1‑149 (MQ=255)
                                                                                                                     aaGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCt                                >  1:1013620/1‑149 (MQ=255)
                                                                                                                       gagCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCTGCACTCTGGCGGCAATCGGATGATGAAACTCGACAATCAACCCCTTCGACTGCGCATTAGCAATCACCGCGCttt                              >  1:580688/1‑149 (MQ=255)
                                                                                                                       gagCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACAGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACACGACAATCAACACCTTCGACTGCGCGTTAGCGATCACCGCGCttt                              >  1:328525/1‑149 (MQ=255)
                                                                                                                        agCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCtttt                             <  1:475946/149‑1 (MQ=255)
                                                                                                                          cggcggGATTTCATATATCGACTCTGCACCCACCGGGTAAGGCACGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCa                           <  2:949476/149‑1 (MQ=255)
                                                                                                                           ggcggGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCACTCACCGCGCTTTTCAc                          >  2:536106/1‑149 (MQ=255)
                                                                                                                                gATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCATCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCCACACCTTCGACTGCGCATTAGCATTCACCGCGCTTTTCACTTACg                     >  1:109197/1‑149 (MQ=255)
                                                                                                                                 aTTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCgg                    <  1:633041/149‑1 (MQ=255)
                                                                                                                                 aTTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCgg                    <  2:356390/149‑1 (MQ=255)
                                                                                                                                   ttCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACCCCTTCGACAGCGCATTAGCAATCACCGCGCTTTTCACTTCCGgca                  >  1:931664/1‑149 (MQ=255)
                                                                                                                                    tCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCGgcac                 <  2:452722/149‑1 (MQ=255)
                                                                                                                                    tCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCGGATGAAACGCGACAATCAACACCTGCGACGGCGCATTAGCCCTCCCCGCTCTTTTCACTTCCGgccc                 >  1:647960/1‑147 (MQ=255)
                                                                                                                                        atTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGATTTTCACTTCCGgcaccgca             >  1:1029484/1‑149 (MQ=255)
                                                                                                                                         tttCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCGgcaccgcac            >  1:805138/1‑149 (MQ=255)
                                                                                                                                             gACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGACACTCGACACTCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCGGCACCGCaccacc        >  2:298246/1‑149 (MQ=255)
                                                                                                                                                   gCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCGTCCAGCACTCTGGCGGCAATCGGCTGATGAACCTCGACAATCAACACCTTCGACTGCGCATTAGCAATCCCCGCGCTTTTCACATCCGGCCCAGCACCACCGTTCAg  >  1:1017368/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGACAATCAACACCTTCGACTGCGCATTAGCAATCACCGCGCTTTTCACTTCCGGCACCGCACCACCGTTAAG  >  NC_000913/3504948‑3505243

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