Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,540,351	(A)6→5	coding (189/228 nt)	feoA →	ferrous iron transporter, protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,540,346	0	A	.	100.0%	341.8 / NA	72	coding (184/228 nt)	feoA	ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (40/32);  total (40/32)

CTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAAT  >  NC_000913/3540247‑3540439                                                                                                    |                                                                                              ctTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCa                                                                                               >  1:51776/1‑100 (MQ=255)  ttGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCaa                                                                                              >  1:1865507/1‑100 (MQ=255)      cATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCaaaaag                                                                                          >  2:1551128/1‑99 (MQ=255)       aTGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCaaaaaga                                                                                         <  1:1548098/100‑3 (MQ=255)            aCCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTgg                                                                                   >  1:442304/1‑100 (MQ=255)             ccTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTgg                                                                                   >  1:1707525/1‑99 (MQ=255)                 gCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCtta                                                                              <  1:762875/100‑1 (MQ=255)                   tcctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCttatt                                                                            <  2:938771/100‑1 (MQ=255)                    cctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCttatta                                                                           >  1:750431/1‑100 (MQ=255)                    cctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCttatta                                                                           <  2:51776/100‑1 (MQ=255)                      tctTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGa                                                                         >  1:317482/1‑100 (MQ=255)                        ttttAATGTGTTGCACGTCGCGCCACTCGGCGACCCCATTCATATCGAAACGCGTCGTGTGAACCGGTTATTACGC‑AAAAAGATCTGGCCgcatta                                                                           >  2:1236014/1‑96 (MQ=255)                         tttAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGTGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGt                                                                      >  1:387551/1‑100 (MQ=255)                           tAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTgg                                                                    >  2:43119/1‑100 (MQ=255)                               gtggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTgccagc                                                                >  2:1564532/1‑95 (MQ=255)                               gtggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGc                                                                >  2:2012656/1‑100 (MQ=255)                               gtggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGc                                                                >  1:1874790/1‑100 (MQ=255)                               gtggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTAGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGc                                                                <  2:2020326/100‑1 (MQ=255)                                tggtgTGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCg                                                               >  1:1606473/1‑100 (MQ=255)                                tggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCg                                                               >  1:1316568/1‑100 (MQ=255)                                tggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCg                                                               >  1:1174616/1‑100 (MQ=255)                                tggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCg                                                               >  1:1317107/1‑100 (MQ=255)                                 ggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCgg                                                              >  1:211669/1‑100 (MQ=255)                                     cgcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTc                                                          <  2:738185/100‑1 (MQ=255)                                     cgcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTc                                                          <  2:750428/100‑1 (MQ=255)                                      gcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTcc                                                         <  2:1874773/100‑1 (MQ=255)                                        gTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTg                                                       <  2:317481/100‑1 (MQ=255)                                         tCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGt                                                      >  2:2016138/1‑100 (MQ=255)                                         tCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGt                                                      >  2:534860/1‑100 (MQ=255)                                          cGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGtt                                                     >  2:608581/1‑100 (MQ=255)                                          cGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGtt                                                     >  1:1687871/1‑100 (MQ=255)                                          cGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGtt                                                     <  1:1789671/100‑1 (MQ=255)                                           gCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGtt                                                     >  2:386307/1‑99 (MQ=255)                                           gCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTa                                                    <  1:347419/100‑1 (MQ=255)                                           gCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTa                                                    >  1:444795/1‑100 (MQ=255)                                             tCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAAt                                                  <  1:1948592/100‑1 (MQ=255)                                              ccACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATa                                                 >  2:820406/1‑100 (MQ=255)                                              ccACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTATCCTGTTAATa                                                 >  2:396557/1‑100 (MQ=255)                                                  tCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGt                                             <  2:546566/100‑1 (MQ=255)                                                   cGGCGACCCCATTCATATCGAAACCCGTTGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTg                                            <  2:105954/100‑1 (MQ=255)                                                   cGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTg                                            <  2:1604981/100‑1 (MQ=255)                                                    ggCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGa                                           <  1:408112/100‑1 (MQ=255)                                                          cccATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATaacaa                                     >  1:1378900/1‑100 (MQ=255)                                                           ccATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATaacaac                                    >  1:870945/1‑100 (MQ=255)                                                            cATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAacaaca                                   <  2:1606064/100‑1 (MQ=255)                                                            cATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAacaaca                                   >  1:217798/1‑100 (MQ=255)                                                            cATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAacaaca                                   >  2:1310053/1‑100 (MQ=255)                                                              ttCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAAt                                 >  1:856804/1‑100 (MQ=255)                                                              ttCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAAt                                 >  1:1392145/1‑100 (MQ=255)                                                                cATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGa                               >  2:1195684/1‑100 (MQ=255)                                                                cATATCGAAACCCGTCGCGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGa                               >  1:334526/1‑100 (MQ=255)                                                                  tatCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTTCTGTTAATACGGTGATAACAACAATGa                               <  1:1479926/98‑1 (MQ=255)                                                                  tatCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaa                             <  1:1564542/100‑1 (MQ=255)                                                                   atCGAAACTCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAATCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaa                            <  2:1243417/100‑1 (MQ=255)                                                                   atCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaa                            <  2:1217187/100‑1 (MQ=255)                                                                     cGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaaaa                          >  1:523909/1‑100 (MQ=255)                                                                      gAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAt                         >  2:1623259/1‑100 (MQ=255)                                                                       aaaCCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAtt                        >  1:1463254/1‑100 (MQ=255)                                                                          cccGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAAc                     <  1:1114868/100‑1 (MQ=255)                                                                           ccGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAAcc                    >  1:1481373/1‑100 (MQ=255)                                                                           ccGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAAcc                    <  1:868851/100‑1 (MQ=255)                                                                             gtcgtGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCAt                  >  1:1565793/1‑100 (MQ=255)                                                                                gtgtGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTgg               <  1:27162/100‑1 (MQ=255)                                                                                 tgtgAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGc              <  2:211668/100‑1 (MQ=255)                                                                                 tgtgAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGc              <  1:1509612/100‑1 (MQ=255)                                                                                  gtgAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCt             <  2:1314894/100‑1 (MQ=255)                                                                                   tgAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTaagccaatggtt            <  1:598470/100‑11 (MQ=255)                                                                                   tgAGCCTGATATTACGC‑AAAAAGATCTGGCCTTATTAGAAGAGGAAACGGTTTCCTGTTAATACGGTGATAAAAAAAAAAAAAAAATTAACCATTGGCt             <  2:444794/99‑1 (MQ=255)                                                                                     aGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTaa          <  2:334525/100‑1 (MQ=255)                                                                                     aGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTaa          <  1:43119/100‑1 (MQ=255)                                                                                      gCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAAt         >  1:305205/1‑100 (MQ=255)                                                                                      gCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAAt         >  1:1071645/1‑100 (MQ=255)                                                                                       ccTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAAtt        <  2:1317099/100‑1 (MQ=255)                                                                                       ccTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAAtt        <  2:1316560/100‑1 (MQ=255)                                                                                       ccTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAAtt        >  2:875895/1‑100 (MQ=255)                                                                                             aTTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAAt  >  1:1887998/1‑100 (MQ=255)                                                                                                    |                                                                                              CTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAAT  >  NC_000913/3540247‑3540439

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