Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 167,843 +T coding (360/2244 nt) fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,8421.T94.4% 116.4 / 0.3 36coding (359/2244 nt)fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (15/19);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

CCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATA  >  NC_000913/167698‑167955
                                                                                                                                                 |                                                                                                                 
ccGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAttga                                                                                                               >  3:117898/1‑146 (MQ=255)
               aaaaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTcg                                                                                                 >  3:50210/1‑149 (MQ=255)
                 aaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAtt                                                                                                                 >  7:22819/1‑129 (MQ=255)
                 aaaGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTAtt                                                                                                                 <  8:22819/130‑2 (MQ=255)
                          caGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCgg                                                                                      <  2:1934/149‑1 (MQ=255)
                           aGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGc                                                                                     <  5:129835/149‑1 (MQ=255)
                                       ttgtGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTATAAGAAGCGCTTAGCTACACGCCGGGTGTGTCTGTTCGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCaaa                                                                          <  2:126634/148‑1 (MQ=255)
                                          tgACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCa                                                                      <  4:50210/149‑1 (MQ=255)
                                              cgccgAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGaa                                                                   >  6:15633/1‑148 (MQ=255)
                                                   aaGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAgg                                                                               >  4:95330/1‑131 (MQ=255)
                                                   aaGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAgg                                                                               <  3:95330/131‑1 (MQ=255)
                                                     gagaTGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATc                                                           <  5:15633/149‑1 (MQ=255)
                                                      agaTGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCt                                                          <  8:17221/149‑1 (MQ=255)
                                                           gcgcTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGaccac                                                                                                            >  3:90299/1‑94 (MQ=255)
                                                           gcgcTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGaccac                                                                                                            <  4:90299/94‑1 (MQ=255)
                                                                 cATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTg                                                >  3:43347/1‑148 (MQ=255)
                                                                          aaGTCGGTAAAAGAAGCGCTTAGCTACATGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGc                                                                                        <  2:32545/99‑1 (MQ=255)
                                                                          aaGTCGGTAAAAGAAGCGCTTAGCTACATGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGc                                                                                        >  1:32545/1‑99 (MQ=255)
                                                                               ggTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAAc                                 >  6:72829/1‑149 (MQ=255)
                                                                                 tAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACtt                               <  1:113104/149‑1 (MQ=255)
                                                                                    aagaagCGCTTAGCTACACGCCGGGTGGCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCCCTGCAACTTcag                            <  5:121367/149‑3 (MQ=255)
                                                                                     agaagCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCAAACAACTATTTACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTGCta                            >  7:56043/1‑148 (MQ=255)
                                                                                       aagCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATaa                         <  4:117898/149‑1 (MQ=255)
                                                                                                  tACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTca                             >  2:52700/1‑133 (MQ=255)
                                                                                                  tACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTca                             <  1:52700/134‑2 (MQ=255)
                                                                                                     acGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCAtt            >  6:2815/1‑148 (MQ=255)
                                                                                                         cGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAccc       <  5:2815/149‑1 (MQ=255)
                                                                                                               tctcTGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGtata  <  3:51456/148‑1 (MQ=255)
                                                                                                                 tctGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGcc                                                                             >  1:151323/1‑71 (MQ=255)
                                                                                                                 tctGTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGcc                                                                             <  2:151323/71‑1 (MQ=255)
                                                                                                                 tctGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTAt                                                            <  4:29109/87‑1 (MQ=255)
                                                                                                                 tctGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTAt                                                            >  3:29109/1‑87 (MQ=255)
                                                                                                                    gTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCa                                                                      >  7:77472/1‑75 (MQ=255)
                                                                                                                    gTTGGTACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCa                                                                      <  8:77472/75‑1 (MQ=255)
                                                                                                                         tACGCGTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCta                            >  2:157770/1‑112 (MQ=255)
                                                                                                                             cgTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCa                                        <  6:138152/96‑1 (MQ=255)
                                                                                                                             cgTGGCGCATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCa                                        >  5:138152/1‑96 (MQ=255)
                                                                                                                                 gcgcATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGc                                         >  6:129442/1‑91 (MQ=255)
                                                                                                                                 gcgcATCCAACACCTATTGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGc                                         <  5:129442/91‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                 
CCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATA  >  NC_000913/167698‑167955

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