Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,540,736	A→G	E110E (GAA→GAG)	feoB →	ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,540,736	0	A	G	100.0%	137.5 / NA	37	E110E (GAA→GAG)	feoB	ferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/21);  total (16/21)

CGCAGACCTCGCTCGATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAACTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTCAACCCGTGGTCGCGGTATTGAAGC  >  NC_000913/3540609‑3540879                                                                                                                                |                                                                                                                                                cGCAGACCTCGCTCGATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTgg                                                                                                                           <  1:314619/150‑1 (MQ=255)      aCCTCGCTCGATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGa                                                                                                                     <  2:268789/151‑1 (MQ=255)          cgctcgATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACa                                                                                                                  <  2:77895/150‑1 (MQ=255)            ctcgATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGc                                                                                                               >  1:112877/1‑151 (MQ=255)              cgATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTc                                                                                                             <  1:216719/151‑1 (MQ=255)                       aaTCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGcc                                                                                                    <  1:603891/151‑1 (MQ=255)                          cGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCgag                                                                                                 <  1:229249/151‑1 (MQ=255)                           gCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCg                                                                                                   >  2:211131/1‑148 (MQ=255)                                tCACTACATTTTGAGTGGCGACGCCGACCTTCTGATTAACGTGGTGGATGCGTCTACCCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGc                                                                                             <  2:223808/149‑1 (MQ=255)                                  aCTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCaaa                                                                                          >  1:460594/1‑150 (MQ=255)                                     aCATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAat                                                                                        <  1:226075/149‑1 (MQ=255)                                      cATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATAtt                                                                                     >  2:378792/1‑151 (MQ=255)                                        ttttGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCg                                                                                   >  2:512834/1‑151 (MQ=255)                                         tttGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGt                                                                                  >  2:380325/1‑151 (MQ=255)                                            gAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTAt                                                                                <  1:23286/150‑1 (MQ=255)                                            gAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTAt                                                                                <  1:221321/150‑1 (MQ=255)                                              gTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTg                                                                              >  1:109323/1‑150 (MQ=255)                                                ggCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGaa                                                                            <  1:555361/150‑1 (MQ=255)                                                ggCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGaa                                                                            <  2:611478/150‑1 (MQ=255)                                                ggCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGaa                                                                            <  2:184108/150‑1 (MQ=255)                                                ggCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGaa                                                                            <  1:277065/150‑1 (MQ=255)                                                       ccGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGc                                                                    <  1:519450/151‑1 (MQ=255)                                                             tgctgATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTc                                                              >  2:78963/1‑151 (MQ=255)                                                              gctgATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTc                                                              <  1:380332/150‑1 (MQ=255)                                                                        gtggtgGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTgg                                                   >  2:231212/1‑151 (MQ=255)                                                                          ggtggATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGc                                                 >  2:339119/1‑151 (MQ=255)                                                                           gtggATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCt                                                <  1:23568/151‑1 (MQ=255)                                                                             ggATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTg                                               <  1:3708/150‑1 (MQ=255)                                                                                    tCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTg                                        >  2:93616/1‑150 (MQ=255)                                                                                       aaCCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTTGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATcc                                    >  1:598073/1‑151 (MQ=255)                                                                                       aaCCTGGAGCGTAGCCTGTACCTGACGGTACAGCGGCTGGAGCTCGGCATGCCCTGCATTGTGGCACTGAACATGCTCGACAGTGCCGAGAAGCAAAATATGCGTATTGAAAGTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGa                                       <  2:329555/148‑1 (MQ=255)                                                                                        aCCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATccc                                   >  2:504758/1‑151 (MQ=255)                                                                                                 gTAACCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTc                          >  1:583799/1‑151 (MQ=255)                                                                                                    aCCTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTCAAc                       >  2:70442/1‑151 (MQ=255)                                                                                                     ccTGTACCTGACGCTACAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTCAAcc                      >  2:419286/1‑151 (MQ=255)                                                                                                                     cAACTGCTGGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTCAACCCGTGGTCGCGGTAtt       <  1:419293/150‑1 (MQ=255)                                                                                                                          gctgGAGCTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTCAACCCGTGGTCGCGGTATTGaatc  <  1:483826/150‑3 (MQ=255)                                                                                                                                |                                                                                                                                                CGCAGACCTCGCTCGATGAGCAAATCGCCTGTCACTACATTTTGAGTGGCGACGCCGACCTGCTGATTAACGTGGTGGATGCGTCTAACCTTGAGCGTAACCTGTACCTGACGCTACAACTGCTGGAACTCGGCATTCCCTGCATTGTGGCACTGAACATGCTCGACATTGCCGAGAAGCAAAATATTCGTATTGAAATTGATGCTCTGTCGGCGCGTCTGGGCTGTCCGGTGATCCCGCTGGTTTCAACCCGTGGTCGCGGTATTGAAGC  >  NC_000913/3540609‑3540879

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