Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AC_000091	4,017,654	C→T	R1190R (CGG→CGA)	rhsB ←	rhsB element core protein RshB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AC_000091	4,017,654	0	C	T	100.0%	113.3 / NA	40	R1190R (CGG→CGA)	rhsB	rhsB element core protein RshB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

GACTCCTCATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCC  >  AC_000091/4017625‑4017738                              |                                                                                     gACTCCTCATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca                          <  2:401215/90‑1 (MQ=14)  aCTCCTCATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTcatcat                         <  2:2251372/90‑1 (MQ=17)  aCTCCTCATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTcatcat                         >  2:1032071/1‑90 (MQ=17)   ctcctcATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTCGCCCATTCATCata                        >  1:3044562/1‑90 (MQ=11)    tcctcATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCatat                       >  2:1782690/1‑90 (MQ=18)    tcctcATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCatat                       >  2:1474265/1‑90 (MQ=18)    tcctcATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGGTGCCCCATTCATCatat                       >  2:388858/1‑90 (MQ=255)      ctcATCATACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTc                     >  1:783382/1‑90 (MQ=18)       tcatcaTACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCt                    <  1:1142248/90‑1 (MQ=18)        catcatACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTg                   <  2:3269731/90‑1 (MQ=18)        catcatACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTg                   >  1:610748/1‑90 (MQ=18)         atcatACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTgc                  <  1:1645564/90‑1 (MQ=18)         atcatACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTgc                  >  1:1656972/1‑90 (MQ=18)          tcatACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTgcg                 <  1:2145403/90‑1 (MQ=18)           catACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTgcgc                >  1:915203/1‑90 (MQ=18)            atACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGca               >  1:2075129/1‑90 (MQ=18)            atACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGca               >  2:2273239/1‑90 (MQ=18)            atACTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGca               >  1:87755/1‑90 (MQ=18)              aCTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacc             >  1:3049621/1‑90 (MQ=18)              aCTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacc             <  1:3089287/90‑1 (MQ=18)              aCTGCTGCCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacc             <  2:3092029/90‑1 (MQ=25)                tgctgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac           >  1:1143103/1‑90 (MQ=25)                tgctgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac           <  1:2159813/90‑1 (MQ=25)                tgctgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac           >  1:2003045/1‑90 (MQ=25)                 gctgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACg          <  2:3268417/90‑1 (MQ=255)                 gctgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACg          <  2:2783411/90‑1 (MQ=255)                  ctgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGc         >  1:1531859/1‑90 (MQ=255)                   tgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt        <  1:339013/90‑1 (MQ=255)                   tgcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt        <  2:2094868/90‑1 (MQ=255)                    gcCCCGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTg       >  2:1862252/1‑90 (MQ=255)                     ccccGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGt      <  1:90711/90‑1 (MQ=255)                     ccccGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGt      >  1:1513073/1‑90 (MQ=255)                      cccGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGtt     <  2:167380/90‑1 (MQ=255)                      cccGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGtt     >  1:2346827/1‑90 (MQ=255)                      cccGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGtt     <  1:1300404/90‑1 (MQ=255)                      cccGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGtt     <  2:749555/90‑1 (MQ=255)                        cGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGc   >  1:1751417/1‑90 (MQ=255)                        cGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGc   >  1:1285816/1‑90 (MQ=255)                        cGGCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGc   <  2:527556/90‑1 (MQ=255)                         ggCAGTCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGcc  <  2:371354/90‑1 (MQ=255)                              |                                                                                     GACTCCTCATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCC  >  AC_000091/4017625‑4017738

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