Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,645,063 G→A intergenic (‑161/+56) rem ← / ← hokD Qin prophage, uncharacterized protein/Qin prophage, small toxic polypeptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,645,0630GA100.0% 78.2 / NA 26intergenic (‑161/+56)rem/hokDQin prophage, uncharacterized protein/Qin prophage, small toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/12);  total (14/12)

CCCCTTTAAGGGACCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGTTCGGATTCGTACCTCGCAGA  >  NC_000913/1644929‑1645189
                                                                                                                                      |                                                                                                                              
ccccTTTAGGGAACCTTCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGa                                                                                                                            <  1:503021/139‑1 (MQ=255)
 cccTTTAAGGGACCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAg                                                                                                                           >  1:1237152/1‑139 (MQ=255)
              cATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCcaca                                                                                                              >  1:296294/1‑139 (MQ=255)
              cATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCcaca                                                                                                              >  2:297195/1‑139 (MQ=255)
                      ttATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCATCGTCCAGAACGTCAAACGCGTTGAAAATAAAAACAATAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtg                                                                                                                       <  2:191106/122‑1 (MQ=255)
                      ttATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtg                                                                                                                       >  1:191106/1‑122 (MQ=255)
                                gcgcTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGAt                                                                                            >  2:373335/1‑139 (MQ=255)
                                gcgcTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGAt                                                                                            >  1:340741/1‑139 (MQ=255)
                                                tttGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtgcc                                                                                                                     <  2:31989/98‑1 (MQ=255)
                                                tttGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtgcc                                                                                                                     >  1:31989/1‑98 (MQ=255)
                                                 ttGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGtc                                                                           <  2:747344/139‑1 (MQ=255)
                                                    aTTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtgcctgc                                                                                                                  <  1:1064705/97‑1 (MQ=255)
                                                    aTTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtgcctgc                                                                                                                  >  2:1064705/1‑97 (MQ=255)
                                                                    ggccAGAACTGCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtg                                                                                                                       <  2:1157679/74‑1 (MQ=38)
                                                                    gTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGa                                                                                                                         >  2:208770/1‑74 (MQ=255)
                                                                    gTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGa                                                                                                                         <  1:208770/74‑1 (MQ=255)
                                                                    gTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGAtg                                                                                                                       >  1:1157679/1‑76 (MQ=255)
                                                                    gTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGt                                                        >  1:1007949/1‑139 (MQ=255)
                                                                                   cgcgTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAg                                         <  2:1007949/139‑1 (MQ=255)
                                                                                                  aacaaaaacCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACa                                          >  1:511917/1‑123 (MQ=255)
                                                                                                  aacaaaaacCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACa                                          <  2:511917/123‑1 (MQ=255)
                                                                                                       aaacCCGCCGAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGtt                     <  1:373335/139‑1 (MQ=255)
                                                                                                             gccgAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGTTCg                   >  1:449757/1‑135 (MQ=255)
                                                                                                             gccgAAGCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGTTCg                   <  2:449757/135‑1 (MQ=255)
                                                                                                                   gCGGGTTAAGTGCGGGTGCATTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGTTCGGATTCGTAcc         >  2:246040/1‑139 (MQ=255)
                                                                                                                         tAAGTGCGGGTGCATTGAGGATGCCTGCCCATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGTTCGGATTCGTACCTCGCaga  <  2:1237152/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                              
CCCCTTTAAGGGACCATCCCTCTTATCCCTGCGCGCTACTTAAGTATTTTTGATTCTATTCCGGCACCGTCCAGAACTTCAAACGCGTTGAAAATAAAAACAAAAACCCGCCGAAGCGGGTTAAGTGCGGGTGCGTTGAGGATGCCTGCCACATCAGAGGTGGCGAGGGATTTCTCCCCCGCCGGGTCTCTTACTCCTCAGGTTCGTAAGCTGTGAAGACAGCGACCTCCGTCTGGCCGGTTCGGATTCGTACCTCGCAGA  >  NC_000913/1644929‑1645189

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