Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,900,559 C→A E343* (GAA→TAA)  yoaE ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,5590CA100.0% 107.5 / NA 30E343* (GAA→TAA) yoaEputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/20);  total (10/20)

GAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCC  >  NC_000913/1900483‑1900634
                                                                            |                                                                           
gAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATcc                                                                         <  2:347687/81‑1 (MQ=255)
gAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATcc                                                                         <  2:347686/81‑1 (MQ=255)
  aGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAg                                                                       <  1:347688/81‑1 (MQ=255)
   gCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGt                                                                      <  1:347689/81‑1 (MQ=255)
      aTCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCa                                                                   >  2:347747/1‑81 (MQ=255)
         gccgccACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACcg                                                                <  1:347690/81‑1 (MQ=255)
            gccACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCTCGa                                                             >  2:347749/1‑81 (MQ=255)
               acaTCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACAt                                                          <  2:347691/81‑1 (MQ=255)
                    aaCGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCg                                                      <  1:347692/80‑1 (MQ=255)
                      cGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGa                                                   <  1:347693/81‑1 (MQ=255)
                            ctTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGAc                                             >  1:347729/1‑81 (MQ=255)
                               ccAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACt                                            <  2:347694/79‑1 (MQ=255)
                                cAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACtgtg                                         >  2:347731/1‑81 (MQ=255)
                                cAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACtgt                                          <  1:347695/80‑1 (MQ=255)
                                  gcgcCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCg                                       <  1:347696/81‑1 (MQ=255)
                                     ccaccaGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTg                                    <  2:347697/81‑1 (MQ=255)
                                       accaGCAGTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGaa                                  <  2:347698/81‑1 (MQ=255)
                                        ccggcagTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAag                                 <  2:347699/78‑1 (MQ=255)
                                             agTTCTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAagagag                             >  1:347739/1‑80 (MQ=255)
                                                 cTTTAGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGtt                        >  2:347757/1‑81 (MQ=255)
                                                     aGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCt                     <  2:347700/80‑1 (MQ=255)
                                                     aGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTc                    >  1:347751/1‑81 (MQ=255)
                                                     aGCACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTc                    >  2:347750/1‑81 (MQ=255)
                                                       cACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTcgc                  <  2:347702/81‑1 (MQ=255)
                                                       cACGTACAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTcgc                  <  1:347701/81‑1 (MQ=255)
                                                            aCAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGAtt             >  2:347752/1‑81 (MQ=255)
                                                             cAATACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGAtt             <  1:347703/80‑1 (MQ=255)
                                                                tACCGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCAt         <  1:347704/81‑1 (MQ=255)
                                                                  ccGATGATTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATc        <  1:347705/80‑1 (MQ=255)
                                                                       gatTTAATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGAcccc  >  1:347760/1‑81 (MQ=255)
                                                                            |                                                                           
GAAGCAATCGCCGCCACATCAACGCCCTCTTCCAGCGCCACCAGCAGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCC  >  NC_000913/1900483‑1900634

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