Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 168,066 A→T 76.7% K195* (AAA→TAA)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913168,0660AT76.7% 101.4 / 33.7 56K195* (AAA→TAA) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base T (21/22);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01

GCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGT  >  NC_000913/167920‑168204
                                                                                                                                                  |                                                                                                                                          
gCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGaa                                                                                                                                           <  4:37179/148‑1 (MQ=255)
    cTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAaagaa                                                                                                                                       <  5:125708/148‑1 (MQ=255)
      tCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGt                                                                                                                                     >  3:89234/1‑148 (MQ=255)
       cTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGt                                                                                                                                     >  1:232573/1‑147 (MQ=255)
          taACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTc                                                                                                                                   <  2:56304/146‑1 (MQ=255)
           aaCGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGt                                                                                                                                >  8:234853/1‑148 (MQ=255)
            aCGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGt                                                                                                                                <  4:278695/147‑1 (MQ=255)
               aTGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTaa                                                                                                                            <  5:187922/148‑1 (MQ=255)
                tGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTa                                                                                                                             >  1:309061/1‑146 (MQ=255)
                tGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTa                                                                                                                             <  2:309061/146‑1 (MQ=255)
                  cgggctttGACCCGTATATGCTGGATCGCGCTGAAATTATGCGTGGCCCGTTTTCCGTGCTTTACGGTAAAAGCAGGCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAg                                                                                                                          <  6:104422/141‑1 (MQ=255)
                          gACCCGTATATGCTGGAACGCGCTGAAATTATGCGTTGCCCGGTTTCCGTGCGTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCTAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACt                                                                                                                  >  2:221567/1‑147 (MQ=255)
                              tgTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAg                                                                                                             >  3:21815/2‑148 (MQ=255)
                                 atatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCt                                                                                                          >  6:233549/1‑148 (MQ=255)
                                 atatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCt                                                                                                          >  5:232475/1‑148 (MQ=255)
                                 atatGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCt                                                                                                          >  3:99166/1‑148 (MQ=255)
                                      cTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTc                                                                                                      <  4:99166/147‑1 (MQ=255)
                                      cTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGGCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                                                                                     <  2:289527/148‑1 (MQ=255)
                                      cTGGAACGCGCTGAAATTATGCGTGGCCCGGATTCCGTGCTTTACGGTAAAAGCAGTCCTGCCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                                                                                     <  2:285034/148‑1 (MQ=255)
                                                   aaaTTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTg                                                                                         <  5:233549/147‑1 (MQ=255)
                                                   aaaTTATGCGTGGCACGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGCTTTTGa                                                                                        >  2:37241/1‑148 (MQ=255)
                                                      ttATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGGCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAGGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACtt                                                                                     <  8:124041/148‑1 (MQ=255)
                                                         tGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTa                                                                                   <  6:232475/147‑1 (MQ=255)
                                                             tGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGttt                                                                                           <  3:140303/135‑1 (MQ=255)
                                                             tGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGttt                                                                                           <  3:145107/135‑1 (MQ=255)
                                                             tGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGttt                                                                                           >  4:140303/1‑135 (MQ=255)
                                                             tGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGttt                                                                                           >  4:145107/1‑135 (MQ=255)
                                                             tGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGAt                                                                              <  7:234853/148‑1 (MQ=255)
                                                                 ccGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGt                                                                          <  4:89234/148‑1 (MQ=255)
                                                                   ggTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTg                                                                        <  4:21815/148‑1 (MQ=255)
                                                                     tttCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACt                                                                                      <  7:162593/132‑1 (MQ=255)
                                                                     tttCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACAGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACt                                                                                      >  8:162593/1‑132 (MQ=255)
                                                                       tCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGgat                                                                     >  1:149970/1‑147 (MQ=255)
                                                                           tGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTc                                                                            <  7:42979/136‑1 (MQ=255)
                                                                           tGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTc                                                                            >  8:42979/1‑136 (MQ=255)
                                                                             cTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCCCCGAACCGCTGTAAGAAGTTCCGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGAc                                                               <  2:239472/147‑1 (MQ=255)
                                                                             cTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACg                                                              >  4:76982/1‑148 (MQ=255)
                                                                                   ggTAAAAGCAGTCCTGTCGGCCTGTTGAATATGGTCAGCAAGCGTCAGACCAACGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGGGttt                                                        <  1:221567/148‑1 (MQ=255)
                                                                                       aaaGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTAct                                                     <  2:149970/147‑1 (MQ=255)
                                                                                               ccTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                                                                                     >  5:210515/1‑91 (MQ=255)
                                                                                               ccTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTcc                                                                                                     <  6:210515/91‑1 (MQ=255)
                                                                                                          gTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCtt                                 >  5:188473/1‑148 (MQ=255)
                                                                                                          gTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCtt                                 >  7:259548/1‑148 (MQ=255)
                                                                                                          gTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCtt                                 >  3:67156/1‑148 (MQ=255)
                                                                                                           ttGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTg                                <  2:232573/148‑1 (MQ=255)
                                                                                                              aaTATGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTgcgc                             >  7:191514/1‑148 (MQ=255)
                                                                                                                 atGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGtt                          >  2:297829/1‑148 (MQ=255)
                                                                                                                  tGGTCAGCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTc                         <  4:67156/148‑1 (MQ=255)
                                                                                                                  tGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGtt                          <  1:297829/147‑1 (MQ=255)
                                                                                                                   ggTCAGCAAGCGTCCGACCACCGACCCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGAGTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGCCCGGTATTGCGCGTTCt                        >  2:102348/1‑148 (MQ=255)
                                                                                                                        gCAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCaa                   >  2:36128/1‑148 (MQ=255)
                                                                                                                         cAAGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAAt                  >  5:195194/1‑148 (MQ=255)
                                                                                                                           aGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGttt                                                        >  2:175386/1‑108 (MQ=255)
                                                                                                                           aGCGTCCGACCACCGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGttt                                                        <  1:175386/108‑1 (MQ=255)
                                                                                                                                   accaccGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCcagcag         >  6:18766/1‑147 (MQ=255)
                                                                                                                                        cGAACCGCTGTAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAggg   <  3:76982/148‑1 (MQ=255)
                                                                                                                                         gAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCAAATGCCCAGCAGAAAGGGt  <  1:37241/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACGGTGTTTACTCTTATCGCCTGACCGGTCTTGCGCGTTCTGCCAATGCCCAGCAGAAAGGGT  >  NC_000913/167920‑168204

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