Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 766,621 T→G 8.8% V38G (GTG→GGG)  metB → cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947766,6210TG 8.8% 121.3 / 2.5 44V38G (GTG→GGG) metBcystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base T (23/17);  new base G (2/2);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAGCAAGCCGAGTC  >  NC_002947/766473‑766768
                                                                                                                                                     |                                                                                                                                                    
aTCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGt                                                                                                                                                      >  5:229716/1‑149 (MQ=255)
    gtgtACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGt                                                                                                                                                   >  4:107161/1‑148 (MQ=255)
      gtACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCTGACCACGATCGGGGTctg                                                                                                                                                >  6:180195/1‑147 (MQ=255)
        acacAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAgcg                                                                                                                                              >  6:275872/1‑149 (MQ=255)
         cacaACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCTGACCACGATCGGGGTCTgcgt                                                                                                                                             <  5:180195/149‑2 (MQ=255)
         cacaACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAgcgc                                                                                                                                             <  5:107231/149‑1 (MQ=255)
          acaACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGtgcc                                                                                                                                            <  1:38839/149‑1 (MQ=255)
            aaCTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACACGTACAACGCCACGCAGACGCCGATCGTGGTCAGcgccgc                                                                                                                                          >  4:65603/1‑149 (MQ=255)
                 aGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTAt                                                                                                                                      <  6:276737/148‑1 (MQ=255)
                  gCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATgg                                                                                                                                    <  6:176393/149‑1 (MQ=255)
                    tGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGc                                                                                                                                   >  2:262676/1‑148 (MQ=255)
                      ggCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCt                                                                                                                                  <  8:217840/147‑1 (MQ=255)
                         gggAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGTGAGCAGGTACAGCAAACCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACaa                                                                                                                              <  1:85188/149‑1 (MQ=255)
                                    atGAACAACAACGAAAAATTAGAGGGAAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGc                                                                                                                 >  2:6887/1‑149 (MQ=255)
                                      gAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCt                                                                                                               >  5:13391/1‑149 (MQ=255)
                                                 aaaaaTTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTa                                                                                                           <  7:253728/142‑1 (MQ=255)
                                                 aaaaaTTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTa                                                                                                           >  8:253728/1‑142 (MQ=255)
                                                      ttCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGgcg                                                                                                 >  2:315359/1‑147 (MQ=255)
                                                                   ctgcGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGTGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTCGCGCAGGGCAAGCAGcc                                                                                   <  1:6887/148‑1 (MQ=255)
                                                                        gCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAggg                                                                                            >  1:241774/1‑134 (MQ=255)
                                                                        gCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAggg                                                                                            <  2:241774/134‑1 (MQ=255)
                                                                         ccGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                            <  1:262676/149‑1 (MQ=255)
                                                                                 gggACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAg                                                                    >  2:262832/1‑149 (MQ=255)
                                                                                  ggACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGGCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTGCATCTACAGc                                                                   >  4:284819/1‑149 (MQ=255)
                                                                                   gACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGcc                                                                  <  6:229716/149‑1 (MQ=255)
                                                                                      gcgGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑a                                                                                                                        >  6:298573/1‑92 (MQ=255)
                                                                                      gcgGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑a                                                                                                                        <  5:298573/92‑1 (MQ=255)
                                                                                        gggTCGTGTGGAGCGGGGAGCAGGTACAGC‑TCCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                              >  1:28121/1‑148 (MQ=255)
                                                                                        gggTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                              >  4:113307/1‑148 (MQ=255)
                                                                                              tgtgGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAAc                                                       <  3:107161/149‑1 (MQ=255)
                                                                                              tgtgGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAAc                                                       >  1:62834/1‑149 (MQ=255)
                                                                                               gtgGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                            >  1:131255/1‑127 (MQ=255)
                                                                                                  gAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCa                                                   >  6:2557/1‑149 (MQ=255)
                                                                                                  gAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCa                                                   >  7:157016/1‑149 (MQ=255)
                                                                                                    gCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCGCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACg                                                 <  5:275872/149‑1 (MQ=255)
                                                                                                     cGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACgg                                                >  6:124081/1‑149 (MQ=255)
                                                                                                                     c‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTTCGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCCCGGTCGCCACGCTGGAAg                                 >  6:67865/1‑148 (MQ=255)
                                                                                                                           gTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCCCGGTCGCCACGCTGGAAGCAAAAc                           >  2:200274/1‑149 (MQ=255)
                                                                                                                                   cggcGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACGAACGATGCATGGTACGACGTGGCGCAGGCCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCACCCCCACGGACGCCa                                          >  1:71499/4‑127 (MQ=255)
                                                                                                                                   cggcGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCa                                          <  2:71499/123‑1 (MQ=255)
                                                                                                                                     aCGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCt                 >  2:108034/1‑149 (MQ=255)
                                                                                                                                        cAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGa              <  1:315359/149‑1 (MQ=255)
                                                                                                                                             cccGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAGCAAg         >  5:288606/1‑149 (MQ=255)
                                                                                                                                               cGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAGCAAGcc       <  1:78319/149‑1 (MQ=255)
                                                                                                                                                    gTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAGCAAGCCGAGTc  <  6:13391/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAGCAAGCCGAGTC  >  NC_002947/766473‑766768

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