Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 722,155 A→T 11.1% K42N (AAA→AAT fhuC → iron (Fe3+) ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730722,1550AT11.1% 101.6 / 7.2 36K42N (AAA→AATfhuCiron (Fe3+) ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base A (14/18);  new base T (2/2);  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.41e-01

TTTTATTGAAAGTAGGGAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGC  >  CP000730/722008‑722302
                                                                                                                                                   |                                                                                                                                                   
ttttATTGAAAGTAGGGAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAAt                                                                                                                                                    <  2:190318/149‑1 (MQ=255)
 tttATTGAAAGTAGGGAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATc                                                                                                                                                   >  1:167736/1‑149 (MQ=255)
  ttATTGAAAGTAGGGAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCt                                                                                                                                                  <  1:267134/149‑1 (MQ=255)
      tGAAAGTAGGGAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACtt                                                                                                                                              <  2:92190/149‑1 (MQ=255)
             aGGTAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTaaa                                                                                                                                       <  1:152947/149‑1 (MQ=255)
               ggAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAATTCTACTTTGCTAAAgg                                                                                                                                     >  1:130833/1‑149 (MQ=255)
               ggAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAATTCTACTTTGCTAAAgg                                                                                                                                     >  1:130837/1‑149 (MQ=255)
               ggAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAgg                                                                                                                                     >  2:146629/1‑149 (MQ=255)
               ggAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAgg                                                                                                                                     <  2:111728/149‑1 (MQ=255)
               ggAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAgg                                                                                                                                     >  1:210713/1‑149 (MQ=255)
                     tATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGt                                                                                                                               >  1:167208/1‑149 (MQ=255)
                          aTCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGAGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGt                                                                                                                          >  1:90754/1‑149 (MQ=255)
                                    tGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAg                                                                                                                <  1:48174/149‑1 (MQ=255)
                                    tGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAg                                                                                                                <  1:54557/149‑1 (MQ=255)
                                     ggACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGt                                                                                                               <  2:147792/149‑1 (MQ=255)
                                      gACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGtt                                                                                                              >  2:222737/1‑149 (MQ=255)
                                                   aaTTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAg                                                                                                 >  2:111081/1‑149 (MQ=255)
                                                                          ttATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATAtt                                                                          <  1:143876/149‑1 (MQ=255)
                                                                              aaataaatTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATa                                                                      <  2:202003/149‑1 (MQ=255)
                                                                                 taaatTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATacac                                                                   >  2:42765/1‑149 (MQ=255)
                                                                                  aaatGAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTACCGGCGGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATAcaca                                                                  <  1:48786/149‑1 (MQ=255)
                                                                                                aaTACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAgaga                                                    >  2:261074/1‑149 (MQ=255)
                                                                                                aaTACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAgaga                                                    >  1:129255/1‑149 (MQ=255)
                                                                                                   aCCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTg                                                 <  2:136461/149‑1 (MQ=255)
                                                                                                       gATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCaaa                                             >  2:151977/1‑149 (MQ=255)
                                                                                                             aaaGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAAt                                       <  1:84008/149‑1 (MQ=255)
                                                                                                                         aTCATTTGTCCTAACGGCTGCGGGAAAGCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTAcc                           <  2:8330/149‑1 (MQ=255)
                                                                                                                         aTCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTAcc                           >  1:123736/1‑149 (MQ=255)
                                                                                                                                 tCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCAc                   >  1:78595/1‑149 (MQ=255)
                                                                                                                                     aaCGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGa               <  1:78015/149‑1 (MQ=255)
                                                                                                                                      aCGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGaa              <  2:256839/149‑1 (MQ=255)
                                                                                                                                        ggCTGCGGGAATTCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGt            <  2:130833/149‑1 (MQ=255)
                                                                                                                                        ggCTGCGGGAATTCTACTTTGCTAAAGGCATTGTCACGTTTATTGGAAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAAAAAAAAAAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGt            <  2:130837/149‑1 (MQ=255)
                                                                                                                                        ggCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGAAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGt            <  2:210713/149‑1 (MQ=255)
                                                                                                                                        ggCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGGGTAGTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATATAAGAAAAAAAAAAAAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGt            <  1:146629/149‑1 (MQ=255)
                                                                                                                                              gggAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGa      >  2:47082/1‑149 (MQ=255)
                                                                                                                                                  aaTCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAAATAAAAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGc  <  2:141513/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTTATTGAAAGTAGGGAAATTATGAATCGTTTGCATGGACAACAAGTTAAAATTGGTTACGGGGATAACACGATTATAAATAAATTAGATGTTGAAATACCAGATGGCAAAGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGC  >  CP000730/722008‑722302

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