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,8890CA97.0% 112.1 / NA 33D388Y (GAT→TAT) folCtetrahydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (23/9);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTC  >  CP000730/1761754‑1762026
                                                                                                                                       |                                                                                                                                         
cttcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATatacta                                                                                                                                      >  1:243957/1‑141 (MQ=255)
cttcATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATatacta                                                                                                                                      >  1:382647/1‑141 (MQ=255)
    aTGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTg                                                                                                                                  >  1:247629/1‑141 (MQ=255)
      ggCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTa                                                                                                                                >  2:284004/1‑141 (MQ=255)
         tATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTAttt                                                                                                                             <  1:362715/141‑1 (MQ=255)
              tagttttagtttTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTtcatca                                                                                                                        >  1:187982/1‑141 (MQ=255)
                gttttagttttaTAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCAtt                                                                                                                      >  2:201398/1‑141 (MQ=255)
                          tataATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAAtacta                                                                                                            >  2:233630/1‑141 (MQ=255)
                          tataATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAAtacta                                                                                                            >  1:296619/1‑141 (MQ=255)
                            taATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATActacc                                                                                                          <  1:247090/141‑1 (MQ=255)
                            taATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATActacc                                                                                                          <  2:335468/141‑1 (MQ=255)
                                gCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACa                                                                                                      >  1:199456/1‑141 (MQ=255)
                                 cTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTCGTCATCTACTCCTTCTATTTCATCATTTTCAATACTACCTACaa                                                                                                     <  2:116697/141‑1 (MQ=255)
                                     aaaGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttc                                                                                                 >  2:273284/1‑141 (MQ=255)
                                     aaaGTCTAATTTTGATTTAACTGCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATttc                                                                                                 >  2:379512/1‑141 (MQ=255)
                                          cTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTcttctt                                                                                            >  1:148013/1‑141 (MQ=255)
                                                       aaCTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCttt                                                                               <  1:214958/141‑1 (MQ=255)
                                                                tATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATc                                                                      >  2:194352/1‑141 (MQ=255)
                                                                   gAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTCCAATTTCTTCTTTCCGTAACGCTTTCCGACAACCaaa                                                                   >  2:366399/1‑141 (MQ=255)
                                                                             cTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCa                                                         >  1:119245/1‑141 (MQ=255)
                                                                             cTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCa                                                         <  2:247629/141‑1 (MQ=255)
                                                                                  ggTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATaaaa                                                    <  2:296619/141‑1 (MQ=255)
                                                                                          cTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATgcg                                            >  2:209564/1‑141 (MQ=255)
                                                                                             aTGTATCACCTTGATAATTTTTTTTATATTCACCGTAGTCATATACTAATTGGATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCaa                                         <  1:174526/141‑1 (MQ=255)
                                                                                                  tCACCTTGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTaaa                                    >  1:354686/1‑141 (MQ=255)
                                                                                                    aCCTTGATAATTTTTTATTTTTTCTACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACTCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACt                                  >  2:250508/1‑141 (MQ=255)
                                                                                                        tGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTaaaa                              >  2:371167/1‑141 (MQ=255)
                                                                                                        tGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaaa                              >  2:371238/1‑141 (MQ=255)
                                                                                                        tGATAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTaaa                              >  2:348165/1‑141 (MQ=255)
                                                                                                           tAATTTTTTATAAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGt                           >  2:173693/1‑141 (MQ=255)
                                                                                                                   tattattttaaCGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAt                   >  2:353763/10‑141 (MQ=255)
                                                                                                                    ataAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAtc                  <  1:49912/141‑1 (MQ=255)
                                                                                                                    ataAATTCAACGTAGTCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCAtc                  >  2:135176/1‑141 (MQ=255)
                                                                                                                                   tCATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGtttt   <  1:201398/141‑1 (MQ=255)
                                                                                                                                    cATATACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTc  >  1:41873/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CTTCATGGCTATTTTAGTTTTAGTTTTATAATGCTTCAAAGTCTAATTTTGATTTAACTTCACTTATGAAATACAGACTACCGGTAATTACTAATGTATCACCTTGATAATTTTTTATAAATTCAACGTAGTCATCTACTAATTGTATTTCATCATTTTCAATACTACCTACAATTTCTTCTTTGCGTAACGCTTTCGGAAAATCAAATTCAGTTGCATAAAACGTATGCGCAATTAAACTTAAATGTTTGACCATCTCGTTAATCGGTTTTC  >  CP000730/1761754‑1762026

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