Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,392,826 T→C Q18Q (CAA→CAG uspE ← universal stress protein UspE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,392,8260TC100.0% 55.9 / NA 18Q18Q (CAA→CAGuspEuniversal stress protein UspE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/16);  total (2/16)

TTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGC  >  NZ_CP009273/1392739‑1392909
                                                                                       |                                                                                   
ttCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATgctggctg                                                                                   >  2:340324/1‑90 (MQ=255)
       gagaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTcc                                                                            <  1:579769/90‑1 (MQ=255)
                      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                                             >  1:273694/1‑90 (MQ=255)
                           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                                        <  1:233391/90‑1 (MQ=255)
                                cTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCa                                                   <  2:273694/90‑1 (MQ=255)
                                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTa                                                                      >  2:655536/1‑70 (MQ=255)
                                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTa                                                                      <  1:655536/70‑1 (MQ=255)
                                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                                                  <  1:218617/90‑1 (MQ=255)
                                   tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                <  2:313651/90‑1 (MQ=255)
                                      ttttGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCt                                             <  2:546801/90‑1 (MQ=255)
                                           ccaccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATAc                                        <  2:433371/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:488801/90‑1 (MQ=255)
                                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                                      <  1:95044/90‑1 (MQ=255)
                                              ccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATa                                     <  1:468373/90‑1 (MQ=255)
                                               caATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAg                                    <  2:124661/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCGTACAGGGtc                        <  2:46810/90‑1 (MQ=255)
                                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc                        <  1:342061/90‑1 (MQ=255)
                                                                                 tGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGc  <  1:444316/90‑1 (MQ=255)
                                                                                 tGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGc  <  1:421530/90‑1 (MQ=255)
                                                                                       |                                                                                   
TTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGC  >  NZ_CP009273/1392739‑1392909

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