Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,993,856	G→A	*73* (TAG→TAA)	ygeI →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,993,856	0	G	A	100.0%	210.5 / NA	66	*73* (TAG→TAA)	ygeI	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/33);  total (33/33)

ACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993726‑2993987                                                                                                                                   |                                                                                                                                    acaTCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAATTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             >  1:810403/1‑139 (MQ=255)               gCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATa                                                                                                               <  2:599669/139‑1 (MQ=255)                   cTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaaa                                                                                                                        <  2:463018/126‑1 (MQ=255)                   cTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaaa                                                                                                                        >  1:463018/1‑126 (MQ=255)                     gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTCTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         <  2:810403/139‑1 (MQ=255)                      aTACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:1332779/139‑1 (MQ=255)                        aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgag                                                                                                          >  1:1012213/1‑135 (MQ=255)                        aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgag                                                                                                          <  2:1012213/135‑1 (MQ=255)                        aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  1:736994/139‑1 (MQ=255)                           gATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                   >  2:20679/1‑139 (MQ=255)                             tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             >  2:689464/1‑111 (MQ=255)                             tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             <  1:689464/111‑1 (MQ=255)                             tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  2:66556/134‑1 (MQ=255)                             tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      >  1:66556/1‑134 (MQ=255)                                gaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCACACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                              >  1:826941/2‑139 (MQ=255)                                aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  1:299399/131‑1 (MQ=255)                                aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      >  2:299399/1‑131 (MQ=255)                                           ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             <  2:838094/97‑1 (MQ=255)                                           ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             >  1:838094/1‑97 (MQ=255)                                               aGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCt                                                                               <  1:929669/139‑1 (MQ=255)                                                gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt                                                                              <  1:310351/139‑1 (MQ=255)                                                     aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         <  2:826941/139‑1 (MQ=255)                                                     aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         >  1:700134/1‑139 (MQ=255)                                                           gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         >  2:1244194/1‑101 (MQ=255)                                                           gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         <  1:1244194/101‑1 (MQ=255)                                                                 gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacac                                                             >  2:1324573/1‑139 (MQ=255)                                                                      ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                        >  1:484500/1‑139 (MQ=255)                                                                      ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                        >  2:680453/1‑139 (MQ=255)                                                                                      ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         <  1:594435/74‑1 (MQ=255)                                                                                      ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         >  2:594435/1‑74 (MQ=255)                                                                                       tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                       >  2:932615/1‑139 (MQ=255)                                                                                         agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  1:1308624/74‑1 (MQ=255)                                                                                         agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      >  2:1308624/1‑74 (MQ=255)                                                                                         agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                          >  1:373102/1‑86 (MQ=255)                                                                                         agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                          <  2:373102/86‑1 (MQ=255)                                                                                          gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTca                                                                <  2:673886/111‑1 (MQ=255)                                                                                          gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTca                                                                >  1:673886/1‑111 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:226670/1‑67 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:547852/1‑67 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:547852/67‑1 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  1:242778/1‑67 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  2:242778/67‑1 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:226670/67‑1 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   <  1:5791/88‑1 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   >  2:5791/1‑88 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                               <  2:1288142/108‑1 (MQ=255)                                                                                              aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                               >  1:1288142/1‑108 (MQ=255)                                                                                               tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaa                               <  1:680453/139‑1 (MQ=255)                                                                                                aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                              >  1:803347/1‑139 (MQ=255)                                                                                                 aTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                             >  1:738082/1‑139 (MQ=255)                                                                                                        ctcACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                      <  2:738082/139‑1 (MQ=255)                                                                                                            ccTTACTATAAAATTGAGTGGTAAACGATTGACCaaata                                                                                                                      >  2:481493/1‑39 (MQ=255)                                                                                                            ccTTACTATAAAATTGAGTGGTAAACGATTGACCaaata                                                                                                                      <  1:481493/39‑1 (MQ=255)                                                                                                            ccTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:786923/53‑1 (MQ=255)                                                                                                            ccTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:786923/1‑53 (MQ=255)                                                                                                              ttACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCtt                <  1:932615/139‑1 (MQ=255)                                                                                                               tACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc               <  2:1216630/139‑1 (MQ=255)                                                                                                                aCTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                >  1:696963/1‑121 (MQ=255)                                                                                                                aCTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                <  2:696963/121‑1 (MQ=255)                                                                                                                aCTATAAAATTGAGTGGTAAACGATTGACCAAATAAATAATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTAAATCCTTATAAATGAAAAAATATGTGCTCACGAAAAGTTAAATTCTTAACCATCTTCa              >  2:1006597/1‑139 (MQ=255)                                                                                                                    taAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt          <  1:20679/139‑1 (MQ=255)                                                                                                                     aaaaTTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                        <  1:446404/60‑1 (MQ=255)                                                                                                                     aaaaTTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                        >  2:446404/1‑60 (MQ=255)                                                                                                                            aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAcc                                        <  2:735516/101‑1 (MQ=255)                                                                                                                            aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAcc                                        >  1:735516/1‑101 (MQ=255)                                                                                                                            aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat  >  1:1229184/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                    ACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993726‑2993987

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