Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,726,942	C→A	Q140H (CAG→CAT)	yiaA ←	conserved inner membrane protein YiaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,726,942	0	C	A	100.0%	133.3 / NA	39	Q140H (CAG→CAT)	yiaA	conserved inner membrane protein YiaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/18);  total (21/18)

TTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATT  >  NC_000913/3726806‑3727076                                                                                                                                         |                                                                                                                                       tttGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATgg                                                                                                                                     >  1:662038/1‑140 (MQ=255)                cAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCt                                                                                                                      >  2:246062/1‑140 (MQ=255)                      gTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTtggt                                                                                                                <  1:298752/140‑1 (MQ=255)                      gTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTtggt                                                                                                                <  1:474141/140‑1 (MQ=255)                       tttGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTtggtg                                                                                                               >  2:572235/1‑140 (MQ=255)                             aCTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTa                                                                                                         >  2:194335/1‑140 (MQ=255)                                   tATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCgg                                                                                                   <  1:246062/140‑1 (MQ=255)                                    aTGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGc                                                                                                  <  2:282473/140‑1 (MQ=255)                                         ggTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCAc                                                                                             >  1:229736/1‑140 (MQ=255)                                          gTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACg                                                                                            >  1:160604/1‑140 (MQ=255)                                           tACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                                                                <  1:440011/135‑1 (MQ=255)                                           tACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                                                                >  2:440011/1‑135 (MQ=255)                                              aaCGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAata                                                                                        <  1:673714/140‑1 (MQ=255)                                                 gAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAttt                                                                                     >  1:8712/1‑140 (MQ=255)                                                    tCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTc                                                                                  >  2:542047/1‑140 (MQ=255)                                                    tCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTc                                                                                  >  1:208082/1‑140 (MQ=255)                                                      gACTGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCCCGAATATTTTTCTg                                                                                <  2:198251/140‑1 (MQ=255)                                                         aGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCAc                                                                             <  2:208082/140‑1 (MQ=255)                                                              cAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTa                                                                        <  2:79459/140‑1 (MQ=255)                                                                     gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCAcc                                                                 >  1:416379/1‑140 (MQ=255)                                                                        aGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                                                              <  1:572235/140‑1 (MQ=255)                                                                          ttAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACa                                                            <  2:8712/140‑1 (MQ=255)                                                                             aGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGc                                                         >  1:752434/1‑140 (MQ=255)                                                                             aGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGc                                                         >  2:675586/1‑140 (MQ=255)                                                                              gTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCt                                                        <  2:229736/140‑1 (MQ=255)                                                                                aaaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                                                      <  1:737673/140‑1 (MQ=255)                                                                                  aaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaag                                                    >  2:681861/1‑140 (MQ=255)                                                                                       cGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAgaagaa                                               >  2:614941/1‑140 (MQ=255)                                                                                          tCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGc                                            <  2:662038/140‑1 (MQ=255)                                                                                              aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGt                                        >  1:366558/1‑140 (MQ=255)                                                                                                    cTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATaa                                  >  2:387039/1‑140 (MQ=255)                                                                                                       aCCGGGCTTACGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATaaaaa                               >  1:84800/1‑140 (MQ=255)                                                                                                                aCGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCg                      <  2:160604/140‑1 (MQ=255)                                                                                                                aCGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCg                      <  1:675586/140‑1 (MQ=255)                                                                                                                aCGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCg                      <  1:614941/140‑1 (MQ=255)                                                                                                                aCGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCg                      <  2:752434/140‑1 (MQ=255)                                                                                                                 cGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGc                     >  1:143853/1‑140 (MQ=255)                                                                                                                 cGTGAGTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGc                     >  1:470598/1‑140 (MQ=255)                                                                                                                      gTTATTCGCTGTATTCTTCATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGt                >  2:22385/1‑140 (MQ=255)                                                                                                                                    tcttcATGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCAtt  >  2:624392/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       TTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATT  >  NC_000913/3726806‑3727076

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