Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,950,420	G→A	*33* (TAG→TAA)	ilvL →	ilvG operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,950,420	0	G	A	100.0%	206.9 / NA	63	*33* (TAG→TAA)	ilvL	ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/34);  total (29/34)

GGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950287‑3950554                                                                                                                                      |                                                                                                                                       ggCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATc                                                                                                                                   <  2:902567/139‑1 (MQ=255)     ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGc                                                                                                                               >  1:216588/1‑139 (MQ=255)          tCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  2:508567/139‑1 (MQ=255)          tCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  2:866413/139‑1 (MQ=255)                  gATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  <  1:1041401/139‑1 (MQ=255)                  gATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  >  1:390267/1‑139 (MQ=255)                  gATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTACTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  >  1:182444/1‑139 (MQ=255)                            aaGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        <  2:291531/139‑1 (MQ=255)                                 aaaaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  2:762377/139‑1 (MQ=255)                                  aaaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaaa                                                                                                  <  1:437208/139‑1 (MQ=255)                                                      gTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          >  2:1135141/1‑95 (MQ=255)                                                      gTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  1:1135141/95‑1 (MQ=255)                                                       tGATTAGCCTGGTCGTGATTAGCGTGGTGGAGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTCACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGAc                                                                             >  1:157529/1‑139 (MQ=255)                                                        gATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTGACCt                                                                           >  2:1054109/1‑139 (MQ=255)                                                         aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCt                                                                           >  2:1357139/1‑139 (MQ=255)                                                              ccTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGttttttt                                                                                 >  2:320195/1‑128 (MQ=255)                                                              ccTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGttttttt                                                                                 <  1:320195/128‑1 (MQ=255)                                                                  gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCa                                                                                                                                  <  2:666918/75‑1 (MQ=255)                                                                  gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCa                                                                                                                                  >  1:666918/1‑75 (MQ=255)                                                                            gCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCgaggag                                                        <  2:182444/139‑1 (MQ=255)                                                                            gCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTACGGGGGTTTTTTTTGACCTTAAAAACATAACCgaggag                                                        <  2:390267/139‑1 (MQ=255)                                                                             cGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       >  2:1020817/1‑139 (MQ=255)                                                                             cGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAAACCCCCGGAcctttcggtgc                                                                                             <  1:877819/101‑10 (MQ=255)                                                                             cGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAAACCCCCGGAcctttcggtgc                                                                                             >  2:877819/1‑92 (MQ=255)                                                                                ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:1357139/139‑1 (MQ=255)                                                                                  tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                                                                         <  2:1173392/116‑1 (MQ=255)                                                                                  tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                                                                         >  1:1173392/1‑116 (MQ=255)                                                                                           ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                         >  2:559169/1‑139 (MQ=255)                                                                                            tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  >  1:116908/1‑65 (MQ=255)                                                                                            tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  <  2:116908/65‑1 (MQ=255)                                                                                              tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        <  2:389810/73‑1 (MQ=255)                                                                                              tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        >  1:389810/1‑73 (MQ=255)                                                                                                 cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                  <  1:1054109/139‑1 (MQ=255)                                                                                                  actgTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACCAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACa                                                                   >  1:851507/4‑106 (MQ=255)                                                                                                  aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACa                                                                   <  2:851507/106‑1 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  1:555611/50‑1 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          >  2:555611/1‑50 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         <  2:966894/51‑1 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         <  2:541840/51‑1 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         >  1:966894/1‑51 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         >  1:541840/1‑51 (MQ=255)                                                                                                    cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                >  1:1301442/1‑139 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                            <  1:1060212/77‑1 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                            >  2:1060212/1‑77 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  1:393677/1‑139 (MQ=255)                                                                                                        cGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                            <  2:23045/75‑1 (MQ=255)                                                                                                        cGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                            >  1:23045/1‑75 (MQ=255)                                                                                                        cGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                            <  2:1301442/139‑1 (MQ=255)                                                                                                             cTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                       >  2:509047/1‑139 (MQ=255)                                                                                                                cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    <  1:559169/139‑1 (MQ=255)                                                                                                                     ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGc                                                                                                         <  1:144216/49‑1 (MQ=255)                                                                                                                     ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGc                                                                                                         >  2:144216/1‑49 (MQ=255)                                                                                                                     ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                            >  1:23214/1‑62 (MQ=255)                                                                                                                     ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                            <  2:23214/62‑1 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                                                                         <  2:131112/76‑1 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                                                                         >  1:131112/1‑76 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa          >  2:1001292/1‑139 (MQ=255)                                                                                                                            gAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg        >  1:885093/1‑139 (MQ=255)                                                                                                                            gAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg        <  1:1324756/139‑1 (MQ=255)                                                                                                                               aGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg     <  2:216588/139‑1 (MQ=255)                                                                                                                                 gCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa   <  1:1020817/139‑1 (MQ=255)                                                                                                                                  cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  1:509047/139‑1 (MQ=255)                                                                                                                                  cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  2:885093/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       GGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950287‑3950554

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