Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,876,374 G→T 20.0% intergenic (‑47/‑205) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,3740GT20.0% 76.5 / 14.7 30intergenic (‑47/‑205)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (3/3);  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 = 7.66e-01

ACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTC  >  NC_000913/2876094‑2876652
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                      
aCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACtata                                                                                                                                                                                                                                                                                  <  2:416675/287‑1 (MQ=255)
    tCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGc                                                                                                                                                                                                                                                                            >  1:334439/1‑289 (MQ=255)
                         cgGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAAtt                                                                                                                                                                                                                                                        >  1:449532/1‑288 (MQ=255)
                             aTTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCaa                                                                                                                                                                                                                                                   <  1:474757/289‑1 (MQ=255)
                                         tATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                                                                                                                                                                                                                                        <  2:449532/288‑1 (MQ=255)
                                                         cGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa                                                                                                                                                                                                                       >  2:313694/1‑289 (MQ=255)
                                                         cGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa                                                                                                                                                                                                                       >  2:429977/1‑289 (MQ=255)
                                                                     cGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                                                                                                                                                           >  2:77629/1‑289 (MQ=255)
                                                                     cGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                                                                                                                                                           >  2:351900/1‑289 (MQ=255)
                                                                      gTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGc                                                                                                                                                                                                             <  1:429977/286‑1 (MQ=255)
                                                                                    cGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                                                                                                                            <  2:334439/289‑1 (MQ=255)
                                                                                        ttGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                                                                                                                                                                                                                                      >  2:490605/1‑243 (MQ=255)
                                                                                        ttGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                                                                                                                                                                                                                                      <  1:490605/243‑1 (MQ=255)
                                                                                                          aTGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                                                                                                                      <  1:77629/289‑1 (MQ=255)
                                                                                                                               aTAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCCCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa                                                                                                                                                 >  2:78763/1‑289 (MQ=255)
                                                                                                                                             cACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                                                                                                                >  1:294467/1‑276 (MQ=255)
                                                                                                                                             cACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                                                                                                                <  2:294467/276‑1 (MQ=255)
                                                                                                                                                                                taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                                                                                                                           >  1:481425/1‑262 (MQ=255)
                                                                                                                                                                                taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                                                                                                                           <  2:481425/262‑1 (MQ=255)
                                                                                                                                                                                 aatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                               <  2:109565/289‑1 (MQ=255)
                                                                                                                                                                                         tGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTAttt                                                                                       <  1:78763/289‑1 (MQ=255)
                                                                                                                                                                                                                   tGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCa                                                                <  1:351900/286‑1 (MQ=255)
                                                                                                                                                                                                                              tATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTa                                                   >  2:366497/1‑288 (MQ=255)
                                                                                                                                                                                                                                         tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                               >  2:417248/1‑233 (MQ=255)
                                                                                                                                                                                                                                         tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                               <  1:417248/233‑1 (MQ=255)
                                                                                                                                                                                                                                               ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCg                                 >  2:210661/1‑289 (MQ=255)
                                                                                                                                                                                                                                                     tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGTAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTAt                           >  1:39524/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               cTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATtc                  >  1:100853/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                           cATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGc     <  1:313694/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTc  <  2:100853/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                      
ACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTC  >  NC_000913/2876094‑2876652

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