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)
     aTGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATg                                                                                         >  1:1103982/1‑100 (MQ=255)
      tGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATgg                                                                                        >  2:942361/1‑100 (MQ=255)
       naaaCCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATggg                                                                                       >  2:266438/2‑100 (MQ=255)
       gAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATggg                                                                                       >  2:405099/1‑100 (MQ=255)
       gAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATggg                                                                                       >  2:720675/1‑100 (MQ=255)
         aaCCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCg                                                                                     >  1:533762/1‑100 (MQ=255)
            cTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGa                                                                                  <  2:86723/100‑1 (MQ=255)
               gCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATgg                                                                               <  1:726186/100‑1 (MQ=255)
                   agCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCa                                                                           <  2:1046233/100‑1 (MQ=255)
                   agCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCa                                                                           >  1:790426/1‑100 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAg                                                                          >  1:1234875/1‑100 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAg                                                                          >  2:569035/1‑100 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAg                                                                          <  1:57897/100‑1 (MQ=255)
                    gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAg                                                                          <  1:57928/100‑1 (MQ=255)
                    gCCTGAAAGACACCTGCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATc                                                                            <  2:1243031/98‑1 (MQ=255)
                     ccTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGACCAGa                                                                         >  2:901754/1‑100 (MQ=255)
                      cTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGaa                                                                        <  1:569039/100‑1 (MQ=255)
                        gAAAGACACCTTCAGCCTCAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGaaa                                                                       <  2:1123868/99‑1 (MQ=255)
                            gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCg                                                                   >  2:569676/1‑100 (MQ=255)
                            gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCgg                                                                  <  2:32537/100‑1 (MQ=255)
                                 cTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAgcg                                                             >  2:58091/1‑100 (MQ=255)
                                   tCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgca                                                           <  1:683972/100‑1 (MQ=255)
                                    cAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcag                                                          >  2:270412/1‑100 (MQ=255)
                                     aGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                         <  1:1152373/100‑1 (MQ=255)
                                     aGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                         <  2:783758/100‑1 (MQ=255)
                                      gCCTTAGCGTACTGATGAACCCCAGTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                         <  2:602736/99‑1 (MQ=255)
                                       ccTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagcag                                                       >  1:1059849/1‑100 (MQ=255)
                                       ccTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagcag                                                       <  1:491304/100‑1 (MQ=255)
                                                 cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                             >  1:431549/1‑100 (MQ=255)
                                                 cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                             <  1:1131956/100‑1 (MQ=255)
                                                    atgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTa                                          <  2:551447/100‑1 (MQ=255)
                                                     tgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTAc                                         <  2:1304288/100‑1 (MQ=255)
                                                        aCCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTg                                      >  1:704618/1‑100 (MQ=255)
                                                           ctaTTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAtt                                   <  2:793251/98‑1 (MQ=255)
                                                             aTTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAttt                                  <  1:270415/99‑1 (MQ=255)
                                                              ttGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTc                                <  2:533758/100‑1 (MQ=255)
                                                                gaggCCAGCAGAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCg                               <  2:819246/96‑1 (MQ=255)
                                                                gAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                              >  1:264791/1‑100 (MQ=255)
                                                                  aGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                            <  2:435620/100‑1 (MQ=255)
                                                                   gCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCct                           >  2:974847/1‑100 (MQ=255)
                                                                      aGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                        >  1:620307/1‑100 (MQ=255)
                                                                         aaaGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCtgta                     <  2:845717/100‑1 (MQ=255)
                                                                              cGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAc                >  1:153104/1‑100 (MQ=255)
                                                                               gACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAcc               >  1:660886/1‑100 (MQ=255)
                                                                                aCGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCa              <  1:106824/100‑1 (MQ=255)
                                                                                 cGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCAt             >  1:1188539/1‑100 (MQ=255)
                                                                                  ggCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATg            <  2:285082/100‑1 (MQ=255)
                                                                                    cGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa          >  1:255730/1‑100 (MQ=255)
                                                                                      aaaTGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTAtc        >  2:952565/1‑100 (MQ=255)
                                                                                          gggTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg    >  2:631207/1‑100 (MQ=255)
                                                                                          gggTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg    <  1:780361/100‑1 (MQ=255)
                                                                                           ggTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg    >  2:108657/1‑100 (MQ=255)
                                                                                            gTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTg  >  2:1061251/1‑100 (MQ=255)
                                                                                            gTACC‑GGGCGATGGGCGTGATGGATCAGAAACTCGGTAGCGCAGAAGTAGCTTACAGCTACCTGATTTTCCTCCTACTGTATGTACCATGTATCTCGGTg  <  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: ‑