Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,649,516	A→G	L194L (TTA→TTG)	arsB →	arsenite/antimonite:H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,649,516	0	A	G	100.0%	179.5 / NA	50	L194L (TTA→TTG)	arsB	arsenite/antimonite:H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/30);  total (20/30)

GCCTGCCGCTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTACATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGATTTTT  >  NC_000913/3649386‑3649650                                                                                                                                   |                                                                                                                                       gCCTGCCGCTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTAt                                                                                                                               <  2:814670/140‑1 (MQ=255)  ccTGCCGCTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTAtt                                                                                                                              >  1:767754/1‑140 (MQ=255)    tGCCGCTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTAtttt                                                                                                                            <  1:46113/140‑1 (MQ=255)      ccGCTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCg                                                                                                                          >  1:327087/1‑140 (MQ=255)         cTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCaa                                                                                                                       <  2:25801/140‑1 (MQ=255)            aTTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGa                                                                                                                    >  1:417095/1‑140 (MQ=255)              tGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGata                                                                                                                  <  1:454051/140‑1 (MQ=255)                   ccAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCg                                                                                                             <  1:817956/140‑1 (MQ=255)                       ccTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGa                                                                                                         <  2:66722/140‑1 (MQ=255)                          ggTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGaa                                                                                                        >  2:448367/1‑138 (MQ=255)                          ggTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGaa                                                                                                        <  1:448367/138‑1 (MQ=255)                            tGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTAc                                                                                                    >  2:666706/1‑140 (MQ=255)                                 atCGTTTCCGCTGATTTCTTTGGCCTCGGCTCTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTa                                                                                                     >  2:715992/1‑134 (MQ=255)                                 atCGTTTCCGCTGATTTCTTTGGCCTCGGCTCTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTa                                                                                                     <  1:715992/134‑1 (MQ=255)                                     tttCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGgcg                                                                                           >  2:647966/1‑140 (MQ=255)                                     tttCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGgcg                                                                                           <  2:17331/140‑1 (MQ=255)                                     tttCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGgcg                                                                                           >  1:743987/1‑140 (MQ=255)                                      ttCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGgcgc                                                                                          >  1:359838/1‑140 (MQ=255)                                             gATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGaa                                                                                   >  2:434903/1‑140 (MQ=255)                                                          tCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGaa                                                                      <  1:446137/140‑1 (MQ=255)                                                           cGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAg                                                                     <  1:143107/140‑1 (MQ=255)                                                            ggCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGc                                                                    <  2:417095/140‑1 (MQ=255)                                                                ttCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATc                                                                <  1:647966/140‑1 (MQ=255)                                                                 tCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCa                                                               <  2:351118/140‑1 (MQ=255)                                                                  cgcgAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCaa                                                              <  2:767754/140‑1 (MQ=255)                                                                       aTACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATc                                                         <  1:898788/140‑1 (MQ=255)                                                                            ccTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCa                                                               >  2:571215/1‑129 (MQ=255)                                                                            ccTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCa                                                               <  1:571215/129‑1 (MQ=255)                                                                            ccTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCt                                                    <  2:165478/140‑1 (MQ=255)                                                                                 gTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGtt                                               >  2:918311/1‑140 (MQ=255)                                                                                    aTGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCaa                                            >  1:110905/1‑140 (MQ=255)                                                                                           cGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAActggc                                     >  1:189559/1‑140 (MQ=255)                                                                                             gTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAActggctg                                   <  1:492876/140‑1 (MQ=255)                                                                                               ggATATCGCCGCGATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTggg                                 <  2:189559/140‑1 (MQ=255)                                                                                                        cgcgATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTAc                        >  2:349221/1‑140 (MQ=255)                                                                                                        cgcgATTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTAc                        >  1:858398/1‑140 (MQ=255)                                                                                                            aTTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGtttt                          <  2:57350/134‑1 (MQ=255)                                                                                                            aTTGTTGCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGtttt                          >  1:57350/1‑134 (MQ=255)                                                                                                                 tgCCACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTctgc               <  2:835050/140‑1 (MQ=255)                                                                                                                   ccACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTctgctg             <  1:278823/140‑1 (MQ=255)                                                                                                                    cACGCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCtgg            <  2:110905/140‑1 (MQ=255)                                                                                                                       gCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCt                    <  2:641734/129‑1 (MQ=255)                                                                                                                       gCTGGTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTAGGTTGTTTTACTGCt                    >  1:641734/1‑129 (MQ=255)                                                                                                                          ggTGATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGAt      >  1:215966/1‑140 (MQ=255)                                                                                                                             gATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTc                  >  1:902243/1‑125 (MQ=255)                                                                                                                             gATGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTc                  <  2:902243/125‑1 (MQ=255)                                                                                                                              aTGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGAttttt  <  1:349221/140‑1 (MQ=255)                                                                                                                              aTGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGAttttt  <  1:918311/140‑1 (MQ=255)                                                                                                                              aTGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGAttttt  <  2:858398/140‑1 (MQ=255)                                                                                                                              aTGTTGCATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCCGAAATCTCCCGCAGAAGCGATCAACGACCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGAttttt  <  2:633424/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       GCCTGCCGCTTATTGTCTCCAACCTGGTGAATATCGTTTCCGCTGATTTCTTTGGCCTCGGCTTTCGCGAATACGCCTCGGTGATGGTGCCGGTGGATATCGCCGCGATTGTTGCCACGCTGGTGATGTTACATCTCTATTTTCGCAAAGATATTCCGCAGAACTACGATATGGCGCTGCTGAAATCTCCCGCAGAAGCGATCAAAGATCCTGCTACGTTCAAAACTGGCTGGGTTGTTTTACTGCTTCTGCTGGTGGGATTTTT  >  NC_000913/3649386‑3649650

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