Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,412,000 C→T *412* (TAG→TAA intR ← Rac prophage, integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,412,0000CT100.0% 81.0 / NA 26*412* (TAG→TAAintRRac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

CGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTG  >  NC_000913/1411873‑1412133
                                                                                                                               |                                                                                                                                     
cGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                            <  2:989478/139‑1 (MQ=255)
                           tGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATgg                                                                                                              >  2:707690/1‑126 (MQ=255)
                           tGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATgg                                                                                                              <  1:707690/126‑1 (MQ=255)
                            gTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTtcgtcg                                                                                                <  2:187763/139‑1 (MQ=255)
                                            tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATcc                                                                                >  2:803295/1‑139 (MQ=255)
                                                                                      ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                                  >  1:21269/1‑63 (MQ=255)
                                                                                      ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                                  <  2:21269/63‑1 (MQ=255)
                                                                                        ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                                       <  2:725977/120‑1 (MQ=255)
                                                                                        ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                                       >  1:725977/1‑120 (MQ=255)
                                                                                        ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGcccc                                    >  2:803369/1‑139 (MQ=255)
                                                                                         tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                              >  2:1435/1‑80 (MQ=255)
                                                                                         tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                              >  2:195780/1‑80 (MQ=255)
                                                                                         tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                              <  1:195780/80‑1 (MQ=255)
                                                                                         tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                              <  1:1435/80‑1 (MQ=255)
                                                                                              cTAAATTGTGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCCCATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                              <  2:916970/139‑1 (MQ=255)
                                                                                              cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                              <  1:803369/139‑1 (MQ=255)
                                                                                              cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                              >  2:34104/1‑139 (MQ=255)
                                                                                                 aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATc                                                                                 >  2:418688/1‑85 (MQ=255)
                                                                                                 aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATc                                                                                 <  1:418688/85‑1 (MQ=255)
                                                                                                 aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                              <  1:1109369/136‑1 (MQ=255)
                                                                                                 aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                              >  2:1109369/1‑136 (MQ=255)
                                                                                                     gggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                       >  2:210322/1‑75 (MQ=255)
                                                                                                     gggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                       <  1:210322/75‑1 (MQ=255)
                                                                                                                gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                       >  2:1117695/1‑64 (MQ=255)
                                                                                                                gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                       <  1:1117695/64‑1 (MQ=255)
                                                                                                                          tGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTg  <  1:697751/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                     
CGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTCCTG  >  NC_000913/1411873‑1412133

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