Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,477,321	Δ1 bp	coding (12/1182 nt)	yjgN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,477,321	0	A	.	100.0%	207.6 / NA	46	coding (12/1182 nt)	yjgN	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (27/19);  total (27/19)

TATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGT  >  NC_000913/4477250‑4477389                                                                        |                                                                     tatataAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCg                                                                <  2:602583/76‑1 (MQ=255)      aaGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCa                                                           >  2:2431225/1‑76 (MQ=255)       aGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCAt                                                          <  2:3989531/76‑1 (MQ=255)       aGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCAt                                                          <  1:709517/76‑1 (MQ=255)         ttATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTc                                                        >  2:1134517/1‑76 (MQ=255)         ttATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTc                                                        >  1:4357762/1‑76 (MQ=255)          tataTCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                                                       >  1:1867464/1‑76 (MQ=255)          tataTCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCg                                                       >  1:1888757/1‑76 (MQ=255)                  tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt                                               <  1:2719885/76‑1 (MQ=255)                      ttttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCAt                                           >  2:984744/1‑76 (MQ=255)                       tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg                                          <  1:2171652/76‑1 (MQ=255)                       tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg                                          >  1:3662328/1‑76 (MQ=255)                        ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATgg                                         >  1:1240482/1‑76 (MQ=255)                        ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATgg                                         >  1:1139031/1‑76 (MQ=255)                         tttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt                                        <  2:2480564/76‑1 (MQ=255)                         tttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt                                        >  2:2154756/1‑76 (MQ=255)                             aattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTAcg                                     <  2:2310276/74‑2 (MQ=255)                             aattCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACa                                     <  1:2203572/74‑1 (MQ=255)                              attCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg                                    >  1:2197756/2‑76 (MQ=255)                              attCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAg                                    >  2:1561482/2‑76 (MQ=255)                               ttCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTAGTTTGTTTTTCATGGTACAgg                                   >  1:2062593/1‑76 (MQ=255)                                  aGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTga                                <  2:2040690/76‑1 (MQ=255)                                    ggAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgaga                              >  2:3871036/1‑76 (MQ=255)                                    ggAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTgaga                              <  1:2016122/76‑1 (MQ=255)                                        tttttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGatt                          >  2:454552/1‑75 (MQ=255)                                          tttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtt                        >  2:883881/1‑76 (MQ=255)                                          tttATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtt                        >  2:4771410/1‑76 (MQ=255)                                          tttAGGGCTCAAGATATTAATGAAATAGCT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGCACAGGTGAGAGCTAtt                        >  2:4772836/1‑76 (MQ=255)                                            tATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtttt                      >  2:1959976/1‑76 (MQ=255)                                            tATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAtttt                      >  1:615621/1‑76 (MQ=255)                                             aTGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATAttttt                     >  2:3191314/1‑75 (MQ=255)                                              tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCt                    <  2:2781826/76‑1 (MQ=255)                                              tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCt                    <  2:2804865/76‑1 (MQ=255)                                              tGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCt                    >  2:3388594/1‑76 (MQ=255)                                                gCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCttt                  >  1:4590198/1‑75 (MQ=255)                                                  tCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAtt                >  2:676221/1‑76 (MQ=255)                                                  tCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTAtt                <  2:1623012/76‑1 (MQ=255)                                                      gTTATTAATGCCATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGCGAGATATTTTCTTCtttttg            <  1:280417/76‑4 (MQ=255)                                                        tattaATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTtgtgtg          <  2:3817009/76‑1 (MQ=255)                                                            aaTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGa      <  1:3871036/76‑1 (MQ=255)                                                             aTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGaa     <  1:3424331/76‑1 (MQ=255)                                                             aTGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGaa     >  1:609210/1‑76 (MQ=255)                                                              tGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAt    >  1:3184810/1‑76 (MQ=255)                                                               gAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtg   <  1:3531453/76‑1 (MQ=255)                                                               gAAAGGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtg   <  1:3540678/76‑1 (MQ=255)                                                                aaaTGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAAtgt  >  2:3752325/1‑76 (MQ=255)                                                                        |                                                                     TATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGTACAGGTGAGAGATATTTTCTTATTTGTGTGGTGAATGT  >  NC_000913/4477250‑4477389

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