Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,633,139 A→G 7.7% D23G (GAT→GGT)  ydfK → cold shock protein, function unknown, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,633,1390AG 7.7% 137.7 / 3.7 52D23G (GAT→GGT) ydfKcold shock protein, function unknown, Qin prophage
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base G (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

AATCAAAAGACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGACAACCAAAAT  >  NC_000913/1633076‑1633287
                                                               |                                                                                                                                                    
aaTCAAAAGACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGaaaaaa                                                                 <  7:327847/149‑1 (MQ=31)
   cAAAAGACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCa                                                              >  3:177314/1‑149 (MQ=32)
    aaaaGACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGGTGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCAt                                                             >  1:504238/1‑149 (MQ=17)
    aaaaGACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGGTGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCAt                                                             >  7:34166/1‑149 (MQ=17)
        gACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCt                                                         <  3:473877/149‑1 (MQ=32)
             ccTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGAt                                                     >  5:502204/1‑148 (MQ=32)
             ccTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGAt                                                     >  4:69697/1‑148 (MQ=32)
             ccTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGAt                                                     >  8:197542/1‑148 (MQ=32)
               tAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATgg                                                           >  8:460554/1‑140 (MQ=32)
               tAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATgg                                                           <  7:460554/140‑1 (MQ=32)
                 taGTGGTTCCCTTCGAAGCTTCC‑GAAGTAAGAATTATCCTAGGTGATGCTGTAGTGGAAGTAGCAAAACAGGTAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACaaa                                               <  3:161958/148‑1 (MQ=11)
                 aaGTTGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAAc                                                  <  1:295107/147‑1 (MQ=17)
                 aaGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACa                                                 >  5:60649/1‑148 (MQ=32)
                   gaggTTCCCTGCGCAGCTT‑CAGAAGTAAGAATTATCCTAGGGGAAGCAGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATAACATTGAAGCAAACATAAAGAAAAAATCATTGCAGGATAACaaa                                               >  4:161958/3‑147 (MQ=11)
                   gTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACa                                                 >  2:116848/1‑146 (MQ=32)
                       ttCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGCat                                           >  8:552232/1‑148 (MQ=17)
                         cccTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGGTGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAatta                                         <  2:504238/148‑1 (MQ=17)
                         cccTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGGTGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAatta                                         <  8:34166/148‑1 (MQ=17)
                         cccTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAatta                                         <  3:115549/148‑1 (MQ=32)
                           cTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAAttatta                                      <  5:158936/149‑1 (MQ=32)
                            tGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAAttatta                                      >  5:543836/1‑148 (MQ=32)
                            tGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAAttatta                                      >  3:233345/1‑148 (MQ=32)
                            tGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAAttatta                                      >  2:354395/1‑148 (MQ=32)
                              cgcAGCTTCCTGAAGTAAGATTTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACTTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTAc                                     <  2:534350/147‑1 (MQ=17)
                              cgcAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTAc                                     <  6:543836/147‑1 (MQ=32)
                               gcAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAattatt                                       <  2:361546/144‑1 (MQ=32)
                               gcAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAattatt                                       >  1:361546/1‑144 (MQ=32)
                               gcAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACaaa                                  <  4:177314/149‑1 (MQ=32)
                                 aGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACa                                <  4:436981/149‑1 (MQ=32)
                                  gCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAg                               >  1:449835/1‑149 (MQ=32)
                                   cTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGc                              <  6:60649/149‑1 (MQ=32)
                                   cTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGc                              >  5:28847/1‑149 (MQ=32)
                                      ccTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGa                            >  8:213607/1‑148 (MQ=32)
                                      ccTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGa                            <  5:524275/148‑1 (MQ=32)
                                       cTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTTCATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGc                              >  3:55113/1‑145 (MQ=17)
                                       cTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGc                              <  4:55113/145‑1 (MQ=32)
                                        tGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGatat                         >  6:47309/1‑149 (MQ=32)
                                         gAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAg                                              <  5:551381/127‑1 (MQ=32)
                                         gAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAg                                              >  6:551381/1‑127 (MQ=32)
                                          aaGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGatat                         <  6:8059/147‑1 (MQ=32)
                                          aaGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCa                       <  7:197542/149‑1 (MQ=32)
                                             tAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGtt                    >  3:75272/1‑149 (MQ=32)
                                               aGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGCCCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGGTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCt                  >  2:318554/1‑149 (MQ=17)
                                               aGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTc                   <  4:233345/148‑1 (MQ=32)
                                               aGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTc                   <  1:116848/148‑1 (MQ=32)
                                               aGAATTATCATATTGGATGCTGTAGTGTAATTATCAAAAAAGTGAATACCTATAAATACCAGAACATTGCTTGATTACATTGAAGTAAACATAAAGAAAAAATCATGGCTGTATAACAAATAATTATTACAAACAGCGATATCAGTTc                   <  8:125851/148‑1 (MQ=11)
                                                   ttATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTaaa              <  3:225691/149‑1 (MQ=32)
                                                   ttATCCTAGGGGATACTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTaaa              <  4:215790/149‑1 (MQ=17)
                                                      tCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGAc           <  1:187387/149‑1 (MQ=32)
                                                             ggATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGACAACCaaa    <  2:449835/149‑1 (MQ=32)
                                                              gATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGACAACCaaaa   <  6:28847/149‑1 (MQ=32)
                                                              gATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGACAACCaaaa   <  2:397518/149‑1 (MQ=32)
                                                               aTGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGACAACCAAAAt  <  1:354395/149‑1 (MQ=32)
                                                               |                                                                                                                                                    
AATCAAAAGACACCCTAAAGTGGTTCCCTGCGCAGCTTCCTGAAGTAAGAATTATCCTAGGGGATGCTGTAGTGGAAGTAGCAAAACAGGGAAGACCTATCAATACCAGAACATTGCTTGATTACATTGAAGGAAACATAAAGAAAAAATCATGGCTGGATAACAAAGAATTATTACAAACAGCGATATCAGTTCTTAAAGACAACCAAAAT  >  NC_000913/1633076‑1633287

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