Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,216,855 G→A 8.9% E1034K (GAA→AAA)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,216,8550GA 8.9% 132.7 / 3.7 45E1034K (GAA→AAA) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGT  >  CP009273/4216705‑4217059
                                                                                                                                                      |                                                                                                                                                                                                            
cTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTAAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTaa                                                                                                                                                 >  2:313174/1‑212 (MQ=255)
    aTAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCa                                                                                                                                             <  1:313179/212‑1 (MQ=255)
            atctccGAGAAAACGCTGAATCCGCTTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTTGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCgca                                                                      <  1:554990/271‑1 (MQ=255)
                             gaaggcgcgTGGTGTGGTGGGCCAGTTCCCGGAAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGAACGAAAAAACAGGATTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTacaa                                                       <  1:257922/268‑1 (MQ=255)
                                             gtggGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTAGCTGACTTCGTTGCGCCGAAGCTTTCTGGt                                                                                                                                                   >  2:457725/1‑165 (MQ=255)
                                              tggGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCAAAAAAACAGGCTTCGCTAactac                                                                                                                                                                                       >  2:316290/1‑128 (MQ=255)
                                                 gCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTAGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAg                                                                                                                                               <  1:457731/165‑1 (MQ=255)
                                                  ccTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCAAAAAAACAGGCTTCGCTAACTACTGtc                                                                                                                                                                                   <  1:316295/128‑1 (MQ=255)
                                                          cGGCAAACCGTGTGGGCGATGACATTAAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCTCATTTGCCGTGANTGGCGGGCNGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGt                        >  1:437014/1‑275 (MQ=255)
                                                           ggCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCAAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTgc                                                                                                                                                                   >  2:209622/1‑135 (MQ=255)
                                                               aaCCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCAAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCg                                                                                                                                                               <  1:209626/135‑1 (MQ=255)
                                                                    tgtgGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGtct                                                                                                                                                                                  >  2:206557/1‑111 (MQ=255)
                                                                        ggCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCt                                                                                                                                                                              <  1:206561/111‑1 (MQ=255)
                                                                        ggCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCt                                  >  1:374942/1‑251 (MQ=255)
                                                                            atgtcatgtaaaaCTACCCTGACGAAACGCCTAAACATGTGATCCCCGTTAGCCACAATCCGCGTATACATACCGAAAAAACTGTCTGCGCTGCCTAATGTCTCGTTGAATTCGTAGCGCCGAAGGTATCTGTTAAGGCAGATTACATCGGCGCATTTGCCGAGACGGGCTGGCATGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACTAAATCATGGTTAACGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCg      <  2:646880/266‑1 (MQ=255)
                                                                              gACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACa                                                                                                                                                                                                      >  2:378704/1‑81 (MQ=255)
                                                                              gACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGt                                                                                                                                                   >  2:505925/1‑132 (MQ=255)
                                                                                  ttGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCt                                                                                                                                                                                                  <  1:378709/81‑1 (MQ=255)
                                                                                  ttGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCTAAGCTTTCTGGTAAAg                                                                                                                                               <  1:505935/132‑1 (MQ=255)
                                                                                   tGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACt                                                                                                                  >  2:171429/1‑160 (MQ=255)
                                                                                       aTCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCg                                                                                                              <  1:171432/160‑1 (MQ=255)
                                                                                         cTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGt                                                                                                                                                   >  2:347879/1‑121 (MQ=255)
                                                                                             cGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAg                                                                                                                                               <  1:347884/121‑1 (MQ=255)
                                                                                                      aCGCGTACCCATGTGATCAACGTCAGCCACCATCCGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTTACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTCACTGGCGGGCTGGAAGAGGACGCACTGGCTGACGCCTTTGAAGCGNAGCACGATNATTACAACAAAATCATGGTGAAAGCTCTTGCCGACCGTTTAGCCGAAGCCTTTGCg      >  1:453542/1‑249 (MQ=255)
                                                                                                                 tgtgATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACt                                                                                                                  >  2:380003/1‑130 (MQ=255)
                                                                                                                  gtgATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCgg                                                                                                                                  >  2:276209/1‑113 (MQ=255)
                                                                                                                     aTCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCg                                                                                                              <  1:380008/130‑1 (MQ=255)
                                                                                                                      tCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCa                                                                                                                              <  1:276214/113‑1 (MQ=255)
                                                                                                                        aaCGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCtt                                                                                                                                                                                                 >  1:343225/1‑44 (MQ=255)
                                                                                                                        aaCGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGa                                                                                                                    >  1:41860/1‑121 (MQ=255)
                                                                                                                          cGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGt                                                                                                                                                   >  2:497245/1‑88 (MQ=255)
                                                                                                                          cGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGa                                                                                                       >  2:359170/1‑132 (MQ=255)
                                                                                                                            tCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCt                                                                                                                                                                                             <  2:343220/44‑1 (MQ=255)
                                                                                                                            tCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTgg                                                                                                                <  2:41860/121‑1 (MQ=255)
                                                                                                                              aGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTa                                                                                                                                                                                            >  1:635137/1‑43 (MQ=255)
                                                                                                                              aGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGATTTCGTTGCGCCGAAGCTTTCTGGTAAAg                                                                                                                                               <  1:497253/88‑1 (MQ=255)
                                                                                                                              aGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag                                                                                                   <  1:359175/132‑1 (MQ=255)
                                                                                                                                  accaACTGCGTCAACAGACCGAAAAAACAGGCTTAGCTAacca                                                                                                                                                                                        <  2:635126/43‑3 (MQ=255)
                                                                                                                                          cGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGa                                                                                                                    >  1:548302/1‑103 (MQ=255)
                                                                                                                                            tCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAgcg                                                                        >  1:81788/1‑145 (MQ=255)
                                                                                                                                              aaCAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTgg                                                                                                                <  2:548291/103‑1 (MQ=255)
                                                                                                                                              aaCAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGt  <  1:74467/213‑1 (MQ=255)
                                                                                                                                              aaCAGACCGAAAAAACAGGCTGCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTTTGGTAAAGCAGATTACATCGGCGCATTTGCCGTTACTTGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGa                <  1:547945/199‑1 (MQ=255)
                                                                                                                                                cAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCgcagc                                                                    <  2:81788/145‑1 (MQ=255)
                                                                                                                                                    ccGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGa                                                                                                                    >  2:216111/1‑93 (MQ=255)
                                                                                                                                                      |                                                                                                                                                                                                            
CTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAATCATGGTGAAAGCGCTTGCCGACCGTTTAGCCGAAGCCTTTGCGGAGT  >  CP009273/4216705‑4217059

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