Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,412,027 T→C I403M (ATA→ATG intR ← Rac prophage, integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,412,0270TC100.0% 112.0 / NA 30I403M (ATA→ATGintRRac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/14);  total (16/14)

TCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGCAAAAGTATGC  >  NC_000913/1411903‑1412162
                                                                                                                            |                                                                                                                                       
tCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGtt                                                                                                                           <  1:153724/139‑1 (MQ=255)
              tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATcc                                                                                                             >  2:85424/1‑139 (MQ=255)
                 ttctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCcaca                                                                                                          <  1:151702/139‑1 (MQ=255)
                  tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGtt                                                                                                                           <  1:750615/121‑1 (MQ=255)
                  tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGtt                                                                                                                           >  2:750615/1‑121 (MQ=255)
                     ttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATAcc                                                                                                      >  2:651253/1‑139 (MQ=255)
                       gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCa                                                                                                                         >  2:606422/1‑118 (MQ=255)
                       gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCa                                                                                                                         <  1:606422/118‑1 (MQ=255)
                                               attTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCAt                                                                                                                        <  1:379253/95‑1 (MQ=255)
                                               attTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCAt                                                                                                                        >  2:379253/1‑95 (MQ=255)
                                                        ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCca                                                                                                            >  2:64995/1‑98 (MQ=255)
                                                        ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCca                                                                                                            <  1:64995/98‑1 (MQ=255)
                                                           tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCa                                                                >  2:527057/1‑139 (MQ=255)
                                                           tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCa                                                                >  1:715200/1‑139 (MQ=255)
                                                                cTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTc                                                                        >  2:337651/1‑126 (MQ=255)
                                                                cTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCAGCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTc                                                                        <  1:337651/126‑1 (MQ=255)
                                                                 tAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggc                                                          <  1:439005/139‑1 (MQ=255)
                                                                       gggCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCta                                                    >  1:95022/1‑139 (MQ=255)
                                                                        ggCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAttt                                                                                              <  1:706571/96‑1 (MQ=255)
                                                                        ggCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAttt                                                                                              >  2:706571/1‑96 (MQ=255)
                                                                        ggCAAATTGGGGACAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                                                   >  1:158560/1‑139 (MQ=255)
                                                                                    gCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGTTGGCTATAAATGCCGGGtt                                       <  2:715200/139‑1 (MQ=255)
                                                                                         cTTTGCAACTACGATAACCGCGCATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTc                                  <  1:527057/139‑1 (MQ=255)
                                                                                                          ccGCGCATTCAACAGGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATGTGCGCAGTTTCATGCCCCATTTGGCTCGCTATAAATGCCGGGTTCGCTCCTGCCGTCAAcagccag                 >  1:476095/1‑139 (MQ=255)
                                                                                                             cgcATTCAACATGGCCATCTGTTCGTCGTTCATGtcatc                                                                                                                  <  2:476240/39‑1 (MQ=255)
                                                                                                             cgcATTCAACATGGACATCTGTTCGTCGTTCATGtcatc                                                                                                                  >  1:476240/1‑39 (MQ=255)
                                                                                                               cATTCAACATGGCCATCTGTTCGTCGTTCATGtcatc                                                                                                                  <  2:510159/37‑1 (MQ=255)
                                                                                                               cATTCAACATGGCCATCTGTTCGTCGTTCATGtcatc                                                                                                                  >  1:510159/1‑37 (MQ=255)
                                                                                                               cATTCAACATGGCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGc            >  2:217626/1‑139 (MQ=255)
                                                                                                                         ggCCATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGCAAAAGTATGc  <  2:95022/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                                       
TCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTGCCGTCAACAGCCAGCAGGCAAAAGTATGC  >  NC_000913/1411903‑1412162

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