Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 170,457 T→A 6.8% I227N (ATT→AAT)  fhuC → iron(3+)‑hydroxamate import ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913170,4570TA 6.8% 184.6 / 4.0 59I227N (ATT→AAT) fhuCiron(3+)‑hydroxamate import ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base A (2/2);  total (29/30)
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.50e-01

CTGGATATCGCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCGGCTTACCTC  >  NC_000913/170309‑170590
                                                                                                                                                    |                                                                                                                                     
cTGGATATCGCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGAt                                                                                                                                       <  1:195155/149‑1 (MQ=255)
   gATATCGCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGc                                                                                                                                    >  3:302950/1‑149 (MQ=255)
       tCGCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGTGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAg                                                                                                                                >  3:147279/1‑149 (MQ=255)
       tCGCCCACCAGGTTGATGTGCTGTCGCTGGGGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGGGTTGCACGATAGCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCa                                                                                                                                 <  5:70873/148‑1 (MQ=255)
         gCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAggg                                                                                                                              <  7:206795/149‑1 (MQ=255)
         gCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAgg                                                                                                                               >  3:295929/1‑148 (MQ=255)
         gCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAg                                                                                                                                >  3:274869/1‑147 (MQ=255)
         gCCCACCAGGTTGATGTGCTGTCGCTGGGGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGGTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAg                                                                                                                                <  4:274869/147‑1 (MQ=255)
                  gTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGc                                                                                                                     >  2:142109/1‑149 (MQ=255)
                              tCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATgc                                                                                                          <  3:32261/148‑1 (MQ=255)
                                gCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGcgcg                                                                                                       >  3:254846/1‑149 (MQ=255)
                                    gTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCgcggc                                                                                                     <  8:373064/147‑1 (MQ=255)
                                      gCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGaaa                                                                                                 >  4:233282/1‑149 (MQ=255)
                                       cACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAAc                                                                                                <  4:295929/149‑1 (MQ=255)
                                                 gTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAAt                                                                                       >  7:64620/1‑148 (MQ=255)
                                                 gTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAAt                                                                                       >  8:249399/1‑148 (MQ=255)
                                                    aGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGAtt                                                                                   <  1:142109/149‑1 (MQ=255)
                                                     ggAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTc                                                                                            >  4:253791/1‑139 (MQ=255)
                                                     ggAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTc                                                                                            <  3:253791/139‑1 (MQ=255)
                                                     ggAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGAttt                                                                                  <  4:254846/149‑1 (MQ=255)
                                                       aGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTAt                                                                                <  4:79025/149‑1 (MQ=255)
                                                        gCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATg                                                                               >  7:155867/1‑149 (MQ=255)
                                                             gCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCgcggc                                                                                                     >  1:160351/1‑122 (MQ=255)
                                                             gCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCgcggc                                                                                                     <  2:160351/122‑1 (MQ=255)
                                                             cccTGACGGTCATTGCCGTGTTGCACGATAACAATATGTCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATc                                                                          >  8:277222/2‑149 (MQ=255)
                                                                       cATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTa                                                                >  2:54915/1‑149 (MQ=255)
                                                                        aTTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTa                                                                <  7:249399/148‑1 (MQ=255)
                                                                            ccGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTAtttt                                                            <  1:37996/148‑1 (MQ=255)
                                                                             cGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGc                                                          <  8:155867/149‑1 (MQ=255)
                                                                             cGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGc                                                          >  5:84701/1‑149 (MQ=255)
                                                                                       gATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCcggc                                                >  1:124451/1‑149 (MQ=255)
                                                                                                  ggcggcACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCAc                                     >  8:29295/1‑149 (MQ=255)
                                                                                                   gcggcACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCAc                                     <  2:124451/148‑1 (MQ=255)
                                                                                                      gcACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCtgt                                 <  6:84701/149‑1 (MQ=255)
                                                                                                        aCGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATgg                                                                   <  4:220191/113‑1 (MQ=255)
                                                                                                        aCGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATgg                                                                   >  3:220191/1‑113 (MQ=255)
                                                                                                           cTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGtt                            <  1:231301/149‑1 (MQ=255)
                                                                                                                gtgATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGtt                       >  4:28023/1‑149 (MQ=255)
                                                                                                                 tgATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGtt                       <  7:277222/148‑1 (MQ=255)
                                                                                                                 tgATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGtt                       >  6:175981/1‑148 (MQ=255)
                                                                                                                    ttATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTAtt                   >  3:134900/1‑149 (MQ=255)
                                                                                                                    ttATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTAtt                   <  7:29295/149‑1 (MQ=255)
                                                                                                                    ttATCTGGTCGCCCTGCGCGGCGGTGAAATGAATGCTCAGGGAACGCCt                                                                                                                       >  2:179305/1‑49 (MQ=255)
                                                                                                                    ttATCTGGTCGCCCTGCGCGGCGGTGAAATGAATGCTCAGGGAACGCCt                                                                                                                       >  2:184656/1‑49 (MQ=255)
                                                                                                                    ttATCTGGTCGCCCTGCGCGGCGGTGAAATGAATGCTCAGGGAACGCCt                                                                                                                       <  1:179305/49‑1 (MQ=255)
                                                                                                                    ttATCTGGTCGCCCTGCGCGGCGGTGAAATGAATGCTCAGGGAACGCCt                                                                                                                       <  1:184656/49‑1 (MQ=255)
                                                                                                                       tCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCAc                                     <  2:209401/128‑1 (MQ=255)
                                                                                                                       tCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCAc                                     >  1:209401/1‑128 (MQ=255)
                                                                                                                       tCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTa                     <  6:329387/144‑1 (MQ=255)
                                                                                                                       tCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTa                     >  5:329387/1‑144 (MQ=255)
                                                                                                                           gTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATtgat                <  4:134900/145‑1 (MQ=255)
                                                                                                                            tCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCg           <  1:54915/149‑1 (MQ=255)
                                                                                                                                  tGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCGGCTTAc     >  3:309659/1‑149 (MQ=255)
                                                                                                                                  tGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCGGCTTAc     >  3:304449/1‑149 (MQ=255)
                                                                                                                                     gcggcgGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGTGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCGGCTTACctc  <  4:302950/149‑1 (MQ=255)
                                                                                                                                     gcggcgGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCGGCTTACctc  <  4:147279/149‑1 (MQ=255)
                                                                                                                                          ggTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGt                        <  3:138155/122‑1 (MQ=255)
                                                                                                                                          ggTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGt                        >  4:138155/1‑122 (MQ=255)
                                                                                                                                             gAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGtg                                           >  8:199585/1‑100 (MQ=255)
                                                                                                                                             gAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGtg                                           <  7:199585/100‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
CTGGATATCGCCCACCAGGTTGATGTGCTGTCGCTGGTGCACCGTTTAAGTCAGGAGCGTGGCCTGACGGTCATTGCCGTGTTGCACGATATCAATATGGCGGCACGCTACTGTGATTATCTGGTCGCCCTGCGCGGCGGTGAAATGATTGCTCAGGGAACGCCTGCGGAAATTATGCGCGGCGAAACCCTCGAAATGATTTATGGCATCCCGATGGGTATTTTGCCGCATCCGGCGGGTGCTGCACCTGTGAGTTTTGTTTATTGATGAGCGGCTTACCTC  >  NC_000913/170309‑170590

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: