Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,542,093 G→A V563M (GTG→ATG)  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,542,0930GA98.4% 207.0 / ‑4.1 63V563M (GTG→ATG) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (26/36);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGC  >  NC_000913/3541996‑3542187
                                                                                                 |                                                                                              
cAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTg                                                                                              >  1:153321/1‑100 (MQ=255)
  gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                            <  1:430887/100‑1 (MQ=255)
  gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                            <  1:1666824/100‑1 (MQ=255)
  gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGCCGCTGGCCTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGACATGCa                                                                                            >  2:325977/1‑100 (MQ=255)
    ttCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCTCCCCGTTCTTCAAGCCAATTGGcatgcatg                                                                                          <  1:1082398/100‑1 (MQ=255)
    ttCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGcatgcatg                                                                                          <  2:1180406/100‑1 (MQ=255)
    ttCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGcatgcatg                                                                                          <  2:1513393/100‑1 (MQ=255)
     tctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGa                                                                                         >  2:1308087/1‑100 (MQ=255)
     tctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGa                                                                                         <  1:1291755/100‑1 (MQ=255)
        cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGa                                                                                      <  1:641149/100‑1 (MQ=255)
        cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGa                                                                                      <  1:803165/100‑1 (MQ=255)
              gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCTATTGGCATGCATGAAGATAACt                                                                                 <  1:1442712/99‑1 (MQ=255)
              gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTg                                                                                <  2:181653/100‑1 (MQ=255)
              gCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTg                                                                                <  2:179514/100‑1 (MQ=255)
                gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTagc                                                                              <  2:1625768/100‑1 (MQ=255)
                  gAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcag                                                                            >  1:1108541/1‑100 (MQ=255)
                  gAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcag                                                                            >  1:1036524/1‑100 (MQ=255)
                    aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGGCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGTATGAAGATAACTggcaggc                                                                          <  1:1164076/100‑1 (MQ=255)
                    aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACaggcaggc                                                                          >  2:1759003/1‑100 (MQ=255)
                     aaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcaggca                                                                         >  2:1438834/1‑100 (MQ=255)
                            gATAACATCAACGGCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTg                                                                  <  2:585723/100‑1 (MQ=255)
                            gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTg                                                                  <  2:1271877/100‑1 (MQ=255)
                             aTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGTTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTgg                                                                 <  1:67264/100‑1 (MQ=255)
                               aaCATCAACGTCTCGGCGCTGGCGCCCGTCAGTCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGcc                                                               <  1:903971/100‑1 (MQ=255)
                                 cATCAACGACTCGGCGCTGGCGTCCGTCAGCTGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                             <  2:991462/100‑1 (MQ=255)
                                   ncAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGtt                                                           >  1:281284/2‑100 (MQ=255)
                                    cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATTACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGTTTGGCCTGttt                                                          >  1:1152042/1‑100 (MQ=255)
                                          cTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                                                     <  1:1061223/99‑1 (MQ=255)
                                            cGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGc                                                  >  2:1550470/1‑100 (MQ=255)
                                            cGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGc                                                  >  2:1368726/1‑100 (MQ=255)
                                             ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                 >  2:1136664/1‑100 (MQ=255)
                                                  tGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGc                                            <  2:291692/100‑1 (MQ=255)
                                                  tGGCGTCCGCCAGCCGGGTGATCACTCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGc                                            <  2:1141863/100‑1 (MQ=255)
                                                   ggCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                           >  1:263049/1‑100 (MQ=255)
                                                     cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                         >  2:1205757/1‑100 (MQ=255)
                                                     cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGAATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                         <  2:130022/100‑1 (MQ=255)
                                                       tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaag                                       <  1:846624/100‑1 (MQ=255)
                                                        ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaaga                                      >  1:1229043/1‑100 (MQ=255)
                                                         cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                     >  2:1417750/1‑100 (MQ=255)
                                                              gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                <  1:1683655/100‑1 (MQ=255)
                                                               ccGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTa                               >  2:135654/1‑100 (MQ=255)
                                                               ccGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTa                               >  2:761682/1‑100 (MQ=255)
                                                                cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAg                              <  2:1351284/100‑1 (MQ=255)
                                                                cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAg                              >  1:1798182/1‑100 (MQ=255)
                                                                  ggttgtCACCCCGGTCTTCAAGCCAATTGGCAGGTATGAAGATAAGTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCACGGCGAAAGAAGTGGTAGTg                            <  2:877517/95‑1 (MQ=255)
                                                                  ggTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                            <  1:529993/100‑1 (MQ=255)
                                                                  ggTGATCACCCCGGTCTTCAAGCCAACTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                            <  1:759976/100‑1 (MQ=255)
                                                                   gTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                            <  1:325975/99‑1 (MQ=255)
                                                                    tGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                          <  2:300988/100‑1 (MQ=255)
                                                                     gATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGt                         >  1:539663/1‑100 (MQ=255)
                                                                      aTCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                        <  2:1817740/100‑1 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTaag                      >  1:1810178/1‑96 (MQ=255)
                                                                          ccccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGc                     >  1:742991/1‑99 (MQ=255)
                                                                           cccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                   >  1:1173789/1‑100 (MQ=255)
                                                                           cccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                   >  1:259019/1‑100 (MQ=255)
                                                                             cGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCaa                 <  1:1136653/100‑1 (MQ=255)
                                                                                tCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac              <  2:1609764/100‑1 (MQ=255)
                                                                                  ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc            <  2:263050/100‑1 (MQ=255)
                                                                                  ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc            <  1:462523/100‑1 (MQ=255)
                                                                                   tCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct           <  1:1661779/100‑1 (MQ=255)
                                                                                     aaGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct         <  1:1401897/100‑1 (MQ=255)
                                                                                      aGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa        >  2:1215751/1‑100 (MQ=255)
                                                                                         cAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac     >  1:39231/1‑100 (MQ=255)
                                                                                         cAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac     >  1:1734360/1‑100 (MQ=255)
                                                                                           aTTGGCATGTATGAAGATGACTGGCAGGTAACGGTTGGACTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCg   <  2:1204527/100‑1 (MQ=255)
                                                                                            ttGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc  <  1:90519/100‑1 (MQ=255)
                                                                                                 |                                                                                              
CAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGC  >  NC_000913/3541996‑3542187

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