Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,542,354 Δ1 bp coding (1948/2322 nt) 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,3510G.96.3% 236.5 / 2.9 54coding (1945/2322 nt)feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (25/27);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GATAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTG  >  NC_000913/3542254‑3542446
                                                                                                 |                                                                                               
gATAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCggg                                                                                               <  1:992942/100‑1 (MQ=255)
 atcgATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCgggc                                                                                              >  1:628271/4‑99 (MQ=255)
 aTAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCgggc                                                                                              <  1:1186803/100‑2 (MQ=255)
     aTGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATg                                                                                         >  1:1103982/1‑100 (MQ=255)
      tGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATgg                                                                                        >  2:942361/1‑100 (MQ=255)
       naaaCCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATggg                                                                                       >  2:266438/2‑100 (MQ=255)
       gAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATggg                                                                                       >  2:405099/1‑100 (MQ=255)
       gAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATggg                                                                                       >  2:720675/1‑100 (MQ=255)
         aaCCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCg                                                                                     >  1:533762/1‑100 (MQ=255)
            cTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGa                                                                                  <  2:86723/100‑1 (MQ=255)
               gCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATgg                                                                               <  1:726186/100‑1 (MQ=255)
                   agCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCa                                                                           <  2:1046233/100‑1 (MQ=255)
                   agCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCa                                                                           >  1:790426/1‑100 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAg                                                                          >  1:1234875/1‑100 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAg                                                                          >  2:569035/1‑100 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAg                                                                          <  1:57897/100‑1 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAg                                                                          <  1:57928/100‑1 (MQ=255)
                    gCCTGAAAGACACCTGCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATc                                                                            <  2:1243031/98‑1 (MQ=255)
                     ccTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGACCAGa                                                                         >  2:901754/1‑100 (MQ=255)
                      cTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGaa                                                                        <  1:569039/100‑1 (MQ=255)
                        gAAAGACACCTTCAGCCTCAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGaaa                                                                       <  2:1123868/99‑1 (MQ=255)
                            gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCg                                                                   >  2:569676/1‑100 (MQ=255)
                            gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCgg                                                                  <  2:32537/100‑1 (MQ=255)
                                 cTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAgcg                                                             >  2:58091/1‑100 (MQ=255)
                                   tCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgca                                                           <  1:683972/100‑1 (MQ=255)
                                    cAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcag                                                          >  2:270412/1‑100 (MQ=255)
                                     aGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                         <  1:1152373/100‑1 (MQ=255)
                                     aGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                         <  2:783758/100‑1 (MQ=255)
                                      gCCTTAGCGTACTGATGAACCCCAGTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                         <  2:602736/99‑1 (MQ=255)
                                       ccTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagcag                                                       >  1:1059849/1‑100 (MQ=255)
                                       ccTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagcag                                                       <  1:491304/100‑1 (MQ=255)
                                                 cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                             >  1:431549/1‑100 (MQ=255)
                                                 cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                             <  1:1131956/100‑1 (MQ=255)
                                                    atgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTa                                          <  2:551447/100‑1 (MQ=255)
                                                     tgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTAc                                         <  2:1304288/100‑1 (MQ=255)
                                                        aCCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTg                                      >  1:704618/1‑100 (MQ=255)
                                                           ctaTTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAtt                                   <  2:793251/98‑1 (MQ=255)
                                                             aTTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAttt                                  <  1:270415/99‑1 (MQ=255)
                                                              ttGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTc                                <  2:533758/100‑1 (MQ=255)
                                                                gaggCCAGCAGAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCg                               <  2:819246/96‑1 (MQ=255)
                                                                gAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                              >  1:264791/1‑100 (MQ=255)
                                                                  aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                            <  2:435620/100‑1 (MQ=255)
                                                                   gCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCct                           >  2:974847/1‑100 (MQ=255)
                                                                      aGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                        >  1:620307/1‑100 (MQ=255)
                                                                         aaaGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCtgta                     <  2:845717/100‑1 (MQ=255)
                                                                              cGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAc                >  1:153104/1‑100 (MQ=255)
                                                                               gACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAcc               >  1:660886/1‑100 (MQ=255)
                                                                                aCGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCa              <  1:106824/100‑1 (MQ=255)
                                                                                 cGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCAt             >  1:1188539/1‑100 (MQ=255)
                                                                                  ggCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATg            <  2:285082/100‑1 (MQ=255)
                                                                                    cGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa          >  1:255730/1‑100 (MQ=255)
                                                                                      aaaTGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTAtc        >  2:952565/1‑100 (MQ=255)
                                                                                          gggTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg    >  2:631207/1‑100 (MQ=255)
                                                                                          gggTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg    <  1:780361/100‑1 (MQ=255)
                                                                                           ggTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg    >  2:108657/1‑100 (MQ=255)
                                                                                            gTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTg  >  2:1061251/1‑100 (MQ=255)
                                                                                            gTACCGGGCGATGGGCGTGATGGATCAGAAACTCGGTAGCGCAGAAGTAGCTTACAGCTACCTGATTTTCCTCCTACTGTATGTACCATGTATCTCGGTg  <  2:146922/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GATAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTG  >  NC_000913/3542254‑3542446

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