Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	915,343	C→T	intergenic (‑486/+9)	ybjE ← / ← aqpZ	putative transporter/aquaporin Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	915,343	0	C	T	100.0%	125.3 / NA	42	intergenic (‑486/+9)	ybjE/aqpZ	putative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/23);  total (19/23)

GGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCG  >  NC_000913_3_pae_tpiA/915203‑915488                                                                                                                                             |                                                                                                                                                  ggcggcggcGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGttaa                                                                                                                                           <  2:178424/149‑1 (MQ=255)  ggggcggccAATCGCCGCCGCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGCTAATACGCCAGGCCTTTTGCGttaat                                                                                                                                          >  2:236318/3‑149 (MQ=255)          gATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTCGCTCATGTAAGGAGCGCAACAAGATGATGAATCAACATATAACATGTACTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGttaattatttccg                                                                                                                                  >  2:53414/1‑143 (MQ=38)                ccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTc                                                                                                                            <  1:800741/149‑1 (MQ=255)                  gccCCCTGCGCGCCCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGCAAATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTAAATAAATAACGCTTTTCca                                                                                                                          >  2:415833/1‑149 (MQ=255)                    ccccTGCGCGCGCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAGGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACAACATCTTAAAAAAAGGCCTGACATCACGCCAGGCCTTTTGCGTCGAATAATCACCCTTTTCcagc                                                                                                                        >  2:849504/1‑149 (MQ=255)                                 cTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGCGCTCAACAAGATGATGGATCAACAAATAACAACCTCGTAAAAAAAGACCTAACATTACGCCAGGCCTTTTACGATAATTAATCACGCTTTTCAACCAGGGTCCGGaaaa                                                                                                           >  2:397931/1‑149 (MQ=255)                                 cTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTcc                                                                                                                           <  2:906362/133‑1 (MQ=255)                                 cTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTcc                                                                                                                           >  1:906362/1‑133 (MQ=255)                                  gtgggcCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATCAATCAACATATAACAGCATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTCCCAGCAGGGTCCGGTAAAt                                                                                                          <  2:791005/145‑1 (MQ=255)                                      gcaaGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGa                                                                                                      <  1:615514/146‑1 (MQ=255)                                              aaaTATTTTGCTTAGGTCATCAATGGTGCGCAACCAGATCATGAATCAACATATATCAACATCTTAAAAAAAGGCCTGAGAATACGCCAGGCCTTTTGCGTTAATTAATCACGCTCTTCCCGCCGGCTCCGTTCAATCAGACTACCGAt                                                                                              <  2:768355/148‑1 (MQ=255)                                              aaaTAGGTTGCTTAGCTGATGAAAGTAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGAt                                                                                              <  1:846263/148‑1 (MQ=255)                                                 tatGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAAt                                                                                           >  1:899840/1‑149 (MQ=255)                                                    gTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTCATCACGCTGTTCCAGCAGGGTCCGGTAAATCAGACCACGGATAATgcg                                                                                        >  1:429374/1‑148 (MQ=255)                                                      tGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgc                                                                                      <  2:568685/149‑1 (MQ=255)                                                       gCTTAGTGCATGAAGGGTGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgcc                                                                                     <  1:66346/149‑1 (MQ=255)                                                        cTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGccgccg                                                                                    >  1:699384/1‑149 (MQ=255)                                                          tAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGAc                                                                                  >  1:679064/1‑149 (MQ=255)                                                             cTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAACTAACACCACCGATAATGCCGCCGACAAt                                                                               >  1:663258/1‑149 (MQ=255)                                                              tCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAtt                                                                              <  2:699384/149‑1 (MQ=255)                                                              tCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAtt                                                                              <  2:924882/149‑1 (MQ=255)                                                                aTGAAAGGCGCGCAACAAGATAATGGATCAACATATAACAGCATCTTGAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTgg                                                                            >  2:557067/1‑149 (MQ=255)                                                                           gcAACAAGATGATGAACCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCaga                                                                 <  1:331339/149‑1 (MQ=255)                                                                               caaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAgaagaa                                                             <  1:438100/149‑1 (MQ=255)                                                                                aaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAAc                                                            <  2:844437/149‑1 (MQ=255)                                                                                           tCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTc                                                 >  1:525141/1‑149 (MQ=255)                                                                                               caaataACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAAt                                             <  1:810598/146‑1 (MQ=255)                                                                                                    aacaacATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCa                                        >  1:845232/1‑149 (MQ=255)                                                                                                      caacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgc                                      >  1:835236/1‑149 (MQ=255)                                                                                                              tAAAAAAAGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTgg                              <  1:396265/149‑1 (MQ=255)                                                                                                               aaaaaaaGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCGGCCCTGGa                             >  1:907265/1‑149 (MQ=255)                                                                                                                    aaGGCCTGAGATTACGCCAGGTCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATa                        <  2:22005/149‑1 (MQ=255)                                                                                                                    aaGGCCTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATa                        <  1:656230/149‑1 (MQ=255)                                                                                                                         cTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATa                                                                                             >  1:782503/1‑75 (MQ=255)                                                                                                                         cTGACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATa                                                                                             <  2:782503/75‑1 (MQ=255)                                                                                                                           gACATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGACCCACAGTTGTTCTAATGCCCACCCCCCCTCGAAAAAAACAAccc                 >  2:369277/1‑148 (MQ=255)                                                                                                                            aCATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCt                                >  1:752997/1‑133 (MQ=255)                                                                                                                            aCATTACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCt                                <  2:752997/133‑1 (MQ=255)                                                                                                                                tACGCCAGGCCTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGtg            <  1:473884/149‑1 (MQ=255)                                                                                                                                      aGGCCTTTGGCGGTAATTAAGCACGCTTGTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTgcgc      <  2:834411/149‑1 (MQ=255)                                                                                                                                          cTTTTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGATCACCGATAAGGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCg  <  1:415833/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  GGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCG  >  NC_000913_3_pae_tpiA/915203‑915488

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