Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 275,065 G→A T18I (ACC→ATC)  insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA275,0650GA100.0% 57.4 / NA 24T18I (ACC→ATC) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/11);  total (13/11)

GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTAA  >  NC_000913_3_pae_tpiA/275005‑275210
                                                            |                                                                                                                                                 
gACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGTGCTGGTTCGCCCCCTCCGGAg                                                           >  2:604949/1‑149 (MQ=11)
  cAAGAAAATCTCTTTTCTTGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGAAAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGc                                                         <  1:305660/149‑1 (MQ=255)
   aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACCTGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCg                                                        >  1:15459/1‑149 (MQ=25)
   aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATAATCTCATATCAGGGACTTGTTCGCACATTCCGGAGGCg                                                        >  1:13051/1‑149 (MQ=255)
     gAAAATCTGTGTTCAGGGCTCACGCGGCTTAATGCTGAATTCACTGTCGGCGAAGATAAGGTGATGACTCATGATGATCCCGGTTGTATGTCTCCATATTACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt                                                      <  2:815299/149‑1 (MQ=11)
     gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt                                                      >  1:188927/1‑149 (MQ=37)
          tctcTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAg                                                 >  2:808208/1‑149 (MQ=255)
          tctcTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAg                                                 >  1:750232/1‑149 (MQ=255)
                        tgtcGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTTGTCAATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTTTCTCCTTCCGAAGGCGTTATGAGCTGGCGGCCCtttt                                   >  2:984468/4‑149 (MQ=17)
                        tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCGTGATGAACCCTGTTCTATGGCTCCAGATGACAAAAATGAGCTCATATCAGGGACTTGTTCGCACCTTACGGAGGCGTTATAAGCTGGCGGCCCtttt                                   >  1:346560/1‑149 (MQ=17)
                        tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGAATTGTCAGCATTGTCCGGAGGCGTTATGAGCTGGCGGCCCtttt                                   <  1:774300/149‑1 (MQ=17)
                         gACGGCGCTTACTTGTGAATTCACGGTCGGCGAAGATAAGTTGATGACCCATGCTGAACCCTGTTCTATGGCGCCAGATCACCAACATGATCTCATACCAGGGACTTGTTCGCACCTTCCGGAGGCTTTATGAGCTCGCGGCCCttttt                                  <  2:70402/149‑1 (MQ=17)
                          ccGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGGTATGAGCTGGCGGCCCTTTTTg                                 >  1:564921/2‑149 (MQ=255)
                           cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGt                                >  1:131931/1‑149 (MQ=255)
                             gctgttACTGCTGAATTCACTGGCGGCGAAGATAAGTTGATGACTTATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTAt                              <  2:374079/145‑1 (MQ=255)
                             gcgcTTACTTGTGAATTCACTGTCCGAGAAGATATTTTGGTGTTTAGTCATGAACCCTGTTCTATGGATCCAGGTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTAt                              <  1:236200/149‑1 (MQ=17)
                             gcgcTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGTCTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTAt                              <  1:938509/149‑1 (MQ=255)
                                 ttACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGa                          <  1:152888/149‑1 (MQ=255)
                                    ctgctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatta                       <  1:808208/149‑1 (MQ=255)
                                               aCTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAc            >  1:710188/1‑149 (MQ=255)
                                                 tGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAcca          >  2:745648/1‑148 (MQ=255)
                                                      gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAt     <  2:188927/149‑1 (MQ=255)
                                                        gAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTa   >  1:889659/1‑149 (MQ=255)
                                                         aaGATAAGTTGATGACTCATGATGAACCCTGTGCTATGGTTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTaa  <  2:556504/149‑1 (MQ=255)
                                                            |                                                                                                                                                 
GACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTAA  >  NC_000913_3_pae_tpiA/275005‑275210

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