Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 365,410 G→T P299Q (CCG→CAG)  lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913365,4100GT100.0% 34.0 / NA 12P299Q (CCG→CAG) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (3/9);  total (3/9)

CAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCC  >  NC_000913/365280‑365539
                                                                                                                                  |                                                                                                                                 
cagACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGa                                                                                                                          <  1:96834/140‑1 (MQ=255)
    ccATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGtt                                                                                                                      <  1:318778/140‑1 (MQ=255)
                       cgcgGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGa                                                                                                   <  2:92372/140‑1 (MQ=255)
                       cgcgGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGa                                                                                                   <  1:343597/140‑1 (MQ=255)
                       cgcgGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGa                                                                                                   <  1:370424/140‑1 (MQ=255)
                        gcgGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGAt                                                                                                  >  2:380755/1‑140 (MQ=255)
                                                                                    aaCCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGa                                      >  1:290088/1‑140 (MQ=255)
                                                                                     cccaccGCAGGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCCCGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGAc                                     <  2:67002/139‑1 (MQ=255)
                                                                                      ccaccGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGAcc                                    >  2:299625/1‑140 (MQ=255)
                                                                                       caccGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCt                                   <  1:299625/140‑1 (MQ=255)
                                                                                                 tAGAGATTCGGGATTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCCCGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCAcc                         <  1:187579/140‑1 (MQ=255)
                                                                                                 aagagATTGGGGTTTTCGGCGCTCCACAGTTTCTGTTTTTCGACGTTAAGACGTAGTGTGACGGGATCGGCATACCCCCCCCGCTCATCGATAATTCCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCAcc                         <  2:217707/139‑1 (MQ=255)
                                                                                                  agagATTCGGGATTTCGGCGCTCCACAGTTTCTGTTTTTCGACGTTCAGACGTAGTGTGCCGCGATCGGCATACCCCCCCCGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCAccc                        <  2:268040/140‑1 (MQ=255)
                                                                                                             aTTTCGGCGCTCCACAGTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGa             <  1:346680/140‑1 (MQ=255)
                                                                                                                       cccACATTTTCTGGTTTTCGACGTTCAGACGTAGTGGGACGCGATGGGCATAACCCCCCCCCTCATCGATATTTTCCCCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCCCCCTGCCATAAAGAAACTGTTAcc   <  2:145337/139‑1 (MQ=255)
                                                                                                                       cccACAGTTTCTGTTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTCCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTAcc   <  2:128425/139‑1 (MQ=255)
                                                                                                                        ccACATTTTCTGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCAAACCCCCCCCGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGTGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTAccc  <  2:50829/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
CAGACCATTTTCAATCCGCACCTCGCGGAAACCGACATCGCAGGCTTCTGCTTCAATCAGCGTGCCGTCGGCGGTGTGCAGTTCAACCACCGCACGATAGAGATTCGGGATTTCGGCGCTCCACAGTTTCGGGTTTTCGACGTTCAGACGTAGTGTGACGCGATCGGCATAACCACCACGCTCATCGATAATTTCACCGCCGAAAGGCGCGGTGCCGCTGGCGACCTGCGTTTCACCCTGCCATAAAGAAACTGTTACCC  >  NC_000913/365280‑365539

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