Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,117,133	A→G	intergenic (‑5/+480)	b2974 ← / ← yghK	putative endoglucanase/putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,117,133	0	A	G	100.0%	243.0 / NA	78	intergenic (‑5/+480)	b2974/yghK	putative endoglucanase/putative permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/34);  total (44/34)

CCACCGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATAACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCATTCTTA  >  NC_000913/3117046‑3117221                                                                                        |                                                                                         ccaccGCCATCACACCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTGAAATTTCTTATTCATTAATGACGCa                                                                                     <  2:456574/93‑1 (MQ=255)   cccGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAg                                                                                   <  2:961339/92‑1 (MQ=255)    ccGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGt                                                                                  <  2:1320961/93‑1 (MQ=255)    ccGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGt                                                                                  >  2:2547489/1‑93 (MQ=255)         aTCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACaa                                                                             >  2:183648/1‑93 (MQ=255)         aTCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACaa                                                                             <  1:2036486/93‑1 (MQ=255)          tCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACaaa                                                                            <  1:1999333/93‑1 (MQ=255)           cacaACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACaaaa                                                                           <  1:3802490/93‑1 (MQ=255)            acaACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAAc                                                                          <  2:945984/93‑1 (MQ=255)              aaCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACaa                                                                        <  2:1403909/93‑1 (MQ=255)              aaCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACaa                                                                        <  2:4353927/93‑1 (MQ=255)              aaCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACaa                                                                        <  2:4479757/93‑1 (MQ=255)               aCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAAt                                                                       <  1:2565476/93‑1 (MQ=255)               aCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAAt                                                                       <  2:4308608/93‑1 (MQ=255)                ccGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATg                                                                      <  1:2377678/93‑1 (MQ=255)                ccGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATg                                                                      >  1:606288/1‑93 (MQ=255)                 cGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGt                                                                     <  1:2183439/93‑1 (MQ=255)                           ggTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTc                                                           <  2:4528741/93‑1 (MQ=255)                            gTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCg                                                          >  1:2033521/1‑93 (MQ=255)                            gTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCg                                                          >  1:3191686/1‑93 (MQ=255)                            gTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCg                                                          <  1:4206836/93‑1 (MQ=255)                                     aaaaTAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGga                                                 <  2:759411/93‑1 (MQ=255)                                        aTAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGaa                                              >  1:3891458/1‑93 (MQ=255)                                            ccGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAg                                          <  2:2608983/93‑1 (MQ=255)                                              gCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGtt                                        >  1:2305328/1‑93 (MQ=255)                                                tAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGtttt                                      <  1:915668/93‑1 (MQ=255)                                                tAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGtttt                                      <  1:3043539/93‑1 (MQ=255)                                                  aaaGCGATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTAt                                    <  2:693826/93‑1 (MQ=255)                                                       gATTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCttt                               >  1:3790169/1‑93 (MQ=255)                                                        aTTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTa                              <  2:393258/93‑1 (MQ=255)                                                        aTTTCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTa                              >  2:724087/1‑93 (MQ=255)                                                         tttCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAt                             >  1:4342953/1‑93 (MQ=255)                                                         tttCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAt                             >  1:1099802/1‑93 (MQ=255)                                                           tCTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtct                           >  1:1475401/1‑93 (MQ=255)                                                            cTTATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGATATAGTTTTATGCTTTAtctt                          >  2:3857761/1‑93 (MQ=255)                                                             ttATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttc                         >  1:3242372/1‑93 (MQ=255)                                                             ttATATTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttc                         <  2:3273413/93‑1 (MQ=255)                                                              tataTTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttct                        >  2:4205368/1‑93 (MQ=255)                                                              tataTTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttct                        >  2:2915752/1‑93 (MQ=255)                                                              tataTTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttct                        >  1:4634382/1‑93 (MQ=255)                                                              tataTTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttct                        >  1:2736922/1‑93 (MQ=255)                                                              tataTTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttct                        >  1:2646142/1‑93 (MQ=255)                                                              tataTTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTAtcttct                        >  2:177440/1‑93 (MQ=255)                                                               atatTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTtctc                       <  2:3605573/93‑1 (MQ=255)                                                               atatTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTtctc                       <  2:3968434/93‑1 (MQ=255)                                                                 atTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCta                     >  2:4400102/1‑93 (MQ=255)                                                                 atTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCta                     >  1:799436/1‑93 (MQ=255)                                                                 atTTAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCta                     >  1:663646/1‑93 (MQ=255)                                                                  tttAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGCAATAGTTTTATGCTTTATCTTctctca                    >  1:1935241/1‑91 (MQ=255)                                                                  tttAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCtaa                    >  1:3205833/1‑93 (MQ=255)                                                                   ttAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCtaat                   <  1:3023661/93‑1 (MQ=255)                                                                   ttAAATTTCTTATTCAGTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCtaat                   <  2:1133828/93‑1 (MQ=255)                                                                    tAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCtaata                  <  1:3348583/93‑1 (MQ=255)                                                                    tAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCtaata                  >  1:2785654/1‑93 (MQ=255)                                                                    tAAATTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGGTTTATGCTTTATCTTCTCttata                  >  2:2885165/1‑93 (MQ=255)                                                                     aaaTTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCtaataa                 >  1:4359775/1‑93 (MQ=255)                                                                      aaTTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAAc                <  2:2804002/93‑1 (MQ=255)                                                                      aaTTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAAc                <  2:1548464/93‑1 (MQ=255)                                                                       aTTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACt               >  1:454731/1‑93 (MQ=255)                                                                       aTTTCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACt               >  1:4220744/1‑93 (MQ=255)                                                                        tttCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACtt              >  2:2255133/1‑93 (MQ=255)                                                                         ttCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTtc             >  2:1206771/1‑93 (MQ=255)                                                                          tCTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTtcc            >  1:1900943/1‑93 (MQ=255)                                                                           cTTATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTtcct           >  1:2662854/1‑93 (MQ=255)                                                                             tATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTtcctcc         >  1:2600619/1‑93 (MQ=255)                                                                             tATTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTtcctcc         <  1:507366/93‑1 (MQ=255)                                                                              aTTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCa        <  1:1107975/93‑1 (MQ=255)                                                                              aTTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCa        <  1:4124705/93‑1 (MQ=255)                                                                              aTTCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCa        >  2:2635806/1‑93 (MQ=255)                                                                               ttCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAt       <  2:2841866/93‑1 (MQ=255)                                                                               ttCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAt       >  2:489697/1‑93 (MQ=255)                                                                               ttCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAt       >  2:834700/1‑93 (MQ=255)                                                                               ttCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAt       >  2:1293038/1‑93 (MQ=255)                                                                                tCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAtt      >  1:446668/1‑93 (MQ=255)                                                                                tCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAtt      >  2:358618/1‑93 (MQ=255)                                                                                tCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAtt      >  2:4512906/1‑93 (MQ=255)                                                                                tCATTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCAtt      >  2:3122575/1‑93 (MQ=255)                                                                                  aTTAATGACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCATTCt    >  1:693569/1‑93 (MQ=255)                                                                                    taatgaCGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCATTCTTa  >  1:3486319/1‑93 (MQ=255)                                                                                    taatgaCGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCATTCTTa  >  2:2142943/1‑93 (MQ=255)                                                                                    taatgaCGCAAGTGACAAAACAATGTTAAAAAAGTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCATTCTTa  <  1:2361187/93‑1 (MQ=255)                                                                                        |                                                                                         CCACCGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTTAAATTTCTTATTCATTAATAACGCAAGTGACAAAACAATGTTAAAAAATTCGTAACTGGGAGAAATAGTTTTATGCTTTATCTTCTCTAATAACTTCCTCCATTCTTA  >  NC_000913/3117046‑3117221

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