Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,586	C→T	V15I (GTT→ATT)	iclR ←	transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,586	0	C	T	100.0%	204.8 / NA	51	V15I (GTT→ATT)	iclR	transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/41);  total (10/41)

AGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT  >  NC_000913/4223438‑4223720                                                                                                                                                     |                                                                                                                                       aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAAtgg                                                                                                                                      <  1:23810/151‑1 (MQ=255)                gTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGcc                                                                                                                         >  2:291423/1‑149 (MQ=255)                        gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTcgcg                                                                                                               <  1:471747/151‑1 (MQ=255)                         ggTATCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTGCTGCCGCGTTTCGCgg                                                                                                              <  2:492187/151‑1 (MQ=255)                          gTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCggg                                                                                                             <  2:96788/151‑1 (MQ=255)                              cccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCTATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAAt                                                                                                          <  2:75141/150‑1 (MQ=255)                               ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATgg                                                                                                        <  1:139295/151‑1 (MQ=255)                               ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATgg                                                                                                        <  2:649589/151‑1 (MQ=255)                                cGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATggg                                                                                                       <  1:626193/151‑1 (MQ=255)                                      gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAc                                                                                                 <  1:559230/151‑1 (MQ=255)                                         gcgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAt                                                                                              <  2:592348/151‑1 (MQ=255)                                          cgcCAGTTCCGTGAGTGCCACCCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCTTTTCGCGGGAATGGGTGCGACCAt                                                                                              <  2:113203/150‑1 (MQ=255)                                            ccAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc                                                                                           <  1:638724/151‑1 (MQ=255)                                                  ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct                                                                                      <  1:67526/150‑1 (MQ=255)                                                        gTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtttttt                                                                               <  1:172591/151‑1 (MQ=255)                                                         tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGGTTTCTGCCGCGGTTCGCGGGGATGGGTGCGACCATGACAGTCTCtttttt                                                                               >  2:394647/1‑150 (MQ=255)                                                          gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt                                                                             <  1:511700/151‑1 (MQ=255)                                                           ccACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTg                                                                            <  2:175938/151‑1 (MQ=255)                                                            cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGt                                                                           <  1:641754/151‑1 (MQ=255)                                                            cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGt                                                                           <  2:710920/151‑1 (MQ=255)                                                               acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATc                                                                        <  1:86092/151‑1 (MQ=255)                                                                cTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATc                                                                        <  2:190054/150‑1 (MQ=255)                                                                  gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCg                                                                       <  1:665959/149‑1 (MQ=255)                                                                  gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCg                                                                       <  2:761540/149‑1 (MQ=255)                                                                  gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTg                                                                     <  1:184159/151‑1 (MQ=255)                                                                   ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGACACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGTGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCg                                                                       <  2:133920/148‑1 (MQ=255)                                                                    cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGa                                                                   <  1:665008/151‑1 (MQ=255)                                                                        ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAt                                                                <  1:285303/150‑1 (MQ=255)                                                                        ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAt                                                                <  1:733944/150‑1 (MQ=255)                                                                        ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGGTTTTTGCCCCGTTTCGCGGGAAAGGGGGGGACCATGACAGTCTCTTTTTTTTGTATCGTGGgaa                                                                 >  2:546344/1‑149 (MQ=255)                                                                         gATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtca                                                              <  2:655429/151‑1 (MQ=255)                                                                         gATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtca                                                              <  2:560614/151‑1 (MQ=255)                                                                          aTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGAAAtcat                                                             >  2:429643/1‑151 (MQ=255)                                                                          aTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcat                                                             <  1:671541/151‑1 (MQ=255)                                                                             cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcatttt                                                          <  1:645306/151‑1 (MQ=255)                                                                                    ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt                                                   <  2:352704/151‑1 (MQ=255)                                                                                         ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt                                              <  1:337341/151‑1 (MQ=255)                                                                                           aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATtttttttttt                                             >  1:466572/1‑150 (MQ=255)                                                                                            gTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                                           <  2:615331/151‑1 (MQ=255)                                                                                              aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGc                                          <  2:558186/150‑1 (MQ=255)                                                                                                 ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGGTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGGAATCATTTTCATTTTTATTGTTAGCTaa                                       >  2:16173/1‑150 (MQ=255)                                                                                                     ggCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCaa                                  <  2:236741/151‑1 (MQ=255)                                                                                                     ggCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCaa                                  <  2:372428/151‑1 (MQ=255)                                                                                                                aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                       <  1:240214/151‑1 (MQ=255)                                                                                                                aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga                        <  2:161066/150‑1 (MQ=255)                                                                                                                 aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCCTGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTTTTTTTAGCTAATGCAATAGTTactgaa                       >  1:724688/1‑150 (MQ=255)                                                                                                                    aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactga                   <  2:721830/151‑1 (MQ=255)                                                                                                                     cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGCGGCCATGACAGGCTCTTTTTTTTGTATTGGGGGAATCATTTTTATTTTTATTTTTAGCTAATGCAATAGTTtttgaa                       >  2:234387/1‑146 (MQ=255)                                                                                                                       gAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATcc                <  2:632376/151‑1 (MQ=255)                                                                                                                               ccAGTCGCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGtt        <  1:349079/151‑1 (MQ=255)                                                                                                                                     gCTGGTGCGGTGGCAATGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGGGGGACCATGAGAATCTCTTTTTTTTGTATCGTGGAAATCCTTTTTATTTTTTTTGTTAAATAATGCAATAGTTACTGAACTGATCCGGTGAGTTAATGt   >  1:91990/1‑150 (MQ=255)                                                                                                                                       tGGTGCGGTGGCAATGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGTGCGGCCATGACAGGCTCTTTTTTCTGTATCGTGGGAATCATTTTCATTTTTTTTGTTAGCTAATGCAATAATTACTGAAATGATCCGATTAGTTAATGtt  >  2:209877/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                       AGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT  >  NC_000913/4223438‑4223720

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