Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,761,889 C→A D388Y (GAT→TAT)  folC ← tetrahydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,761,8890CA100.0% 79.8 / NA 25D388Y (GAT→TAT) folCtetrahydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/6);  total (19/6)

TAGATTCTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTCCGT  >  CP000730/1761748‑1762029
                                                                                                                                             |                                                                                                                                            
tAGATTCTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTaa                                                                                                                                        >  1:481/1‑148 (MQ=255)
tAGATTCTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTaa                                                                                                                                        >  1:244809/1‑148 (MQ=255)
    ttcttcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGt                                                                                                                                    >  1:34767/1‑148 (MQ=255)
        tcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTAttt                                                                                                                                >  2:163721/1‑148 (MQ=255)
            ggCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTtcatc                                                                                                                            >  1:60596/1‑148 (MQ=255)
                      gttttagttttaTAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTTAAtt                                                                                                                  >  1:64604/1‑147 (MQ=255)
                               ttaTAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACa                                                                                                         <  2:219618/148‑1 (MQ=255)
                                     tGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttct                                                                                                   >  2:222770/1‑148 (MQ=255)
                                     tGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttct                                                                                                   >  2:187379/1‑148 (MQ=255)
                                     tGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttct                                                                                                   >  2:222754/1‑148 (MQ=255)
                                      gCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttctt                                                                                                  <  2:259608/148‑1 (MQ=255)
                                                   aTTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGc                                                                                     >  2:299955/1‑148 (MQ=255)
                                                     tttGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCtt                                                                                   <  2:15092/148‑1 (MQ=255)
                                                      ttGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCttt                                                                                  >  2:135757/1‑148 (MQ=255)
                                                         aTTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCgg                                                                               >  1:285182/1‑148 (MQ=255)
                                                               cTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATc                                                                         >  1:120802/1‑148 (MQ=255)
                                                                          aaaTACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTGTTCTTTGCGTAACGCTTTTGGACAAGCAACTTCAGTTg                                                              >  2:177104/1‑148 (MQ=255)
                                                                               cAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATaa                                                         <  2:241587/148‑1 (MQ=255)
                                                                                aGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATaaa                                                        >  1:120538/1‑148 (MQ=255)
                                                                                          tAATTACTAATGTATCACCTTGATAATTTTTTATATATTTAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATgcgc                                              >  2:163818/1‑148 (MQ=255)
                                                                                                         cACCTTGATAATTTTTTTTTTTTTTTTCGTAGTCATATACTAATTGTATTTTATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAATATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATg                               >  2:151380/1‑148 (MQ=255)
                                                                                                         cACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATg                               >  2:209046/1‑148 (MQ=255)
                                                                                                                            aaaTTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATc            >  1:132150/1‑148 (MQ=255)
                                                                                                                                    cGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTcc    <  1:300682/148‑1 (MQ=255)
                                                                                                                                      tAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAGACTTAAATGTTTGACCATCTCGTTAATCGGTTTTCCGt  <  2:61701/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
TAGATTCTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTCCGT  >  CP000730/1761748‑1762029

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