Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 387,351 C→A 12.9% A163E (GCG→GAG)  glyA‑I → serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947387,3510CA12.9% 79.7 / 2.6 31A163E (GCG→GAG) glyA‑Iserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/12);  new base A (2/2);  total (17/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.62e-01

ACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGCCCACGTTGCCGGCCTGGTTGCCGCTGGC  >  NC_002947/387206‑387499
                                                                                                                                                 |                                                                                                                                                    
acCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTc                                                                                                                                                  >  1:442879/1‑150 (MQ=255)
  cATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGGCGa                                                                                                                                                >  1:370971/1‑150 (MQ=255)
   aTCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGCCAGGCGCTCGac                                                                                                                                               <  1:644855/150‑2 (MQ=255)
      cTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGcac                                                                                                                                            >  2:213964/1‑150 (MQ=255)
         ggCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAg                                                                                                                                         <  1:213964/150‑1 (MQ=255)
            aTGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGAGGTCGAg                                                                                                                                               >  1:784901/1‑141 (MQ=255)
            aTGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGAGGTCGAg                                                                                                                                               <  2:784901/141‑1 (MQ=255)
                     gCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATc                                                                                                                             <  1:82048/150‑1 (MQ=255)
                               gTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGtt                                                                                                                   <  1:404643/150‑1 (MQ=255)
                                      gACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGc                                                                                                            <  1:295029/150‑1 (MQ=255)
                                      gACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGAGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGc                                                                                                            >  2:553867/1‑150 (MQ=255)
                                         ccACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTa                                                                                                         >  1:355232/1‑150 (MQ=255)
                                         ccACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTa                                                                                                         >  1:227345/1‑150 (MQ=255)
                                                  cAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGCACGGCATCGAGACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGAc                                                                                                >  1:217351/1‑150 (MQ=255)
                                                           gtcCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACtt                                                                                       >  2:621602/1‑150 (MQ=255)
                                                              cTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGc                                                                                                            >  2:685048/1‑126 (MQ=255)
                                                              cTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGc                                                                                                            <  1:685048/126‑1 (MQ=255)
                                                                        cTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCcgc                                                                          <  1:621602/150‑1 (MQ=255)
                                                                         tGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCcgcg                                                                         >  2:724056/1‑150 (MQ=255)
                                                                              aaCGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGCCAGGCGCTCGac                                                                                                                                               <  1:685020/75‑2 (MQ=37)
                                                                              aaCGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGCCAGGCGCTCGac                                                                                                                                               >  2:685020/1‑74 (MQ=37)
                                                                                         gTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGt                                                         >  2:437592/1‑150 (MQ=255)
                                                                                          tACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTg                                                        <  2:227345/150‑1 (MQ=255)
                                                                                                       ccAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGAGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGt                                           <  1:553867/150‑1 (MQ=255)
                                                                                                        cAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGtt                                          >  2:635720/1‑150 (MQ=255)
                                                                                                                     aTCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGcc                             >  2:311338/1‑150 (MQ=255)
                                                                                                                             cgacgaAGTCGAGCGCCTGGCGGGCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGCCCACGTTGc                     >  2:228354/1‑150 (MQ=255)
                                                                                                                              gacgaAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGCCCACGTTGcc                    >  1:309275/1‑150 (MQ=255)
                                                                                                                                       gcgcgcCTGGCGGGCGAGCACAAGCCGAAAATGAGCGGTGCCCGTAGCTCGGCCTACTGAAAGACCCTCGACCTCCCCCGTTTCCGCGCCATAGCCCACAAGGTGGGCGCACTGCTGGTCGTCGACATGGCCCACGTTGCCGGCCTGGtt           <  2:154180/148‑1 (MQ=21)
                                                                                                                                          cgcCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGCCTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGCCCACGTTGCCGGCCTGGTTGct        <  1:437592/150‑2 (MQ=255)
                                                                                                                                                gCGGTCGAGAACAAGACGACGATGACCGTGGCCTGTTTCTCGGCTTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGCCCACGTTGCCGGCCTGGTTGCCGCTGGc  <  2:370971/150‑1 (MQ=25)
                                                                                                                                                 |                                                                                                                                                    
ACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTTCGTCGACATGGCCCACGTTGCCGGCCTGGTTGCCGCTGGC  >  NC_002947/387206‑387499

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