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%	213.7 / NA	65	*33* (TAG→TAA)	ilvL	ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/33);  total (32/33)

TAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950285‑3950554                                                                                                                                        |                                                                                                                                        tAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAaaga                                                                                                                                      >  2:657235/1‑139 (MQ=255)   ggCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATc                                                                                                                                    <  2:546857/139‑1 (MQ=255)         ccTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCt                                                                                                                              >  1:421346/1‑139 (MQ=255)             cGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                          <  2:644880/139‑1 (MQ=255)                                         aCAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                     <  2:78267/100‑1 (MQ=255)                                         aCAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                     >  1:78267/1‑100 (MQ=255)                                                       aGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                    <  1:144869/119‑1 (MQ=255)                                                       aGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                    >  2:144869/1‑119 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                           >  1:134104/1‑92 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                           <  2:134104/92‑1 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                   <  2:453562/100‑1 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                   >  1:453562/1‑100 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTg                                                                               <  2:983846/135‑1 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTg                                                                               >  1:983846/1‑135 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCt                                                                           >  2:326677/1‑139 (MQ=255)                                                           aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCt                                                                           >  1:44597/1‑139 (MQ=255)                                                            ttAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCtt                                                                          <  2:421346/139‑1 (MQ=255)                                                                cttGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGAc                                                                             <  1:673776/130‑1 (MQ=255)                                                                ccTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTaaaa                                                                      >  1:1117260/1‑139 (MQ=255)                                                                ccTGGTCGTGATTAGCGTGGTGGTGATTATTAACCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGAc                                                                             >  2:673776/1‑132 (MQ=255)                                                                 cTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTaaaaa                                                                     <  2:1117260/139‑1 (MQ=255)                                                                        tGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAAcc                                                              <  2:35249/139‑1 (MQ=255)                                                                          aTTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCg                                                                                                           >  1:786167/1‑93 (MQ=255)                                                                          aTTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCg                                                                                                           <  2:786167/93‑1 (MQ=255)                                                                              gCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCgaggag                                                        <  2:313933/139‑1 (MQ=255)                                                                                   gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                           >  1:792350/1‑68 (MQ=255)                                                                                   gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                           <  2:792350/68‑1 (MQ=255)                                                                                   gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCgaggag                                                        <  2:468089/134‑1 (MQ=255)                                                                                   gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCgaggag                                                        >  1:468089/1‑134 (MQ=255)                                                                                              tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:76429/127‑1 (MQ=255)                                                                                              tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  2:76429/1‑127 (MQ=255)                                                                                                   cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                   >  2:543335/1‑139 (MQ=255)                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                          >  1:1131920/1‑51 (MQ=255)                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                          <  2:1131920/51‑1 (MQ=255)                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      >  1:621400/1‑118 (MQ=255)                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      <  2:621400/118‑1 (MQ=255)                                                                                                     ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                 >  2:224660/1‑139 (MQ=255)                                                                                                      cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                             <  2:723218/79‑1 (MQ=255)                                                                                                      cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                                                                             >  1:723218/1‑79 (MQ=255)                                                                                                      cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      >  1:845278/1‑117 (MQ=255)                                                                                                      cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      <  2:845278/117‑1 (MQ=255)                                                                                                       gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                  >  1:788655/1‑120 (MQ=255)                                                                                                       gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                  <  2:788655/120‑1 (MQ=255)                                                                                                       gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                               >  2:396215/1‑139 (MQ=255)                                                                                                        tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  2:857404/1‑139 (MQ=255)                                                                                                        tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  1:781003/1‑139 (MQ=255)                                                                                                          cGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑‑TTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:836119/1‑114 (MQ=255)                                                                                                          cGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑‑TTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:836119/114‑1 (MQ=255)                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTa                                                                         >  2:205282/1‑86 (MQ=255)                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTa                                                                         <  1:205282/86‑1 (MQ=255)                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAAcc                                                              <  1:633229/97‑1 (MQ=255)                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAAcc                                                              >  2:633229/1‑97 (MQ=255)                                                                                                                  cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    <  2:44597/139‑1 (MQ=255)                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                  <  2:790917/105‑1 (MQ=255)                                                                                                                       ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                  >  1:790917/1‑105 (MQ=255)                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGa                                                                                                     >  1:395462/1‑49 (MQ=255)                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGa                                                                                                     <  2:395462/49‑1 (MQ=255)                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTg                                                                               >  2:2298/1‑70 (MQ=255)                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTg                                                                               <  1:2298/70‑1 (MQ=255)                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:533964/97‑1 (MQ=255)                                                                                                                            aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  2:533964/1‑97 (MQ=255)                                                                                                                              gAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg        <  1:543335/139‑1 (MQ=255)                                                                                                                               aaaGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg       <  1:224660/139‑1 (MQ=255)                                                                                                                                  ggCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑‑TTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       <  1:197144/87‑1 (MQ=255)                                                                                                                                  ggCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑‑TTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       >  2:197144/1‑87 (MQ=255)                                                                                                                                    cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  1:857404/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950285‑3950554

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