Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	169,686	Δ1 bp	coding (2203/2244 nt)	fhuA →	ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	169,686	0	G	.	100.0%	194.4 / NA	43	coding (2203/2244 nt)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (22/21);  total (22/21)

GTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAA  >  NC_000913/169543‑169824                                                                                                                                                |                                                                                                                                           gTTATACGGGCGAGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGGGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGTCTGCTTCTGGGGC‑CAGAAc                                                                                                                                      <  3:39614/149‑1 (MQ=255)    ataCGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGt                                                                                                                                    >  4:86660/1‑148 (MQ=255)     taCGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTc                                                                                                                                   <  2:118328/148‑1 (MQ=255)      aCGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAg                                                                                                                                 >  1:69653/1‑149 (MQ=255)        ggTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGt                                                                                                                               <  3:86660/149‑1 (MQ=255)         gTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTc                                                                                                                              >  8:4298/1‑149 (MQ=255)                                       gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGt                                                                                                                                    <  1:1841/113‑1 (MQ=255)                                       gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGt                                                                                                                                    >  2:1841/1‑113 (MQ=255)                                                 gCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTggg                                                                                      <  2:69653/149‑1 (MQ=255)                                                 gCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTggg                                                                                      >  2:86775/1‑149 (MQ=255)                                                      gctggctCCAACGTGGCGCGGCATGTGAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCAc                                                                                  <  7:4298/148‑1 (MQ=255)                                                            tCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTc                                                                           >  5:113083/1‑149 (MQ=255)                                                             ccAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTc                                                                                                                                   >  7:40193/1‑92 (MQ=255)                                                             ccAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTc                                                                                                                                   <  8:40193/92‑1 (MQ=255)                                                                   tGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGccc                                                                    <  6:113083/149‑1 (MQ=255)                                                                      cgcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAtt                                                                 <  1:92778/149‑1 (MQ=255)                                                                       gcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTc                                                                                                                                   >  7:114144/1‑82 (MQ=255)                                                                       gcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTc                                                                                                                                   <  8:114144/82‑1 (MQ=255)                                                                       gcTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAttt                                                                <  1:86775/149‑1 (MQ=255)                                                                                          ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                      >  5:10941/1‑140 (MQ=255)                                                                                          ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                      <  6:10941/140‑1 (MQ=255)                                                                                              aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg                                              >  7:92296/1‑144 (MQ=255)                                                                                              aTCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg                                              <  8:92296/144‑1 (MQ=255)                                                                                                  tGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTACAACCGCAACCTTCCGTTTCTAATTTCTCtttt                                                                                         >  7:59584/1‑97 (MQ=255)                                                                                                  tGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTACAACCGCAACCTTCCGTTTCTAATTTCTCtttt                                                                                         <  8:59584/97‑1 (MQ=255)                                                                                                         gTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGtt                                                 >  2:12777/1‑130 (MQ=255)                                                                                                         gTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGtt                                                 <  1:12777/130‑1 (MQ=255)                                                                                                         gTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAAt                               >  3:92335/1‑148 (MQ=255)                                                                                                              cAGCTGCTTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTcc                         >  5:54693/1‑149 (MQ=255)                                                                                                                gctgctTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCg                        >  1:178674/1‑148 (MQ=255)                                                                                                                  tgctTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACaa                                                           <  6:47861/111‑1 (MQ=255)                                                                                                                  tgctTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACaa                                                           >  5:47861/1‑111 (MQ=255)                                                                                                                  tgctTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATa                     >  3:60899/1‑149 (MQ=255)                                                                                                                    ctTTAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATa                     <  2:178674/147‑1 (MQ=255)                                                                                                                      ttAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                      <  2:72430/112‑1 (MQ=255)                                                                                                                      ttAACACTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg                                                      >  1:72430/1‑112 (MQ=255)                                                                                                                           ccgagctgctgctTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACtttt             <  7:53021/141‑1 (MQ=255)                                                                                                                           acTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATaca                                    <  4:107948/125‑1 (MQ=255)                                                                                                                           acTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATaca                                    >  3:107948/1‑125 (MQ=255)                                                                                                                           acTTATGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTg            >  6:155039/1‑149 (MQ=255)                                                                                                                                tGGCTGCTTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACt        <  4:60899/148‑1 (MQ=255)                                                                                                                                   ctgctTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACCCGAATCATTCCGATACCACtttt             >  8:53021/1‑140 (MQ=255)                                                                                                                                      ctTCTGGGGC‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTaa  >  2:137073/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                           GTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAA  >  NC_000913/169543‑169824

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