Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,876,376 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,3760GT20.0% 77.7 / 14.8 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/2876096‑2876654
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                      
aCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACtata                                                                                                                                                                                                                                                                                  <  1:416675/287‑1 (MQ=255)
    tCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGc                                                                                                                                                                                                                                                                            >  2:334439/1‑289 (MQ=255)
                         cgGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAAtt                                                                                                                                                                                                                                                        >  2:449532/1‑288 (MQ=255)
                             aTTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCaa                                                                                                                                                                                                                                                   <  2:474757/289‑1 (MQ=255)
                                         tATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                                                                                                                                                                                                                                        <  1:449532/288‑1 (MQ=255)
                                                         cGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa                                                                                                                                                                                                                       >  1:429977/1‑289 (MQ=255)
                                                         cGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa                                                                                                                                                                                                                       >  1:313694/1‑289 (MQ=255)
                                                                     cGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                                                                                                                                                           >  1:351900/1‑289 (MQ=255)
                                                                     cGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                                                                                                                                                           >  1:77629/1‑289 (MQ=255)
                                                                      gTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGc                                                                                                                                                                                                             <  2:429977/286‑1 (MQ=255)
                                                                                    cGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                                                                                                                            <  1:334439/289‑1 (MQ=255)
                                                                                        ttGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                                                                                                                                                                                                                                      >  1:490605/1‑243 (MQ=255)
                                                                                        ttGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                                                                                                                                                                                                                                      <  2:490605/243‑1 (MQ=255)
                                                                                                          aTGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                                                                                                                      <  2:77629/289‑1 (MQ=255)
                                                                                                                               aTAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCCCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGaaa                                                                                                                                                 >  1:78763/1‑289 (MQ=255)
                                                                                                                                             cACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                                                                                                                <  1:294467/276‑1 (MQ=255)
                                                                                                                                             cACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                                                                                                                >  2:294467/1‑276 (MQ=255)
                                                                                                                                                                                taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                                                                                                                           >  2:481425/1‑262 (MQ=255)
                                                                                                                                                                                taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                                                                                                                           <  1:481425/262‑1 (MQ=255)
                                                                                                                                                                                 aatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                               <  1:109565/289‑1 (MQ=255)
                                                                                                                                                                                         tGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTAttt                                                                                       <  2:78763/289‑1 (MQ=255)
                                                                                                                                                                                                                   tGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCa                                                                <  2:351900/286‑1 (MQ=255)
                                                                                                                                                                                                                              tATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTa                                                   >  1:366497/1‑288 (MQ=255)
                                                                                                                                                                                                                                         tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                               <  2:417248/233‑1 (MQ=255)
                                                                                                                                                                                                                                         tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa                                                                                               >  1:417248/1‑233 (MQ=255)
                                                                                                                                                                                                                                               ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCg                                 >  1:210661/1‑289 (MQ=255)
                                                                                                                                                                                                                                                     tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGTAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTAt                           >  2:39524/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               cTATTTTCTTGCCATCAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATtc                  >  2:100853/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                           cATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGc     <  2:313694/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cAGATGTTTCTACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTc  <  1:100853/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                      
ACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCAGATGCAGCATGACGCTGGAATCTTTACCGATAGAGTAGAGCATCACCGGATTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATGGTATTTAAGGACTCACTATGTTTTCCGCATTGCGCCACCGTACCGCTGCCCTGGCGCTCGGCGTATGCTTTATTCTCCCCGTACACGCCTC  >  NC_000913/2876096‑2876654

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: