Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,042,793 T→C 9.1% E47E (GAA→GAG)  folE_1 ← GTP cyclohydrolase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,042,7930TC 9.1% 132.0 / 4.6 44E47E (GAA→GAG) folE_1GTP cyclohydrolase 1
Reads supporting (aligned to +/- strand):  ref base T (15/25);  new base C (2/2);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CCTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTAGTTCTGTT  >  ACB143/2042660‑2042923
                                                                                                                                     |                                                                                                                                  
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGCGCTGAACAAAGCGGTGTTGGTGActtc                                                                                                                                  >  3:4518/1‑136 (MQ=255)
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtc                                                                                                                                  >  1:158393/1‑136 (MQ=255)
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtc                                                                                                                                  >  8:148069/1‑136 (MQ=255)
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtc                                                                                                                                  >  7:270736/1‑136 (MQ=255)
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtc                                                                                                                                  >  4:154212/1‑136 (MQ=255)
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtc                                                                                                                                  >  2:64672/1‑136 (MQ=255)
ccTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtc                                                                                                                                  >  3:180883/1‑136 (MQ=255)
ccTTAGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCACGGACCAGCACCATTTCGCTGTTGTCAGAGCTTAACAAAGCGTTGTTGGTGAgttc                                                                                                                                  >  3:199344/1‑136 (MQ=255)
 cTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGActtct                                                                                                                                 <  2:81657/136‑1 (MQ=255)
      gCAAGTAGGCAACGTGTGCCTTGCGAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTTGTGACTTCTTCCAg                                                                                                                            <  5:109346/136‑1 (MQ=255)
      gCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAg                                                                                                                            <  5:63337/136‑1 (MQ=255)
       cAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAgt                                                                                                                           >  4:52098/1‑136 (MQ=255)
                       gCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACcgc                                                                                                           >  7:265944/1‑136 (MQ=255)
                        ccTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACcgcg                                                                                                          >  5:176639/1‑136 (MQ=255)
                                               tggtgTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCttt                                                                                    <  8:265944/135‑1 (MQ=255)
                                               tggtgTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTg                                                                                   <  7:217157/136‑1 (MQ=255)
                                                 gtgtTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGc                                                                                  <  7:175816/135‑1 (MQ=255)
                                                 gtgtTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGc                                                                                  <  2:297387/135‑1 (MQ=255)
                                                       gCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGc                                                                           >  5:151278/1‑136 (MQ=255)
                                                                     gCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAg                                                             <  5:164855/136‑1 (MQ=255)
                                                                        cGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTGATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAg                                                          <  1:76063/136‑1 (MQ=255)
                                                                        cGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAg                                                          >  3:243132/1‑136 (MQ=255)
                                                                         gATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAgg                                                         >  5:239646/1‑136 (MQ=255)
                                                                         gATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAgg                                                         >  7:155204/1‑136 (MQ=255)
                                                                         gATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAgg                                                         >  3:52754/1‑136 (MQ=255)
                                                                         gATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAg                                                          <  3:276295/135‑1 (MQ=255)
                                                                          aTGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGc                                                        <  3:173095/136‑1 (MQ=255)
                                                                             tCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGccc                                                      <  4:111724/135‑1 (MQ=255)
                                                                             tCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCt                                                     <  2:2281/136‑1 (MQ=255)
                                                                                  gACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACgg                                                <  3:259850/136‑1 (MQ=255)
                                                                                           cATTTCGCTGTTTTCAGAGATGAACAAAGAGTTGTTGGTGACTTCTTCAAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCCTCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTcc                                       >  7:137263/1‑136 (MQ=255)
                                                                                                cGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCACGCTTGGGCGTGTCGAGCAGGACCGCACGGGAGACGTCCTCGc                                   >  6:251173/1‑135 (MQ=255)
                                                                                                        tCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAg                                                                                                                            <  2:9442/38‑1 (MQ=255)
                                                                                                        tCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAg                                                                                                                            >  1:9442/1‑38 (MQ=255)
                                                                                                        tCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGc                          >  1:219434/1‑136 (MQ=255)
                                                                                                          agagCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCt                         >  4:207893/1‑135 (MQ=255)
                                                                                                            agCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAg                      <  4:215849/136‑1 (MQ=255)
                                                                                                              cTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGa                    <  6:288109/136‑1 (MQ=255)
                                                                                                              cTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGa                    <  8:143246/136‑1 (MQ=255)
                                                                                                              cTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGa                    <  8:23419/136‑1 (MQ=255)
                                                                                                                gAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCa                                                                                            >  1:157385/1‑62 (MQ=255)
                                                                                                                gAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCa                                                                                            <  2:157385/62‑1 (MQ=255)
                                                                                                                   cAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCg               <  3:2147/136‑1 (MQ=255)
                                                                                                                   cAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCg               <  1:196095/136‑1 (MQ=255)
                                                                                                                   cAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAACCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCg               <  7:47133/136‑1 (MQ=255)
                                                                                                                    aaaGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCgg              >  5:275184/1‑136 (MQ=255)
                                                                                                                    aaaGCGTTGTTGGTGACCTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCgg              >  3:131151/1‑136 (MQ=255)
                                                                                                                    aaaGCGTTGTTGGTGACCTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCgg              >  3:131153/1‑136 (MQ=255)
                                                                                                                        cGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTa          <  2:269124/136‑1 (MQ=255)
                                                                                                                          ttgttgGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTAGt        <  6:104283/136‑1 (MQ=255)
                                                                                                                                gTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTAGTTCTGtt  <  6:275184/136‑1 (MQ=255)
                                                                                                                                gTGACCTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTAGTTCTGtt  <  4:131153/136‑1 (MQ=255)
                                                                                                                                gTGACCTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTAGTTCTGtt  <  4:131151/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
CCTTCGGCAAGTAGGCAACGTGTGCCTTGCCAATGAACGGCAGCATGTGGTGTTCGCACATCGAATACAGCTCGATGTCCCGGACCAGCACCATTTCGCTGTTGTCAGAGCTGAACAAAGCGTTGTTGGTGACTTCTTCCAGTGTTTGCTCATAACCGCGGCAAAGGTACTTCATCGCCTTTGCAGCCCGCTTGGGCGTGTCGAGCAGGCCCTCACGGGAGACGTCCTCGCCAATCTGGCTGAGGATCTCGGTGTAGTTCTGTT  >  ACB143/2042660‑2042923

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