Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	916,431	A→G	intergenic (‑384/‑42)	aqpZ ← / → ybjD	water channel AqpZ/DUF2813 domain‑containing protein YbjD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	916,431	0	A	G	100.0%	135.7 / NA	38	intergenic (‑384/‑42)	aqpZ/ybjD	water channel AqpZ/DUF2813 domain‑containing protein YbjD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/18);  total (20/18)

GCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAACG  >  NC_000913/916296‑916566                                                                                                                                        |                                                                                                                                        gCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTAtgctgc                                                                                                                                     <  2:827488/140‑1 (MQ=255)      ttattTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCg                                                                                                                                >  2:482851/1‑140 (MQ=255)           tGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCa                                                                                                                                   <  1:692311/132‑1 (MQ=255)           tGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCa                                                                                                                                   >  2:692311/1‑132 (MQ=255)               tttGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTTATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATa                                                                                                                       >  2:537500/1‑140 (MQ=255)                ttgttgTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGcc                                                                                                                                    >  1:836162/1‑126 (MQ=255)                ttgttgTTTTTGCAGACTTATCAGCAAGAGGGAGTATAAAGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGcc                                                                                                                                    <  2:836162/126‑1 (MQ=255)                  gttgttTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGt                                                                                                                    <  2:234582/140‑1 (MQ=255)                      tttttGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTCCCCGTTATGCTGCCAGCGTAAAGGATAAGTcaca                                                                                                                <  2:852105/140‑1 (MQ=255)                             aGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCt                                                                                                         >  2:873358/1‑140 (MQ=255)                                              ggAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCt                                                                                        >  2:780833/1‑140 (MQ=255)                                               gAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGa                                                                                                 >  1:203007/1‑130 (MQ=255)                                               gAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGa                                                                                                 <  2:203007/130‑1 (MQ=255)                                                  tataACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAg                                                                                    >  1:820309/1‑140 (MQ=255)                                                  tataACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAg                                                                                    >  2:183081/1‑140 (MQ=255)                                                  tataACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAg                                                                                    >  1:776444/1‑140 (MQ=255)                                                    taACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAgcg                                                                                  <  2:820309/140‑1 (MQ=255)                                                            ttattCGCTCATTTTTCAGACATTTG‑CATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAAt                                                                         >  1:124672/1‑140 (MQ=255)                                                                         tttCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTcg                                                             <  1:873358/140‑1 (MQ=255)                                                                                cATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATc                                                      <  2:621980/140‑1 (MQ=255)                                                                                  tttGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCaa                                                    >  1:96011/1‑140 (MQ=255)                                                                                        aTGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGttt                                              <  1:183081/140‑1 (MQ=255)                                                                                                aTGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGa                                      >  1:723859/1‑140 (MQ=255)                                                                                                     aTGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTg                                 >  2:92300/1‑140 (MQ=255)                                                                                                     aTGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTg                                 >  1:476269/1‑140 (MQ=255)                                                                                                     aTGTCATCACGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTg                                 >  2:694706/1‑140 (MQ=255)                                                                                                        tcatcaCGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGaa                              <  1:694706/140‑1 (MQ=255)                                                                                                          atcaCGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGaaca                            <  2:339852/140‑1 (MQ=255)                                                                                                          atcaCGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGaaca                            <  1:537500/140‑1 (MQ=255)                                                                                                          atcaCGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGaaca                            <  1:815126/140‑1 (MQ=255)                                                                                                              cGTATTAGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGaacaaaac                        <  1:92300/140‑1 (MQ=255)                                                                                                                    aGCAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTc                  <  2:96011/140‑1 (MQ=255)                                                                                                                      cAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGt                   >  1:869342/1‑137 (MQ=255)                                                                                                                      cAAGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAAGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGt                   <  2:869342/137‑1 (MQ=255)                                                                                                                        aGGCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGa              >  2:199389/1‑140 (MQ=255)                                                                                                                         ggCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGa                               >  1:63845/1‑122 (MQ=255)                                                                                                                         ggCCTTTCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGa                               <  2:63845/122‑1 (MQ=255)                                                                                                                             tttCCCGTTATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTggg         <  2:476269/140‑1 (MQ=255)                                                                                                                                    ttATGCTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAAcg  >  2:236856/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        GCTAATTATTTGATTTTGTTGTTTTTGCAGACTTATCAGCAAGAGGGAGTATAACGCGATTATTCGCTCATTTTTCAGACATTTGCCATGCTTAAATGTGATGTCATCACGTATTAGCAAGGCCTTTCCCGTTATACTGCCAGCGTAAAGGATAAGTCACATATTTCTGGAGGGGATATGATTCTTGAGCGCGTTGAAATTGTGGGTTTTCGCGGTATCAACCGTTTGTCGTTGATGCTGGAACAAAACAACGTCCTGATTGGGGAGAACG  >  NC_000913/916296‑916566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑