Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,600,978 C→T intergenic (+30/‑223) bcp → / → hyfA thiol peroxidase/hydrogenase 4 component A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,600,9780CT90.0% 68.9 / ‑1.1 30intergenic (+30/‑223)bcp/hyfAthiol peroxidase/hydrogenase 4 component A
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (14/13);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

GCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGCTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCC  >  NC_000913/2600732‑2601221
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGc                                                                                                                                                                                                                                                 <  1:104119/251‑1 (MQ=255)
            tACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATAcc                                                                                                                                                                                                                                     >  1:126268/1‑251 (MQ=255)
            tACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATAcc                                                                                                                                                                                                                                     >  1:362636/1‑251 (MQ=255)
               gCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCACGCATACCTc                                                                                                                                                                                                                                   >  2:146243/1‑250 (MQ=255)
                cTCCTGTCTGATGAGGACCACCAGGAGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGCCGATTTCAAAACCAGCAATCACCCCGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCAATTACGTGCTGGCTGCGTTTGCGGCCAGCATACCTCAc                                                                                                                                                                                                                                 >  2:15089/1‑251 (MQ=255)
                            taGGACAACCAGGTGTGAGAACAATTAGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGCTTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGACTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCCATTCTCGtactcc                                                                                                                                                                                                                     >  2:341816/2‑246 (MQ=255)
                               gACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTTAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCTCTTCTCGTGATCAAGa                                                                                                                                                                                                                  <  1:70205/251‑1 (MQ=255)
                                    ccaGGTGTGCGAACAATTCGGCGTCAGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAACTTCCTGATTGACGCTGATGGCAAAATCGCACATGACTTTGACGCTTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCAGGCTGCGTTTGCGACCAGCCTACCTCACTTCTCTTGATCAAgtgcac                                                                                                                                                                                                              >  2:106900/1‑245 (MQ=255)
                                                      cGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGc                                                                                                                                                                                            <  2:105043/250‑1 (MQ=255)
                                                                         gccTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCAt                                                                                                                                                                        <  1:15089/250‑1 (MQ=255)
                                                                              cATGGGCAAAACCTACGATGACATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCAAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGAAGATGTGCTGAACAGGCTGAAAGAACACGCCGGCTTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGAGGCCCACAGACCTCCATTCTCGTGATCAAGAACACATTATCGCTTGCCCCTGCGCCACGGGTGTCGGATCCAtttttt                                                                                                                                                                   >  2:143280/1‑251 (MQ=255)
                                                                                    cAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGaa                                                                                                                                                             <  1:233675/251‑1 (MQ=255)
                                                                                      aaaCCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACg                                                                                                                                                           >  1:9245/1‑251 (MQ=255)
                                                                                       aacctaagtaggcattcaccccagtaGCTTCCTGATTGACGCTGATGGCAAAAACGAAAATTGCTTTGACGATTTCAAAACCAGCAATCCCCACGACGTTGTGCTTAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCATCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAGCGt                                                                                                                                                          <  2:126268/226‑1 (MQ=255)
                                                                                          cTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCCGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTa                                                                                                                                                        >  1:339626/1‑250 (MQ=255)
                                                                                            aCGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTaa                                                                                                                                                       <  1:106131/249‑1 (MQ=255)
                                                                                               aTGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGAcc                                                                                                                                                  >  1:213800/1‑250 (MQ=255)
                                                                                                        atcccATCAGCATCCGGATTGACGCTGATGGCAAAACCAACCATGTCCTTGACGAATTCAAAACCAGCAATCCCCTCGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCGCGCTTTCCCCTGCTACATGTGTGTCGAATCCATTTTTTGCTCAACGTTAATGACCATCAttttt                                                                                                                                         <  2:9245/247‑1 (MQ=255)
                                                                                                              tCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCg                                                                                                                                   >  1:240089/1‑251 (MQ=255)
                                                                                                              tCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTATCGAATCCATTTTTTGCTGAACGTTAATGACCTTCATTTTTGTACCg                                                                                                                                   >  1:208774/1‑251 (MQ=255)
                                                                                                                 gCTTCCTGATTGACACTGATGGCAAAATCGAACATGTCTTTGACTATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCTCTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTc                                                                                                                                <  1:43920/251‑1 (MQ=255)
                                                                                                                            gACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCACAACCAGCAATCACAAAGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTGCTTCGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAACATCACATTCTCGCTTTACCCTGCGACACGGGCGTCGAATCCATTTTTTGCTAAACGATAATGACCATCATTTTTGTACCATTCAGattcacgtt                                                                                                                     >  2:373716/1‑243 (MQ=255)
                                                                                                                            gACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGt                                                                                                                      <  1:115333/250‑1 (MQ=255)
                                                                                                                             acgctgtaggCAAAATCGAACATGTCTTTGACGATTTCAAAACCATCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTGTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGtga                                                                                                                    <  2:78880/243‑3 (MQ=255)
                                                                                                                                                              aTTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCa                                                                                   <  1:146243/251‑1 (MQ=255)
                                                                                                                                                                 tCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGCTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGa                                                                                 <  1:7522/250‑1 (MQ=255)
                                                                                                                                                                                  cacGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGCTTGCGGCAAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGgtg                                                                                                                                                                                  <  1:64454/136‑1 (MQ=255)
                                                                                                                                                                                  cacGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTACTCCATTCCGTGCTGGCTGCGCTTGCGGCAAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGgtg                                                                                                                                                                                  >  2:64454/1‑136 (MQ=255)
                                                                                                                                                                                                  aCTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGtt                                               >  1:144320/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggctgCGTTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGcc  <  1:39048/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCTTCCTGATTGACGCTGATGGCAAAATCGAACATGTCTTTGACGATTTCAAAACCAGCAATCACCACGACGTTGTGCTGAACTGGCTGAAAGAACACGCCTGATTACTTTGCTCCATTCCGTGCTGGCTGCGCTTGCGGCCAGCATACCTCACTTCTCGTGATCAAGATCACATTCTCGCTTTCCCCTGCGACACGGGTGTCGAATCCATTTTTTGCTGAACGTTAATGACCATCATTTTTGTACCGTTCAGAATCCAGTTAATACATAACTTATTGAATATATTGAGTTAATCAGAATGGCATCCTTTATGCAATATGAAATGCAATGTTTCATATCATTTTCAAGGAGCCGACATGAACCGCTTTGTGGTGGCC  >  NC_000913/2600732‑2601221

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: