Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,813,442 T→C 16.2% G293G (GGT→GGC)  frmA → glutathione‑dependent formaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,813,4420TC16.2% 66.8 / 4.5 31G293G (GGT→GGC) frmAglutathione‑dependent formaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (13/12);  new base C (3/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 = 1.00e+00

AGGAAGTGATCGTCGACCTCACCGATGGTGGCGTGGACTTCTCCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGACACCT  >  NC_002947/1813306‑1813587
                                                                                                                                        |                                                                                                                                                 
aGGAAGTGATCGTCGACCTCACCGATGGTGGCGTGGACTTCTCCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGaaa                                                                                                                                      >  1:601518/1‑150 (MQ=255)
         tcgtcgACCTCACCGATGGTGGCGTGGACTTCTCCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCAcgc                                                                                                                             <  2:582784/150‑1 (MQ=255)
                    aCCGATGGTGGCGTGGACTTCTCCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAg                                                                                                                  <  2:588814/150‑1 (MQ=255)
                         cgggggCGTGTATTTCTCCTTCGATTGCATCGGCAACGTTCAACTGATGCGCGCCGACCTGGGATGCCACGACAAGCGCTGCGCACAGTCCCTCCCCTTCCGTGTGGCCGTCGTCGGCCATGAAATGTCCACTCGTCCATTCCAGCTGGt                                                                                                             <  2:291396/146‑1 (MQ=255)
                                gTGGACTTCTCCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTACCGCAGGTTCCCAGt                                                                                                      >  2:44562/1‑150 (MQ=255)
                                        ctcCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGACGGCCAGGAAATTTCGACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTg                                                                                              >  1:100505/1‑150 (MQ=255)
                                         tcCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCGGGTGCCGGCCAGAAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGGTTgg                                                                                             >  1:683535/1‑150 (MQ=255)
                                           cTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTAcc                                                                                                         >  1:208304/1‑136 (MQ=255)
                                           cTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTAcc                                                                                                         <  2:208304/136‑1 (MQ=255)
                                                     aTCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCt                                                                                                                <  1:598196/119‑1 (MQ=255)
                                                     aTCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCt                                                                                                                >  2:598196/1‑119 (MQ=255)
                                                           aaCGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGc                                                                                            >  1:35757/1‑133 (MQ=255)
                                                           aaCGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGc                                                                                            <  2:35757/133‑1 (MQ=255)
                                                            aCGTTCCACCGCTGCGCTCTGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGTTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGttcggcgttcggcg                                                                          <  2:100505/150‑1 (MQ=255)
                                                                            cgccgcCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGAAGGGGACGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAg                                                          >  2:308068/1‑150 (MQ=255)
                                                                               cgccCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGTGTGGCCGCAGCGa                                                       <  2:630507/150‑1 (MQ=255)
                                                                                     ggAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGcc                                                 >  2:530920/1‑150 (MQ=255)
                                                                                            tgcCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTAc                                          <  1:44562/150‑1 (MQ=255)
                                                                                               cacaAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTg                                       <  2:245974/150‑1 (MQ=255)
                                                                                                 caAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGa                                     >  1:203343/1‑150 (MQ=255)
                                                                                                                 tCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGc                     >  2:376298/1‑150 (MQ=255)
                                                                                                                   ggTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCga                   >  1:122323/1‑150 (MQ=255)
                                                                                                                   ggTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCga                   >  2:321635/1‑150 (MQ=255)
                                                                                                                       atcatcGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTg                                       <  2:156126/126‑1 (MQ=255)
                                                                                                                       atcatcGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTg                                       >  1:156126/1‑126 (MQ=255)
                                                                                                                            cGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCt          <  1:321635/150‑1 (MQ=255)
                                                                                                                               cctgGCCGGCACCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGa       >  2:698996/3‑150 (MQ=255)
                                                                                                                               cctgGCCGGCACCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGa       >  2:763359/3‑150 (MQ=255)
                                                                                                                                ctgGCCGGCACCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGac      >  1:796435/2‑150 (MQ=255)
                                                                                                                                ctgGCCGGCACCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGac      <  2:796435/149‑1 (MQ=255)
                                                                                                                                    ccGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGACACCt  <  2:199846/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
AGGAAGTGATCGTCGACCTCACCGATGGTGGCGTGGACTTCTCCTTCGAGTGCATCGGCAACGTTCAACTGATGCGCGCCGCCCTGGAATGCTGCCACAAGGGCTGGGGAGAGTCGGTCATCATCGGTGTGGCCGGTGCCGGCCAGGAAATTTCCACGCGTCCATTCCAGCTGGTTACCGGTCGGGTTTGGCGTGGTTCGGCGTTCGGCGGCGTGCGTGGCCGCAGCGAGCTGCCGAGCTACGTGGAAATGTCCGAGAAAGGCGAGATCCCGCTGGACACCT  >  NC_002947/1813306‑1813587

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