Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,087 Δ2 bp coding (942‑943/990 nt) rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,0860C.100.0% 76.0 / NA 17coding (944/990 nt)rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/8);  total (9/8)
*NC_000913_3_pae_pgi3,440,0870C.100.0% 76.0 / NA 17coding (943/990 nt)rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/8);  total (9/8)

TAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATA  >  NC_000913_3_pae_pgi/3439942‑3440208
                                                                                                                                                ||                                                                                                                         
tAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGcagaga                                                                                                                       <  3:47139/150‑6 (MQ=255)
                   gTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGc                                                                                                  >  1:111949/1‑150 (MQ=255)
                     tGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCAc                                                                                                <  2:111949/150‑1 (MQ=255)
                               ttGACGACCACTCCTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                      >  2:18327/1‑150 (MQ=255)
                                  acgacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                   <  3:26228/150‑1 (MQ=255)
                                    gacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGt                                                                                 >  4:74451/1‑150 (MQ=255)
                                          tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAgaga                                                                           <  3:74451/150‑1 (MQ=255)
                                                                     tctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                <  1:18327/150‑1 (MQ=255)
                                                                         aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCt                                                                                            >  4:19443/1‑102 (MQ=255)
                                                                         aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCt                                                                                            <  3:19443/102‑1 (MQ=255)
                                                                         aGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGc                                            >  3:115773/1‑150 (MQ=255)
                                                                                      ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                   >  3:36702/1‑98 (MQ=255)
                                                                                      ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                   <  4:36702/98‑1 (MQ=255)
                                                                                      ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                >  2:55830/1‑133 (MQ=255)
                                                                                      ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                <  1:55830/133‑1 (MQ=255)
                                                                                                                   ttGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATa  >  1:126928/1‑150 (MQ=255)
                                                                                                                               aGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa             >  2:176157/1‑127 (MQ=255)
                                                                                                                               aGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa             <  1:176157/127‑1 (MQ=255)
                                                                                                                                                ||                                                                                                                         
TAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATA  >  NC_000913_3_pae_pgi/3439942‑3440208

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