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.1 / NA 33Y721C (TAC→TGC) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/19);  total (14/19)

ATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATG  >  NC_000913/169516‑169780
                                                                                                                                 |                                                                                                                                       
aTCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                                >  2:84781/1‑139 (MQ=255)
aTCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                                <  1:576557/139‑1 (MQ=255)
 tCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                               <  2:460575/139‑1 (MQ=255)
 tCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                               <  2:80237/139‑1 (MQ=255)
        aaCTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                        <  1:296615/139‑1 (MQ=255)
              tttAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACtt                                                                                                                  <  1:469286/139‑1 (MQ=255)
                   aGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                             <  2:552022/139‑1 (MQ=255)
                       ggAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctgct                                                                                                         >  2:29064/1‑139 (MQ=255)
                            ttATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCAAGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTgg                                                                                                    <  2:790369/139‑1 (MQ=255)
                                aCGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcg                                                                                                <  1:151725/139‑1 (MQ=255)
                                   gTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTTGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAg                                                                                             >  1:498879/1‑139 (MQ=255)
                                      gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGACGCAGAAc                                                                                          >  1:351552/1‑139 (MQ=255)
                                                       ttATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaacc                                                                         >  2:667127/1‑139 (MQ=255)
                                                        tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                        >  2:608419/1‑139 (MQ=255)
                                                          tGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccgca                                                                      >  1:255884/1‑139 (MQ=255)
                                                                 gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCAGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                        >  1:494244/1‑82 (MQ=255)
                                                                 gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGGTAACAACCAGTTCGATCGTGAATGCGTCGCCAGCTGCtttt                                                                                                                        <  2:494244/82‑2 (MQ=255)
                                                                  cgcgAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCg                                                              >  2:138964/1‑139 (MQ=255)
                                                                      aGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTc                                                          >  1:544799/1‑139 (MQ=255)
                                                                       gTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGt                                                                                        <  1:230173/108‑1 (MQ=255)
                                                                       gTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGt                                                                                        >  2:230173/1‑108 (MQ=255)
                                                                            cATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAttt                                                    >  2:1082/1‑139 (MQ=255)
                                                                            cATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAttt                                                    <  2:566783/139‑1 (MQ=255)
                                                                                             tGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGa                                   <  2:244121/139‑1 (MQ=255)
                                                                                              ggCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTgggg                                                                                                  <  1:670173/75‑1 (MQ=255)
                                                                                              ggCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTgggg                                                                                                  >  2:670173/1‑75 (MQ=255)
                                                                                               gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGAtt                                 >  1:631367/1‑139 (MQ=255)
                                                                                                     cATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTg                           <  1:29064/139‑1 (MQ=255)
                                                                                                       tGTTAACAACCTGTTCGATCGTGAATGCGTAGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTggc                         <  1:482709/139‑3 (MQ=255)
                                                                                                              aacCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCCGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTGTGGGGCACGGATTTCCGTGCCCATTTca                  <  2:498879/139‑1 (MQ=255)
                                                                                                                         tCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTg       <  2:544799/139‑1 (MQ=255)
                                                                                                                         tCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTCCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTg       <  1:84781/139‑1 (MQ=255)
                                                                                                                              aaTGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg  <  1:608419/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                       
ATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATG  >  NC_000913/169516‑169780

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: ‑