Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,413,811	Δ1 bp	coding (93/1158 nt)	acrE →	cytoplasmic membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,413,809	0	A	.	100.0%	151.2 / NA	38	coding (91/1158 nt)	acrE	cytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (15/23);  total (15/23)

CTCCTGCCCTCCTTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAGAAAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAGTTAA  >  W3110S.gb/3413743‑3413873                                                                   |                                                                 ctcctGCCCTCCTTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGgagaagagaagct                                                              <  1:4553398/71‑4 (MQ=255) ctcctGCCCTCCTTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGgagaagagaagct                                                              <  1:21276/71‑4 (MQ=255)       cccTCCTTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTc                                                       <  1:4112596/71‑1 (MQ=255)           cctTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGt                                                    >  1:152788/1‑70 (MQ=255)            ctTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACTTCGGTg                                                   <  1:4728/70‑1 (MQ=255)            ctTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGa                                                  <  1:1294814/71‑1 (MQ=255)              ttATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAAc                                                >  1:2042878/1‑71 (MQ=255)              ttATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAAc                                                >  1:1294180/1‑71 (MQ=255)               tATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAAcc                                               >  1:1524509/1‑71 (MQ=255)               tATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAAcc                                               >  1:2070572/1‑71 (MQ=255)                   cTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCa                                            <  1:215493/70‑1 (MQ=255)                   cTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCa                                            <  1:339038/70‑1 (MQ=255)                   cTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCa                                            <  1:4403433/70‑1 (MQ=255)                    tGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCa                                            <  1:3521088/69‑1 (MQ=255)                     gATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCa                                                           >  1:434406/1‑53 (MQ=255)                      aTCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGt                                         >  1:3249632/1‑70 (MQ=255)                        ctcCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTAc                                      >  1:4213961/1‑71 (MQ=255)                           cgcgGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGt                                   <  1:3725339/71‑1 (MQ=255)                              ggCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCa                                <  1:2835531/71‑1 (MQ=255)                              ggCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCa                                <  1:786911/71‑1 (MQ=255)                               gCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCat                               <  1:3963856/71‑1 (MQ=255)                                cTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCa                                            <  1:3701966/57‑1 (MQ=255)                                        gCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:1381085/70‑1 (MQ=255)                                        gCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:3933932/70‑1 (MQ=255)                                        gCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:4197515/70‑1 (MQ=255)                                        gCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:627166/70‑1 (MQ=255)                                        gCCGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:1433245/70‑1 (MQ=255)                                         ccGGTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:923000/69‑1 (MQ=255)                                            gTTGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTg                                                   <  1:2813967/38‑1 (MQ=38)                                              tGTAACGATAAGGGAGAAGAG‑AAGCTCACTTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGcc                >  1:3274389/1‑71 (MQ=255)                                              tGTAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGcc                >  1:2687663/1‑71 (MQ=255)                                                tAACGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTaaa                       <  1:4406630/62‑1 (MQ=255)                                                  aCGATAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAATCGGCCCCGt            >  1:4177738/1‑71 (MQ=255)                                                     aTAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCata                              <  1:1274094/50‑1 (MQ=38)                                                     aTAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGT‑AAAACGGCCCCGTTAGa        <  1:3804975/71‑1 (MQ=255)                                                      tAAGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGa        <  1:294371/71‑1 (MQ=255)                                                        aGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAg      >  1:4591429/1‑71 (MQ=255)                                                        aGGGAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAg      >  1:842492/1‑71 (MQ=255)                                                          ggAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGt                                   >  1:1116675/1‑40 (MQ=37)                                                          ggAGAAGAG‑AAGCTCACGTCGGTGAACCGCAGGTTACCGTTc                                 >  1:1629629/1‑42 (MQ=37)                                                           gagaagag‑aaGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAGtt    >  1:4060955/1‑70 (MQ=255)                                                            agaagag‑aagCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAGTTaa  <  1:4113275/71‑1 (MQ=255)                                                            agaagag‑aagCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAGTTaa  <  1:4223470/71‑1 (MQ=255)                                                            agaagag‑aagCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAGTTaa  <  1:2320749/71‑1 (MQ=255)                                                                   |                                                                 CTCCTGCCCTCCTTTATTCTGATCTCCGCGGCTTTAATCGCCGGTTGTAACGATAAGGGAGAAGAGAAAGCTCACGTCGGTGAACCGCAGGTTACCGTTCATATTGTAAAAACGGCCCCGTTAGAAGTTAA  >  W3110S.gb/3413743‑3413873

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