Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,541,962 T→G I519S (ATC→AGC)  feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,541,9620TG90.7% 249.9 / 14.6 75I519S (ATC→AGC) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (35/33);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

TGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTC  >  NC_000913/3541866‑3542056
                                                                                                |                                                                                              
tGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcat                                                                                             <  2:1530528/100‑1 (MQ=255)
tGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcat                                                                                             >  1:1734311/1‑100 (MQ=255)
    gTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcatcatc                                                                                         >  1:866171/1‑100 (MQ=255)
     ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCg                                                                                        <  1:831354/100‑1 (MQ=255)
     ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCg                                                                                        <  1:478042/100‑1 (MQ=255)
      tGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGt                                                                                       >  2:762670/1‑100 (MQ=255)
       gCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTc                                                                                      <  1:1007274/100‑1 (MQ=255)
       gCCGGTCGATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTc                                                                                      >  2:1023079/1‑100 (MQ=255)
        ccGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCa                                                                                     <  2:876108/100‑1 (MQ=255)
         cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAg                                                                                    <  1:1772985/100‑1 (MQ=255)
         cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAg                                                                                    <  1:1705726/100‑1 (MQ=255)
         cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTCCAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAg                                                                                    >  1:14282/1‑100 (MQ=255)
            tCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCAt                                                                                 <  1:636439/100‑1 (MQ=255)
             cTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAtt                                                                                >  2:437714/1‑100 (MQ=255)
               ntCATGTACCACACGTTAAAAGACTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTCCGTGCTGGTAAAGTGGGCATCATCGTCAGCAtttt                                                                              >  1:909783/2‑100 (MQ=255)
               aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAtttt                                                                              >  2:524339/1‑100 (MQ=255)
               aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAtttt                                                                              <  2:845818/100‑1 (MQ=255)
               aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAtttt                                                                              >  2:729321/1‑100 (MQ=255)
                  aTGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCt                                                                           >  2:1192601/1‑100 (MQ=255)
                   tGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTg                                                                          <  1:119972/100‑1 (MQ=255)
                     tACCAGACGTTAAAGGCATCTTTATCTTGACTTGGCAGCGTCTGAAAGGCTTCGTTCTGTGTGCTGGTAAAGTGAGCATCATCTTCCGCATTTGCCTgn                                                                         <  1:1159449/99‑2 (MQ=255)
                      aCCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgc                                                                       >  1:267031/1‑100 (MQ=255)
                       ccACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcg                                                                      >  1:12679/1‑100 (MQ=255)
                        cacaCGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcgc                                                                     >  2:865770/1‑100 (MQ=255)
                         acacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCt                                                                    >  1:790226/1‑100 (MQ=255)
                         acacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCt                                                                    <  1:1777496/100‑1 (MQ=255)
                          cacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCtt                                                                   >  1:1284797/1‑100 (MQ=255)
                            cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                 >  2:243881/1‑100 (MQ=255)
                            cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                 <  1:92867/100‑1 (MQ=255)
                            cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                 <  1:14247/100‑1 (MQ=255)
                            cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                 <  1:1580062/100‑1 (MQ=255)
                             gTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCa                                                                <  2:1137915/100‑1 (MQ=255)
                             gTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCa                                                                <  1:1811767/100‑1 (MQ=255)
                              ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCaa                                                               >  2:1252735/1‑100 (MQ=255)
                              ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCaa                                                               >  1:1706903/1‑100 (MQ=255)
                               tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAAc                                                              <  2:1619851/100‑1 (MQ=255)
                                 aaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAg                                                            >  1:339608/1‑100 (MQ=255)
                                 aaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAg                                                            <  1:1221274/100‑1 (MQ=255)
                                  aaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGc                                                           <  1:60446/100‑1 (MQ=255)
                                    gCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCtt                                                         >  1:1385995/1‑100 (MQ=255)
                                       tGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTtctc                                                      <  1:566512/100‑1 (MQ=255)
                                        gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCg                                                     >  2:1615492/1‑100 (MQ=255)
                                         aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGc                                                    <  1:602453/100‑1 (MQ=255)
                                          ttATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCt                                                   >  2:1004088/1‑100 (MQ=255)
                                            aTCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                 >  2:1048500/1‑100 (MQ=255)
                                               cAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCg                                              <  1:1458299/100‑1 (MQ=255)
                                                aGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCgg                                             <  2:163405/100‑1 (MQ=255)
                                                 gACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCggg                                            <  1:864789/100‑1 (MQ=255)
                                                   ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                          >  2:1839265/1‑100 (MQ=255)
                                                   ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                          >  2:425783/1‑100 (MQ=255)
                                                       gCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAt                                       <  1:1846701/99‑1 (MQ=255)
                                                        cAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcg                                     >  1:1172729/1‑100 (MQ=255)
                                                        cAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcg                                     <  2:35058/100‑1 (MQ=255)
                                                         aGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgt                                    >  2:437386/1‑100 (MQ=255)
                                                          gCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtc                                   <  1:1107654/100‑1 (MQ=255)
                                                           cGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtcg                                  <  2:866172/100‑1 (MQ=255)
                                                           cGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtcg                                  <  1:1176217/100‑1 (MQ=255)
                                                            gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                 >  2:1302878/1‑100 (MQ=255)
                                                             tCTGAAAGGCTTCGTTCTGCGTGCTCGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGAt                                <  2:530400/100‑1 (MQ=255)
                                                              cTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                               >  2:1801779/1‑100 (MQ=255)
                                                               tGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATaa                              <  1:1252731/100‑1 (MQ=255)
                                                               tGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATaa                              >  1:1447672/1‑100 (MQ=255)
                                                                gAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAc                             >  1:1524001/1‑100 (MQ=255)
                                                                   aGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATc                          <  2:209751/100‑1 (MQ=255)
                                                                              tGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgc               <  2:1677492/100‑1 (MQ=255)
                                                                               gCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcg              >  1:961818/1‑100 (MQ=255)
                                                                               gCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcg              >  2:1684386/1‑100 (MQ=255)
                                                                               gCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcg              >  1:961811/1‑100 (MQ=255)
                                                                                  tGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcgc             >  1:325610/1‑98 (MQ=255)
                                                                                   gCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg          <  1:1801768/100‑1 (MQ=255)
                                                                                   gCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg          <  1:76083/100‑1 (MQ=255)
                                                                                   gCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg          <  1:762669/100‑1 (MQ=255)
                                                                                      ggTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt       >  2:1939626/1‑100 (MQ=255)
                                                                                      ggTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt       >  2:1939478/1‑100 (MQ=255)
                                                                                         aaaGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg    >  1:1619003/1‑100 (MQ=255)
                                                                                         aaaGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg    >  1:507750/1‑100 (MQ=255)
                                                                                         aaaGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg    <  2:1447680/100‑1 (MQ=255)
                                                                                           aGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc  >  2:861127/1‑100 (MQ=255)
                                                                                                |                                                                                              
TGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTC  >  NC_000913/3541866‑3542056

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