Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,639,030 C→T *166* (TAG→TAA rzpQ ← Rz‑like protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,639,0300CT100.0% 77.4 / NA 25*166* (TAG→TAArzpQRz‑like protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/13);  total (12/13)

TTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGT  >  NC_000913/1638898‑1639154
                                                                                                                                    |                                                                                                                            
ttaATATAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGTAATTTTGGCTCATGGAAGATTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACAtt                                                                                                                        <  2:164398/139‑1 (MQ=255)
 taataaAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAGGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACAttt                                                                                                                       >  2:608535/1‑139 (MQ=255)
         cGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCttt                                                                                                               >  2:145595/1‑139 (MQ=255)
             tCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCttttt                                                                                                             >  1:405051/1‑137 (MQ=255)
             tCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCttttt                                                                                                             <  2:405051/137‑1 (MQ=255)
                     aTTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTAt                                                                                                           >  1:454136/1‑131 (MQ=255)
                     aTTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTAt                                                                                                           <  2:454136/131‑1 (MQ=255)
                                    cTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGt                                                                                    >  1:709131/1‑139 (MQ=255)
                                                   tCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCt                                                                     >  1:676341/1‑139 (MQ=255)
                                                    cAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTg                                                                    <  1:68917/139‑1 (MQ=255)
                                                                                                    ccGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTc                                                                       >  1:151399/1‑88 (MQ=255)
                                                                                                    ccGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTc                                                                       <  2:151399/88‑1 (MQ=255)
                                                                                                    ccGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGa                              <  1:141554/129‑1 (MQ=255)
                                                                                                    ccGCGATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGa                              >  2:141554/1‑129 (MQ=255)
                                                                                                      gcgATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGa                                                                                 >  2:113906/1‑76 (MQ=255)
                                                                                                      gcgATGTAAATGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGa                                                                                 <  1:113906/76‑1 (MQ=255)
                                                                                                              aaTGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTcc          <  2:676341/139‑1 (MQ=255)
                                                                                                                tGTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGa        <  1:608535/139‑1 (MQ=255)
                                                                                                                 gTTCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAg       <  2:640991/139‑1 (MQ=255)
                                                                                                                  ttCCGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGc      <  1:145595/139‑1 (MQ=255)
                                                                                                                    ccGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTTTGCTGCAATATCTCCGAGcct    >  2:852022/1‑137 (MQ=255)
                                                                                                                     cGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATg   <  1:802941/139‑1 (MQ=255)
                                                                                                                     cGATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATg   <  2:4355/139‑1 (MQ=255)
                                                                                                                      gATTTAATATTACCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGt  >  1:699020/1‑139 (MQ=255)
                                                                                                                                 aCCTTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTcc          >  2:153528/1‑120 (MQ=255)
                                                                                                                                    |                                                                                                                            
TTAATAAAACGTGTCATTCTGATTAAGACCTTTTATCTTACCCTTAAGATTTCAGGAATTTTGGCTCATGGAAGAGTCCTTTTTATTTAAATTTTACATTCCGCGATGTAAATGTTCCGATTTAATATTACCCTACATTTGATGCTTTTTATCTCTTAAAGATTCATAGATCTGTTGACAAGTCACTCCTGCGATGTAGCGTTCGTCAGCAATTTCAGCATAAAGCTGAGCTTCTGCTGCAATATCTCCGAGCATGT  >  NC_000913/1638898‑1639154

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: