Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,994,150	G→T	D101Y (GAT→TAT)	yzcX →	protein YzcX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,994,150	0	G	T	100.0%	165.2 / NA	47	D101Y (GAT→TAT)	yzcX	protein YzcX
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/31);  total (16/31)

AGCCAGAGTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCAT  >  NC_000913/3994016‑3994277                                                                                                                                       |                                                                                                                                aGCCAGAGTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCa                                                                                                                            <  2:403630/140‑1 (MQ=255) aGCCAGAGTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCa                                                                                                                            >  2:651715/1‑140 (MQ=255)      gagTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGt                                                                                                                       <  2:437136/140‑1 (MQ=255)       agTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTc                                                                                                                      >  1:727031/1‑140 (MQ=255)        gTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCa                                                                                                                     <  1:595024/140‑1 (MQ=255)        gTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCa                                                                                                                     >  1:89738/1‑140 (MQ=255)         ttGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCAc                                                                                                                    <  2:384048/140‑1 (MQ=255)         ttGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCAc                                                                                                                    <  2:89738/140‑1 (MQ=255)         ttGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCAc                                                                                                                    <  2:271924/140‑1 (MQ=255)         ttGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCAc                                                                                                                    <  2:184889/140‑1 (MQ=255)         ttGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCAc                                                                                                                    <  2:695473/140‑1 (MQ=255)         ttGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCAc                                                                                                                    <  2:110490/140‑1 (MQ=255)             tCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGa                                                                                                                <  2:766807/140‑1 (MQ=255)             tCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGa                                                                                                                >  2:651793/1‑140 (MQ=255)                gCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGaata                                                                                                             <  2:58682/140‑1 (MQ=255)                 ccAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGaataa                                                                                                            >  2:93953/1‑140 (MQ=255)                    ggCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAatga                                                                                                         <  1:33199/140‑1 (MQ=255)                              ttCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCtgtg                                                                                               <  2:639564/140‑1 (MQ=255)                                   tGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCg                                                                                          <  1:219090/140‑1 (MQ=255)                                   tGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCg                                                                                          >  2:325327/1‑140 (MQ=255)                                    gTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGa                                                                                         <  1:711571/140‑1 (MQ=255)                                     tCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGaa                                                                                        >  1:294667/1‑140 (MQ=255)                                           aTTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCTCTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGAttt                                                                                  <  2:117413/140‑1 (MQ=255)                                             tgttgtATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGt                                                                                <  2:148428/140‑1 (MQ=255)                                               ttgtATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGtct                                                                              >  1:50939/1‑140 (MQ=255)                                                     ccTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCAt                                                                        <  2:159614/140‑1 (MQ=255)                                                       tGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTa                                                                      >  1:601475/1‑140 (MQ=255)                                                         tttttAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGtct                                                                              >  1:93699/1‑130 (MQ=255)                                                         tttttAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGtct                                                                              <  2:93699/130‑1 (MQ=255)                                                          ttttAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAggg                                                                   <  1:93953/140‑1 (MQ=255)                                                            ttAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGgcg                                                                 >  1:743396/1‑140 (MQ=255)                                                                                 gCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTccc                                            <  2:602175/140‑1 (MQ=255)                                                                                       tGGGGCCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTc                                      >  2:216286/1‑140 (MQ=255)                                                                                        ggggTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCt                                     <  2:381127/140‑1 (MQ=255)                                                                                           gTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGt                                  <  2:601475/140‑1 (MQ=255)                                                                                            tCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTa                                 <  2:294667/140‑1 (MQ=255)                                                                                                     cAAAGTTTGCGCCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATg                        <  1:225225/140‑1 (MQ=255)                                                                                                     cAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATg                        <  1:505639/140‑1 (MQ=255)                                                                                                       aaGTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGAt                      <  2:727031/140‑1 (MQ=255)                                                                                                         gTTTGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCa                    <  1:651715/140‑1 (MQ=255)                                                                                                            tGCACCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGt                 >  2:597425/1‑140 (MQ=255)                                                                                                               aCCTCTGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTa              <  2:50939/140‑1 (MQ=255)                                                                                                                    tGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcgag         >  1:729900/1‑140 (MQ=255)                                                                                                                    tGCGGCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcgag         >  1:517462/1‑140 (MQ=255)                                                                                                                       ggCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcgagcga      <  2:167230/140‑1 (MQ=255)                                                                                                                        gCTGCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACcgagcgag     <  2:687568/140‑1 (MQ=255)                                                                                                                           gCGTTCCGGCACTATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCAt  >  2:619073/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                AGCCAGAGTTGATCAGCCAGGCGATGAAATTCACTGTCGTTCATTGTTGTATCCTGTTTTTAAGTGATGGCGGCAGTATAGCGGCATGGGGTCAGGGCTTCAAAGTTTGCACCTCTGCGGCTGCGTTCCGGCACGATTCATCCGTCACCGGAATAATGATGTCTCTGTGTAGCGAAAGATTTGTCTCTTCATTAGGGCGCAGTTACACCACGTCTTTCCCTGTTTCTGGTAAACATTATGATCAGGTTTACCGAGCGAGCAT  >  NC_000913/3994016‑3994277

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: ‑