Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 169,645 A→G Y721C (TAC→TGC)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913169,6450AG100.0% 117.0 / NA 33Y721C (TAC→TGC) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/16);  total (17/16)

ATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTC  >  NC_000913/169516‑169763
                                                                                                                                 |                                                                                                                      
aTCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                               >  1:193988/1‑139 (MQ=255)
aTCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                               <  1:543607/139‑1 (MQ=255)
                              ataCGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTgggg                                                                                 <  2:248783/139‑1 (MQ=255)
                                   gTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTg                                                                                    >  2:923589/1‑131 (MQ=255)
                                   gTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTg                                                                                    <  1:923589/131‑1 (MQ=255)
                                              gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                            >  2:476345/1‑112 (MQ=255)
                                              gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                            <  1:476345/112‑1 (MQ=255)
                                               ttAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTa                                                                                                >  1:469040/1‑107 (MQ=255)
                                               ttAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTa                                                                                                <  2:469040/107‑1 (MQ=255)
                                                       ttATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaacc                                                        >  2:127437/1‑139 (MQ=255)
                                                               tGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCtt                                                <  2:862833/139‑1 (MQ=255)
                                                                ggCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCtcc                                               >  1:330556/1‑137 (MQ=255)
                                                                                      tCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                            <  2:886896/72‑1 (MQ=255)
                                                                                      tCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                            >  1:886896/1‑72 (MQ=255)
                                                                                           cGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                               >  2:75682/1‑48 (MQ=255)
                                                                                           cGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                               <  1:75682/48‑1 (MQ=255)
                                                                                               gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAgg                                                                   >  1:855547/1‑88 (MQ=255)
                                                                                               gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAgg                                                                   <  2:855547/88‑1 (MQ=255)
                                                                                                  cTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGc                       <  2:19235/129‑1 (MQ=255)
                                                                                                  cTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGc                       >  1:19235/1‑129 (MQ=255)
                                                                                                     cATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTaa                                      >  1:510653/1‑111 (MQ=255)
                                                                                                     cATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTaa                                      <  2:510653/111‑1 (MQ=255)
                                                                                                             caacCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTc  >  1:428418/1‑139 (MQ=255)
                                                                                                                  tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgc                                                            <  1:866971/76‑1 (MQ=255)
                                                                                                                  tGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgc                                                            >  2:866971/1‑76 (MQ=255)
                                                                                                                   gTTCGATCGTGAATGCGTCGTCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTcc             >  1:631118/1‑122 (MQ=255)
                                                                                                                   gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                            <  2:227327/43‑1 (MQ=255)
                                                                                                                   gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                            >  1:227327/1‑43 (MQ=255)
                                                                                                                   gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAGTTCTCTTTTGGGGCACGGATTTcc             <  2:631118/122‑1 (MQ=255)
                                                                                                                    ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCt                                        >  2:384599/1‑94 (MQ=255)
                                                                                                                    ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCt                                        <  1:384599/94‑1 (MQ=255)
                                                                                                                            tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCttt                             >  2:44357/1‑97 (MQ=255)
                                                                                                                            tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCttt                             <  1:44357/97‑1 (MQ=255)
                                                                                                                                 |                                                                                                                      
ATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTC  >  NC_000913/169516‑169763

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