Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,334,835 A→T 19.2% intergenic (+236/‑27) folK → / → panB 2‑amino‑4‑hydroxy‑6‑hydroxymethyldihydropteridine pyrophosphokinase/3‑methyl‑2‑oxobutanoatehydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,334,8350AT19.2% 96.9 / 10.1 51intergenic (+236/‑27)folK/panB2‑amino‑4‑hydroxy‑6‑hydroxymethyldihydropteridine pyrophosphokinase/3‑methyl‑2‑oxobutanoatehydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (17/24);  new base T (8/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 3.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCGCCAGCCAGG  >  ACB122/5334704‑5334967
                                                                                                                                   |                                                                                                                                    
ataGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAgc                                                                                                                                  >  7:249503/1‑136 (MQ=255)
  aGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAg                                                                                                                                   <  3:241232/133‑1 (MQ=255)
    gCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGACCGAACATCACGCGCGATGATGTGCCGCTCAGGAAGATGAATACACGCTTTGAATGCAGTCGTTTTTTAAAGCGcc                                                                                                                               >  6:73445/1‑135 (MQ=255)
       tCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACAGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGGTTTAGAgcggccgc                                                                                                                            >  1:194965/1‑130 (MQ=255)
         ccGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACAGCGccggcac                                                                                                                          >  1:147596/1‑130 (MQ=255)
         acGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGCCGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCAGAACg                                                                                                                         >  8:119771/2‑136 (MQ=255)
           gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCGGAACGa                                                                                                                        >  5:294333/1‑135 (MQ=255)
           gtgtGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATCCAGTCGTGTTTTAAAGCCCCTGAAAGa                                                                                                                        >  3:105824/1‑135 (MQ=255)
                  aCCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTc                                                                                                                <  7:31481/136‑1 (MQ=255)
                   ccAAAGTGCCGCATTCCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGGTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTCAGCGCCTGAACGAGGAATTTc                                                                                                                <  1:305780/135‑1 (MQ=255)
                   ccAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAAttt                                                                                                                 >  4:25140/1‑134 (MQ=255)
                   ccAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTcc                                                                                                               <  8:232070/136‑1 (MQ=255)
                   ccAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTc                                                                                                                <  3:145202/135‑1 (MQ=255)
                    cAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTcc                                                                                                               <  5:177559/135‑1 (MQ=255)
                        gTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACAt                                                                                                          <  7:54210/136‑1 (MQ=255)
                         tGCCGCATACCCGTATTTAGGCCCGGACGGCCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTACATCGCCTGAACGAGGAATTTCCAACaa                                                                                                          >  4:80767/1‑134 (MQ=255)
                         tGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTGTTTTCAGGGCCCGCACGAGGCATTTCCTTCAt                                                                                                          >  6:66001/1‑135 (MQ=255)
                             gCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTg                                                                                                     <  8:338555/136‑1 (MQ=255)
                              cATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTGTTTACAGCGCCTGAACGAGGTATTTACTAAATGCCt                                                                                                      >  2:255408/1‑134 (MQ=255)
                              cATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGCCGGCACGCGGAAATTACTTCATGccc                                                                                                      >  1:149363/1‑133 (MQ=255)
                                tACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGaa                                                                                                   <  5:307409/135‑1 (MQ=255)
                                   ccGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGt                                                                                                 <  3:25140/134‑1 (MQ=255)
                                        tttAGGCCCGGACGGACCTCTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGGCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCt                                                                                           <  8:342758/135‑1 (MQ=255)
                                                cggacggaCCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCAccc                                                                                   >  4:179119/1‑135 (MQ=255)
                                                cggacggaCCTGTGCCCGAACATAACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGAGTTGAATGAAGTCGTTTTTTTAATCGACGGAACGAGGAATTTCCTACATGGCTGAAGTAACCCTTACCAccc                                                                                   >  4:132011/1‑135 (MQ=255)
                                                   acggacCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCATCGCCGGGACGAGGGAATTGCTAACAGGCGGTAGTATCCCGTACCACCCGTa                                                                                >  1:370543/1‑135 (MQ=255)
                                                     ggacCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACCGCGCCGGCACGCGGAAGTTCCGAACTGCCTGAAGTTACCCGTACCAcccggaaa                                                                              >  4:206136/1‑130 (MQ=255)
                                                               ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTCAg                                                                                                                                   >  8:156986/1‑72 (MQ=255)
                                                               ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTTCAg                                                                                                                                   <  7:156986/72‑1 (MQ=255)
                                                               ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCGGAACGACGAAATTACTAAAAGCCTGAACTTACCCTTAACACCCATAACCGCCACAAgg                                                                    >  1:233878/1‑135 (MQ=255)
                                                               ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTAAGCGCCTGAACGAGGAATTTCCTACATGCCGGAAGTATACCTTAACCCCCATAACGGCCTCAAgg                                                                    >  3:372329/1‑135 (MQ=255)
                                                               ccGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGGTTTACAGCGCCCCAACGAGGTATTTCAAAAATGCCTGAAGTAACCCATAAAAACCCAAACCGCCACAAgg                                                                    >  1:340360/1‑135 (MQ=255)
                                                                  aaCATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTTAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCa                                                                 <  3:392342/135‑1 (MQ=255)
                                                                    cATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTCAGCGGCTGCACGCGGAATTTGCTACATGCCGGAAGTAACCCATACCACCCATTACGGCCGCAAGGCCCAg                                                               >  5:75072/1‑135 (MQ=255)
                                                                         cgcgcgATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGa                                                          <  6:242788/135‑1 (MQ=255)
                                                                               atgatgTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTGTTTTTAATCGCCGGAACGACGAATATCATAAATGCCGCAACTAACCCATACCACCCCAAACAGCCCCAACGCCCACGGTGaaaaaa                                                     >  2:156279/1‑134 (MQ=255)
                                                                                   tgtgCCTCTCCGGAAGATGACTACTCGCGTTGTATGGAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCAcc                                                <  4:105824/135‑1 (MQ=255)
                                                                                        cGCTCCGGAAGATGACTACACTCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCt                                           <  2:194965/135‑1 (MQ=255)
                                                                                        cGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCt                                           <  5:206619/135‑1 (MQ=255)
                                                                                               gAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGt                                                                                                 <  4:169079/74‑1 (MQ=255)
                                                                                               gAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGt                                                                                                 >  3:169079/1‑74 (MQ=255)
                                                                                                   aTGACTACACGCGTTGAATGCAGTCGTTTTTTAGAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTacg                                <  1:176138/135‑1 (MQ=255)
                                                                                                             gcgTTGAATGCAGTCGTTTTTTACCGCGCCAGCACGAGGAATTTACTACATTACTGAATTTACCCATAACACCCATAAACGCCGCAAGGCCAACGGTGAAAAAATCACCATGCTGACCAGCTTCGACGCGAcca                       >  8:43642/1‑133 (MQ=255)
                                                                                                             gcgTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGa                                                                                                                     <  7:223350/40‑1 (MQ=255)
                                                                                                             gcgTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGa                                                                                                                     >  8:223350/1‑40 (MQ=255)
                                                                                                                     tGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAgg              <  2:147596/135‑1 (MQ=255)
                                                                                                                          tCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGgccgcc         <  2:153163/135‑1 (MQ=255)
                                                                                                                               ttttACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCgccagcca    <  8:10213/135‑1 (MQ=255)
                                                                                                                                tttACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGTAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCgccagcca    <  4:75683/134‑1 (MQ=255)
                                                                                                                                tttACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAACTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCgccagcca    <  4:348339/134‑1 (MQ=255)
                                                                                                                                  tACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTATCGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCGCCAGCCAgg  <  3:206136/134‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
ATAGGCGTCCCGTGTGGCACCAAAGTGCCGCATACCCGTATTTAGGCCCGGACGGACCTGTGCCCGAACATCACGCGCGATGATGTGCCGCTCCGGAAGATGACTACACGCGTTGAATGCAGTCGTTTTTTACAGCGCCTGAACGAGGAATTTCCTACATGCCTGAAGTAACCCTTACCACCCTTAACGGCCTCAAGGCCAAGGGTGAAAAAATCACCATGCTGACCTGCTACGACGCGACCTTCGCCAAGGCCGCCAGCCAGG  >  ACB122/5334704‑5334967

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