Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,477,281	Δ1 bp	intergenic (‑391/‑29)	b4256 ← / → yjgN	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,477,279	0	T	.	100.0%	227.0 / NA	50	intergenic (‑389/‑31)	b4256/yjgN	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (23/27);  total (23/27)

TACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGT  >  NC_000913/4477204‑4477351                                                                            |                                                                         tACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAt                                                                          >  2:4730833/1‑76 (MQ=255)  aCTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAtt                                                                         >  1:2748120/1‑76 (MQ=255)   cTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAttc                                                                        <  1:4643880/76‑2 (MQ=255)     gATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAAttcag                                                                      <  2:1117532/76‑4 (MQ=255)          gagTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAt                                                                <  2:2995055/76‑1 (MQ=255)           agTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAtt                                                               <  1:479593/76‑1 (MQ=255)           agTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAtt                                                               <  1:4013871/76‑1 (MQ=255)             tATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAAtttt                                                             <  1:319836/76‑1 (MQ=255)               ttGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTa                                                           >  1:863833/1‑76 (MQ=255)                tGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTAt                                                          <  2:1321804/76‑1 (MQ=255)                tGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTAt                                                          <  2:1315426/76‑1 (MQ=255)                tGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTAt                                                          <  1:4833214/76‑1 (MQ=255)                   aaCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGc                                                       >  2:4726970/1‑76 (MQ=255)                   aaCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGc                                                       >  2:3621332/1‑76 (MQ=255)                    aCGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCt                                                      >  1:2032261/1‑76 (MQ=255)                     cGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTc                                                     <  2:1927652/76‑1 (MQ=255)                       tGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCaa                                                   <  1:3228051/76‑1 (MQ=255)                       tGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCaa                                                   <  2:2936992/76‑1 (MQ=255)                            gAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGttat                                              >  2:64151/1‑76 (MQ=255)                            gAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGttat                                              >  1:2257573/1‑76 (MQ=255)                             aaCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGttatt                                             <  2:5456304/76‑1 (MQ=255)                              aCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGttatta                                            <  2:2067122/76‑1 (MQ=255)                              aCATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGttatta                                            <  2:3160936/76‑1 (MQ=255)                               cATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTaa                                           <  1:1349474/76‑1 (MQ=255)                                  tttACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGa                                        >  1:5368791/1‑76 (MQ=255)                                   ttACTGTCTGGATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGaa                                       >  2:3211826/1‑76 (MQ=255)                                     aCTGTCTGTATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAAt                                     <  1:22959/76‑1 (MQ=255)                                      ctgtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATg                                    <  1:4861387/76‑1 (MQ=255)                                        gtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGa                                  <  1:5658204/76‑1 (MQ=255)                                        gtctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGa                                  >  2:2448829/1‑76 (MQ=255)                                         tctgtATATATCAGTTATATCAATGGCTTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAt                                 >  2:676782/1‑76 (MQ=255)                                         tctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAt                                 <  2:863833/76‑1 (MQ=255)                                         tctgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAt                                 <  1:557133/76‑1 (MQ=255)                                           tgtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtgt                               >  2:3435580/1‑74 (MQ=255)                                            gtATATATAAGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAtgtt                              >  1:3548194/1‑73 (MQ=255)                                                   taAGTTATATCAATGTATTTTTTAA‑TTCCGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTccc                      <  2:3317196/76‑1 (MQ=255)                                                    aaGTTATATCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCa                     >  1:2984113/1‑76 (MQ=255)                                                         atatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgt                >  2:4407521/1‑76 (MQ=255)                                                         atatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgt                >  2:4365006/1‑76 (MQ=255)                                                          tatCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtt               <  1:2341833/76‑1 (MQ=255)                                                            tCAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttg             >  2:2092250/1‑76 (MQ=255)                                                             cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgt            >  2:1782507/1‑76 (MQ=255)                                                             cAATGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATGAATGAAATGGAT‑GTTCCGTCCCATTCgtttgt            <  2:2506631/76‑1 (MQ=255)                                                              aaTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt           >  1:1750205/1‑76 (MQ=255)                                                              aaTGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCgtttgtt           >  1:549859/1‑76 (MQ=255)                                                                tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt         >  1:3121426/1‑76 (MQ=255)                                                                tGGATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGtttt         <  1:5810702/76‑1 (MQ=255)                                                                 ggATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGttttt        <  1:64151/76‑1 (MQ=255)                                                                  gATTTTTTAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTc       >  1:4293487/1‑76 (MQ=255)                                                                    ttttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTTGTTTGTTTTTCAt     <  2:2257573/76‑1 (MQ=255)                                                                     tttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATg    <  1:3643951/76‑1 (MQ=255)                                                                      ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCCTgg   >  2:4989531/1‑76 (MQ=255)                                                                      ttttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTCTGTTTTTCATgg   >  1:4791917/1‑76 (MQ=255)                                                                       tttAA‑TTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGAT‑GTTCCGTCCCATTCGTTTGTTTTTCATGGt  <  1:5695099/76‑1 (MQ=255)                                                                            |                                                                         TACTGATAAGAGTATTGTAACGTGGCAGAACATTTTACTGTCTGTATATATAAGTTATATCAATGGATTTTTTAATTTCAGGGAATTTTTATGGCTCAAGTTATTAATGAAATGGATAGTTCCGTCCCATTCGTTTGTTTTTCATGGT  >  NC_000913/4477204‑4477351

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