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

GGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950287‑3950554                                                                                                                                      |                                                                                                                                       ggCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATc                                                                                                                                   <  2:1004747/139‑1 (MQ=255)    aTTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCATTCTACGAGGGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCAACGGGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAg                                                                                                                                <  1:409552/139‑1 (MQ=255)       ccTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGCGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCt                                                                                                                             >  2:906564/1‑139 (MQ=255)                  gATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGa                                                                                                                       <  2:353194/134‑1 (MQ=255)                  gATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGa                                                                                                                       >  1:353194/1‑134 (MQ=255)                                 aaaaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGGGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCAATTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  1:906564/139‑1 (MQ=255)                                   aaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATc                                                                                                                                   >  2:440262/1‑105 (MQ=255)                                   aaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATc                                                                                                                                   <  1:440262/105‑1 (MQ=255)                                                   cGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGttttttt                                                                                 >  1:250264/1‑139 (MQ=255)                                                      gTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  2:978761/95‑1 (MQ=255)                                                      gTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          >  1:978761/1‑95 (MQ=255)                                                         aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  <  1:990847/100‑1 (MQ=255)                                                         aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  >  2:990847/1‑100 (MQ=255)                                                         aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCt                                                                           >  2:214063/1‑139 (MQ=255)                                                           tAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                                                                         <  1:214063/139‑1 (MQ=255)                                                               cTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa                                                                     >  1:447758/1‑139 (MQ=255)                                                                tGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAAc                                                                    >  1:443320/1‑139 (MQ=255)                                                                 ggTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACa                                                                   <  1:611122/139‑1 (MQ=255)                                                                  gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    <  1:50714/73‑1 (MQ=255)                                                                  gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    >  2:50714/1‑73 (MQ=255)                                                                                ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:447758/139‑1 (MQ=255)                                                                                ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:293231/139‑1 (MQ=255)                                                                                ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:33778/139‑1 (MQ=255)                                                                                           ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  2:563613/58‑1 (MQ=255)                                                                                           ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          >  1:563613/1‑58 (MQ=255)                                                                                                 cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         <  2:343474/53‑1 (MQ=255)                                                                                                 cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         >  1:343474/1‑53 (MQ=255)                                                                                                 cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                   >  1:734546/1‑139 (MQ=255)                                                                                                 cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                   >  2:185150/1‑139 (MQ=255)                                                                                                 cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                   <  2:250264/139‑1 (MQ=255)                                                                                                  aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAAc                                                                                                                     >  2:288553/1‑56 (MQ=255)                                                                                                  aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAAc                                                                                                                     <  1:288553/56‑1 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  <  2:38489/58‑1 (MQ=255)                                                                                                   ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  >  1:38489/1‑58 (MQ=255)                                                                                                    cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:706129/119‑1 (MQ=255)                                                                                                    cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  2:706129/1‑119 (MQ=255)                                                                                                    cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                >  1:274197/1‑139 (MQ=255)                                                                                                     gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:325321/1‑118 (MQ=255)                                                                                                     gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:325321/118‑1 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       <  1:566867/114‑1 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       >  2:566867/1‑114 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              <  2:734546/139‑1 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  2:369239/1‑139 (MQ=255)                                                                                                      tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  2:271149/1‑139 (MQ=255)                                                                                                           ggCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                         >  2:256946/1‑139 (MQ=255)                                                                                                             cTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                       >  1:882332/1‑139 (MQ=255)                                                                                                                cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:218468/107‑1 (MQ=255)                                                                                                                cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:218468/1‑107 (MQ=255)                                                                                                                cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    <  2:443320/139‑1 (MQ=255)                                                                                                                       aCGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa             <  2:882332/139‑1 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          >  2:108507/1‑75 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          <  1:108507/75‑1 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:54863/97‑1 (MQ=255)                                                                                                                          aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:54863/1‑97 (MQ=255)                                                                                                                                 gCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:822744/90‑1 (MQ=255)                                                                                                                                 gCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  2:822744/1‑90 (MQ=255)                                                                                                                                  cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  1:369239/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       GGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950287‑3950554

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