Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,089,887 T→G 11.1% K109T (AAA→ACA)  metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,089,8870TG11.1% 90.6 / 5.3 36K109T (AAA→ACA) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/12);  new base G (2/2);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GCGAATTCCCAGGCCAGTTTCTCGCGGCTTTTGTAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGTTCGCCGTC  >  NC_002947/3089743‑3090032
                                                                                                                                                |                                                                                                                                                  
gCGAATTCCCAGGCCAGTTTCTCGCGGCTTTTGTAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGttttt                                                                                                                                                 <  1:161063/148‑1 (MQ=255)
gCGAATTCCCAGGCCAGTTTCTCGCGGCTTTTGTAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAAAACCGGTGTTTg                                                                                                                                                >  5:229331/1‑149 (MQ=255)
     ttCCCAGGCCAGTTTCTCGCGGCTTTTGTAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAAAACCGGTGTTTGTCGaa                                                                                                                                           <  6:229331/149‑2 (MQ=255)
       cccAGGCCAGTTTCTCGCGGCTTTTGTAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGc                                                                                                                                          <  7:196141/148‑1 (MQ=255)
                                 tAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGc                                                                                                               <  6:180057/148‑1 (MQ=255)
                                    tGGTTGTCTTACAGTGTGTCGAACATGGTCAAGGGGCCGGGCAGCGCCGAGTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACGGCACCGACCACCGGCGGT‑ACGCTCg                                                                                                            >  2:195791/1‑148 (MQ=255)
                                            aTACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGc                                                                                                   >  8:28400/1‑149 (MQ=255)
                                                         aTCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCt                                                                                      >  7:188846/1‑149 (MQ=255)
                                                            atggtcatggGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGAc                                                                                   <  3:140571/149‑1 (MQ=255)
                                                                  atggGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGtttt                                                                             <  6:125744/149‑1 (MQ=255)
                                                                       gCCGGGCAGCGCCCACCTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGGGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTCGCGCTCACGGAGTCTGACGGTTTTCCgtt                                                                        >  8:57422/1‑149 (MQ=255)
                                                                              aGCGCCCACTTTATCGGGTGTGTTGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGTGTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGaa                                                                 <  8:25951/149‑1 (MQ=255)
                                                                                 gcCCACTTTATCGGCGGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCGCGGATTCTGACGGTTTTCCGTTGTTCGAAATc                                                              <  5:72104/149‑1 (MQ=255)
                                                                                  cccACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATc                                                              >  6:192076/1‑148 (MQ=255)
                                                                                  cccACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTAGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTATGACGGTTTTACGTTGTTCGAAATc                                                              >  7:90652/1‑148 (MQ=255)
                                                                                   ccACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCaa                                                            >  3:183550/1‑149 (MQ=255)
                                                                                          aTCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTGTCCGTTGTTCGAAATCAACGCCGCt                                                     >  4:227335/1‑149 (MQ=255)
                                                                                                 gtgtggtgtgCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAAcgccg                                                       >  4:284620/1‑140 (MQ=255)
                                                                                                 gtgtggtgtgCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAAcgccg                                                       <  3:284620/140‑1 (MQ=255)
                                                                                                     ggtgtgCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGa                                          >  1:138591/1‑149 (MQ=255)
                                                                                                     ggtgtgCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGa                                          >  4:288447/1‑149 (MQ=255)
                                                                                                           cTGGCTCAAGAATTTGGCATCTTCGAGAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAAc                                     <  8:225528/148‑1 (MQ=255)
                                                                                                           cTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCTCATCGTAGCGATCACGGATTCTGACGTTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAAc                                     <  1:195791/148‑1 (MQ=255)
                                                                                                           cTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAAc                                     <  7:28400/148‑1 (MQ=255)
                                                                                                              gcgcAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAg                                 >  6:22389/1‑149 (MQ=255)
                                                                                                                gcAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGt                                >  1:46501/1‑148 (MQ=255)
                                                                                                                  aaGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGa                              >  2:16397/1‑148 (MQ=255)
                                                                                                                   aGAATTTGGCATCTTCGACAAAGCCCGGTTTTTGTCGAGCGAAAGCACCCACCACCGACGGT‑ACGCTCGCATCGTAGCGAACAGGGATTCTGACGGTTTTCCGAGGTACGAAAAAAACGCCGCAGAGGTGCTCAATAAACGTAGTGAc                             >  7:49249/1‑148 (MQ=255)
                                                                                                                            cATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGTAACGCTCGCAACGAACCGAACCCGGATTCGAACGGTTTTCCGTTGTTCGAAAACAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTg                    >  1:84327/1‑149 (MQ=255)
                                                                                                                            cATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTgg                   >  5:205920/1‑149 (MQ=255)
                                                                                                                                ctCGACAAAAACCGGTGTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAAc                                                           <  4:97644/104‑1 (MQ=255)
                                                                                                                                ctCGACAAAAACCGGTGTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAAc                                                           >  3:97644/2‑105 (MQ=255)
                                                                                                                                      aaaCACCGGTTTTTGTCGAGCGACAGCACCGACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGtt         >  3:22618/1‑149 (MQ=255)
                                                                                                                                        acacCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGTTCg       >  1:49017/1‑149 (MQ=255)
                                                                                                                                           ccGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGTTCGCCg    >  5:244200/1‑149 (MQ=255)
                                                                                                                                            cGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGTTCGcctt   <  4:22618/149‑3 (MQ=255)
                                                                                                                                             ggTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGTTCGCCGTc  <  5:22389/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCGAATTCCCAGGCCAGTTTCTCGCGGCTTTTGTAGTGGTTGTCATACAGTGTGTCGATCATGGTCATGGGGCCGGGCAGCGCCCACTTTATCGGCTGTGTGGTGTGCTGGCGCAAGAATTTGGCATCTTCGACAAACACCGGTTTTTGTCGAGCGACAGCACCCACCACCGTCGGT‑ACGCTCGCATCGTAGCGATCACGGATTCTGACGGTTTTCCGTTGTTCGAAATCAACGCCGCTGAGGTGCTCGATAAACGTAGTGACAAAGTGCTGGCGGGTTTGTTCGCCGTC  >  NC_002947/3089743‑3090032

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: ‑