Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,446 T→G 7.7% intergenic (‑269/+128) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,4460TG 7.7% 178.2 / 7.3 52intergenic (‑269/+128)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base G (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

TGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTT  >  NC_000913/699308‑699590
                                                                                                                                          |                                                                                                                                                
tGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                           >  2:1541135/1‑146 (MQ=255)
tGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                           <  1:1541135/146‑1 (MQ=255)
    aaaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                              >  3:971906/1‑139 (MQ=255)
    aaaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                              <  4:971906/139‑1 (MQ=255)
                    ggATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                >  3:1158827/1‑121 (MQ=255)
                    ggATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                <  4:1158827/121‑1 (MQ=255)
                         tatttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGAAGCGGCGTGAACGCCt                                                                                                                                                  <  2:694824/114‑1 (MQ=255)
                         tatttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGAAGCGGCGTGAACGCCt                                                                                                                                                  >  1:694824/1‑114 (MQ=255)
                          atttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                          <  3:571448/121‑1 (MQ=255)
                          atttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                          >  4:571448/1‑121 (MQ=255)
                           tttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTc                                                                                                                      <  3:667346/140‑1 (MQ=255)
                           tttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTc                                                                                                                      >  4:667346/1‑140 (MQ=255)
                                  atgaATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGc                                                                                                      >  4:841310/1‑149 (MQ=255)
                                           tAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                                  >  3:1817887/1‑112 (MQ=255)
                                           tAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                                  <  4:1817887/112‑1 (MQ=255)
                                              tttttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCaa                                                                                                                    >  1:1382894/1‑123 (MQ=255)
                                              tttttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCaa                                                                                                                    <  2:1382894/123‑1 (MQ=255)
                                               ttttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATAAATTGCAGGAGCTGCGTAGGCCTGa                                                                                         >  2:1230909/1‑149 (MQ=255)
                                                              ttGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGa                                                                                                        <  2:1273417/119‑1 (MQ=255)
                                                              ttGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGa                                                                                                        >  1:1273417/1‑119 (MQ=255)
                                                                    gTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                                      <  2:580406/131‑1 (MQ=255)
                                                                    gTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                                      >  1:580406/1‑131 (MQ=255)
                                                                      ttATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                  <  4:910754/69‑1 (MQ=255)
                                                                      ttATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                                  >  3:910754/1‑69 (MQ=255)
                                                                       tataTTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                               >  4:379848/1‑71 (MQ=255)
                                                                       tataTTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                               <  3:379848/71‑1 (MQ=255)
                                                                             ttACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATCTATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGtt                                                           <  1:1230909/149‑1 (MQ=255)
                                                                                   gggCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcag                                                     >  2:1816857/1‑148 (MQ=255)
                                                                                     gCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcagca                                                   <  1:1816857/149‑1 (MQ=255)
                                                                                       cAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                                      >  2:685113/1‑112 (MQ=255)
                                                                                       cAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                                      <  1:685113/112‑1 (MQ=255)
                                                                                         aaTTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGa                                                                                         >  4:232293/1‑107 (MQ=255)
                                                                                         aaTTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGa                                                                                         <  3:232293/107‑1 (MQ=255)
                                                                                           ttGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCa                                             >  1:1251203/1‑149 (MQ=255)
                                                                                               tGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATc                                                                          <  4:1335040/116‑1 (MQ=255)
                                                                                               tGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATc                                                                          >  3:1335040/1‑116 (MQ=255)
                                                                                                    aaGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGTCTACAAAAGCATGAAAATTCaa                                                                                                                    <  4:985103/69‑1 (MQ=255)
                                                                                                    aaGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCaa                                                                                                                    >  3:985103/1‑69 (MQ=255)
                                                                                                          gATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTa                                                                                 >  4:1090851/1‑98 (MQ=255)
                                                                                                          gATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTa                                                                                 <  3:1090851/98‑1 (MQ=255)
                                                                                                          gATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATc                                                                          <  4:1211389/105‑1 (MQ=255)
                                                                                                          gATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATc                                                                          >  3:1211389/1‑105 (MQ=255)
                                                                                                            tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATAtt                                                                                                              >  3:1236189/1‑67 (MQ=255)
                                                                                                            tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATAtt                                                                                                              <  4:1236189/67‑1 (MQ=255)
                                                                                                             ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTg                           <  3:841310/149‑1 (MQ=255)
                                                                                                              tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTgg                          <  2:1251203/149‑1 (MQ=255)
                                                                                                                             aggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcatc                                                    <  1:523645/107‑1 (MQ=18)
                                                                                                                             aggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcatc                                                    >  2:523645/2‑108 (MQ=18)
                                                                                                                               gcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGc                                                                                                      <  1:1163959/56‑1 (MQ=255)
                                                                                                                               gcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGc                                                                                                      >  2:1163959/1‑56 (MQ=255)
                                                                                                                                     aCGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGAtt   >  4:1651772/1‑149 (MQ=255)
                                                                                                                                      cGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGAttt  >  3:1575987/1‑149 (MQ=255)
                                                                                                                                       gCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                                      <  2:56644/64‑1 (MQ=255)
                                                                                                                                       gCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATaa                                                                                      >  1:56644/1‑64 (MQ=255)
                                                                                                                                         cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCaga                                                 <  3:1177642/99‑1 (MQ=255)
                                                                                                                                         cTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCaga                                                 >  4:1177642/1‑99 (MQ=255)
                                                                                                                                          |                                                                                                                                                
TGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTT  >  NC_000913/699308‑699590

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