Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,908,423	C→T	Q42Q (CAG→CAA)	yebO ←	putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,908,423	0	C	T	100.0%	188.3 / NA	51	Q42Q (CAG→CAA)	yebO	putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/24);  total (27/24)

CCAGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATCTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTA  >  NC_000913/1908275‑1908565                                                                                                                                                     |                                                                                                                                               ccAGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTg                                                                                                                                              >  2:442870/1‑151 (MQ=255) ccAGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTg                                                                                                                                              >  1:597061/1‑151 (MQ=255)  cAGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGt                                                                                                                                             <  2:347409/151‑1 (MQ=255)    gACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGtt                                                                                                                                            >  2:56141/1‑150 (MQ=255)      cGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGt                                                                                                                                         <  2:243170/151‑1 (MQ=255)        aataatATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGtt                                                                                                                                        >  1:157464/1‑150 (MQ=255)          taataTCTTCGTATTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTAGTACAAACGACGTAACAGTTCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGCTCAATTTGTTCGTTAg                                                                                                                                      >  1:369881/1‑150 (MQ=255)           aataTCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTc                                                                                                                                    >  2:412507/1‑151 (MQ=255)           aataTCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTc                                                                                                                                    >  1:305644/1‑151 (MQ=255)              atCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCg                                                                                                                                  <  1:554978/150‑1 (MQ=255)                      cTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGc                                                                                                                         >  1:452983/1‑151 (MQ=255)                        tCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCa                                                                                                                        <  1:813885/150‑1 (MQ=255)                              gATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTg                                                                                                                  >  1:427183/1‑150 (MQ=255)                              gATTTTTGACCCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGa                                                                                                                 <  2:416442/151‑1 (MQ=255)                                       ctctCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAAcc                                                                                                        >  2:624341/1‑151 (MQ=255)                                        tctcAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAAcc                                                                                                        >  2:681852/1‑150 (MQ=255)                                            aaCGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACcaca                                                                                                     >  2:66008/1‑149 (MQ=255)                                            aaCGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACcaca                                                                                                     >  1:467553/1‑149 (MQ=255)                                                     tttATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGccc                                                                                           <  2:466703/150‑1 (MQ=255)                                                        aTCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGAt                                                                                        >  1:139229/1‑150 (MQ=255)                                                         tCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGAt                                                                                        >  2:688212/1‑149 (MQ=255)                                                               ttCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATaagaag                                                                                  <  1:535514/149‑1 (MQ=255)                                                               ttCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATaagaag                                                                                  <  2:597100/149‑1 (MQ=255)                                                                tctcTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAac                                                                               >  2:301824/1‑151 (MQ=255)                                                                 ctctGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGa                                                                                                  >  1:155070/1‑131 (MQ=255)                                                                 ctctGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGa                                                                                                  <  2:155093/131‑1 (MQ=255)                                                                 ctctGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGAATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAgaagaa                                                                                >  1:27268/1‑149 (MQ=255)                                                                   ctGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAaccac                                                                            >  1:150111/1‑151 (MQ=255)                                                                     ggTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAaccacc                                                                           >  2:467153/1‑150 (MQ=255)                                                                        tCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCg                                                                       >  1:199958/1‑151 (MQ=255)                                                                          gTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGAt                                                                     >  2:162771/1‑151 (MQ=255)                                                                              gCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACa                                                                  <  1:373871/150‑1 (MQ=255)                                                                                  cGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCa                                                              >  1:618606/1‑150 (MQ=255)                                                                                   gCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCaaa                                                            <  1:143944/151‑1 (MQ=255)                                                                                    cacaAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAg                                                           <  2:226241/151‑1 (MQ=255)                                                                                    cacaAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAg                                                           <  1:442836/151‑1 (MQ=255)                                                                                       aaaCGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACg                                                        <  1:663189/151‑1 (MQ=255)                                                                                         acgacgTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCa                                                      <  2:34663/151‑1 (MQ=255)                                                                                         acgacgTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCa                                                      <  2:417681/151‑1 (MQ=255)                                                                                            acgTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATg                                                   <  1:618449/151‑1 (MQ=255)                                                                                             cgTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGt                                                  >  2:178268/1‑151 (MQ=255)                                                                                                aaCAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTtc                                                <  1:294027/150‑1 (MQ=255)                                                                                                 aCAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTtcat                                              <  1:822268/151‑1 (MQ=255)                                                                                                             gACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAAtt                                  >  1:519339/1‑151 (MQ=255)                                                                                                                         gATCCAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGt                       <  1:178245/150‑1 (MQ=255)                                                                                                                            ccAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCAt                   <  2:740938/151‑1 (MQ=255)                                                                                                                            ccAGCAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCAt                   <  1:645043/151‑1 (MQ=255)                                                                                                                                cAACGCCTCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTca               <  1:258775/151‑1 (MQ=255)                                                                                                                                       tCAAGCAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACctc        >  1:112567/1‑151 (MQ=255)                                                                                                                                            cAGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGt   >  1:59415/1‑151 (MQ=255)                                                                                                                                             aGTTCAATTTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTa  <  2:666862/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CCAGACGAATAATATCTTCGTCTTCAACCGATTTTTGACTCTCAACGGTCTTTTTATCTGCTTTCTCTGGTTCGTTTGCTTCGCACAAACGACGTAACAGTGCATTTTGACGTTTTTGCTGATCCAGCAACGCCTCAAGCAGTTCAATCTGTTCGTTAGTCCGTGAACTGGCACGATTGATAAAAAACCACAAGATGAGCCCGATAAGAAGAACCACCACCGATACAACCAAAGACGCAATGTTCATCACGCCTGAATTTACAACTTCGTTCATTTCACCACCTCTGAGTA  >  NC_000913/1908275‑1908565

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