Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,240,670	Δ1 bp	coding (55/441 nt)	BW25113_RS06225 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,240,667	0	G	.	100.0%	163.5 / NA	37	coding (52/441 nt)	BW25113_RS06225	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (19/18);  total (19/18)

CAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACC  >  NZ_CP009273/1240578‑1240750                                                                                          |                                                                                    caGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATg                                                                                     >  2:1350726/1‑90 (MQ=255)     aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATgggcg                                                                                 >  1:1050872/1‑88 (MQ=255)          tttttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAg                                                                           >  2:1161824/1‑90 (MQ=255)             ttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGa                                                                        <  2:2555174/90‑1 (MQ=255)               aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAgg                                                                      <  2:226938/90‑1 (MQ=255)                 ttattGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCg                                                                    >  2:2781811/1‑90 (MQ=255)                        ttCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTgg                                                             >  2:2126185/1‑90 (MQ=255)                          cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCt                                                           <  2:618258/90‑1 (MQ=255)                                   tGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTAtt                                                  <  1:1350726/90‑1 (MQ=255)                                     gtgtGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTGa                                                >  2:546186/1‑90 (MQ=255)                                     gtgtGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGa                                                <  1:2377759/90‑1 (MQ=255)                                      tgtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                               <  1:1649871/90‑1 (MQ=255)                                      tgtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                               <  2:561109/90‑1 (MQ=255)                                      tgtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                               <  2:2956922/90‑1 (MQ=255)                                        tgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTaaaa                                             >  1:1620287/1‑88 (MQ=255)                                           ttCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACAttt                                          <  1:2618702/90‑1 (MQ=255)                                                aaaaTAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAAt                                     <  1:1480189/90‑1 (MQ=255)                                                     aaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATc                                <  2:418445/90‑1 (MQ=255)                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCt                                                           <  2:2287871/55‑1 (MQ=255)                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCt                                                           >  1:2287871/1‑55 (MQ=255)                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                               <  1:2496991/67‑1 (MQ=255)                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGaa                                               >  2:2496991/1‑67 (MQ=255)                                                                 cAGTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGc                    <  2:677503/90‑1 (MQ=255)                                                                   gTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAACATTTTGAATTGATCAATCCTTTTTGCAg                  >  2:1680653/1‑90 (MQ=255)                                                                   gTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTATTTAAAATTTTGAATTGATCAAACCTTCTGGCAg                  >  2:2072555/1‑90 (MQ=255)                                                                   gTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTTTTTAAAATTTAGAAATTATCAATCCTTCTGGcgg                  >  1:1356572/1‑88 (MQ=255)                                                                   gTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTTGCAg                  >  2:1198861/1‑90 (MQ=255)                                                                   gTAGCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAg                  >  2:655441/1‑90 (MQ=255)                                                                      gCCAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGAtt               <  2:1050872/90‑1 (MQ=255)                                                                        cAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAAAATTTAGAATTGATCAATCCTTCTGGCAGATTaa             >  1:521481/1‑90 (MQ=255)                                                                        cAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAACATTTAGAATTGATCAATCCTTCTGGCAGATTaa             >  2:2007162/1‑90 (MQ=255)                                                                        cAGTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTaa             <  1:1561240/90‑1 (MQ=255)                                                                          gTCAGGCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGc           <  2:1700314/90‑1 (MQ=255)                                                                               gCGTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGCTGGGCCTATTTTACATTTAAAATTTTTCAAATCTTCTGGGGGAATAAGTCcaaa      >  1:2122480/1‑90 (MQ=255)                                                                               gCGTTAAAAAT‑GGGCGGATAGTGAGGCGGATTTTGCTGCGCTTATTTAAAATTTAGAATTGATCAATCCTTCTTGGAGATTAAGCCaaaa      >  1:1967870/1‑90 (MQ=255)                                                                                 gTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAAAATTTAGAATTGATCAAACCTTTTGGCAGATTAAGCCaaaaaa    >  1:639868/1‑88 (MQ=255)                                                                                 gTTAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATa    >  2:97225/1‑90 (MQ=255)                                                                                   tAAAAAT‑GGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTAACATTTAGAATTGATCAATCCTTCTGGCAGATTTAGCCAAAAAAcc  >  1:10899/1‑90 (MQ=255)                                                                                   tAAAAAT‑GGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATAcc  <  2:1769827/90‑1 (MQ=255)                                                                                          |                                                                                    CAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACC  >  NZ_CP009273/1240578‑1240750

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