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% 58.6 / NA 20*412* (TAG→TAAintRRac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/10);  total (10/10)

CGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTAT  >  NC_000913/1411873‑1412113
                                                                                                                               |                                                                                                                 
cGTTACGACGTAGACGTTTTTGTAATTTTTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                        <  1:110941/139‑1 (MQ=255)
          cAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                              >  1:507589/1‑139 (MQ=255)
                                                tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                    <  2:738413/95‑1 (MQ=255)
                                                tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                    >  1:738413/1‑95 (MQ=255)
                                                tctttGGAATTTGTTGATTTTCTTGCATTCTTTCAGGCCTCTAGTACTAAATGGGGCAAGTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                        >  2:491026/1‑139 (MQ=255)
                                                        aTTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaaa                                                <  2:719713/139‑1 (MQ=255)
                                                                aTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATac                                        <  2:507589/139‑1 (MQ=255)
                                                                              tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                  <  1:62956/67‑1 (MQ=255)
                                                                              tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                  >  2:62956/1‑67 (MQ=255)
                                                                               ttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                  <  1:577934/66‑1 (MQ=255)
                                                                               ttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                  >  2:577934/1‑66 (MQ=255)
                                                                                        ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAAcc                                                                                                         >  1:220938/1‑50 (MQ=255)
                                                                                        ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAAcc                                                                                                         <  2:220938/50‑1 (MQ=255)
                                                                                                 aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaaa                                                >  2:7707/1‑98 (MQ=255)
                                                                                                 aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaaa                                                <  1:7707/98‑1 (MQ=255)
                                                                                                      ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                   >  2:461155/1‑74 (MQ=255)
                                                                                                      ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                   <  1:461155/74‑1 (MQ=255)
                                                                                                      ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat  >  1:61152/1‑139 (MQ=255)
                                                                                                           aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                    >  2:101220/1‑36 (MQ=255)
                                                                                                           aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                    <  1:101220/36‑1 (MQ=255)
                                                                                                                               |                                                                                                                 
CGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTAT  >  NC_000913/1411873‑1412113

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