Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,537,405 +GC intergenic (+61/‑394) b3407 → / → feoA hypothetical protein/ferrous iron transport protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,537,4041.C98.6% 170.6 / ‑3.8 69intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (21/47);  total (21/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,537,4042.G98.6% 174.8 / ‑3.5 69intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (21/47);  total (21/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTA  >  NC_000913/3537329‑3537468
                                                                              ||                                                                
gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑Ca                                                                    >  1:4002281/1‑76 (MQ=255)
gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑Ca                                                                    >  2:2878587/1‑76 (MQ=255)
 tgCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑c                                                                 >  2:3041947/2‑76 (MQ=255)
 ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑c                                                                 <  2:713511/76‑1 (MQ=255)
   cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cgc                                                               <  2:4377369/76‑3 (MQ=255)
    aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cgcc                                                              <  1:4849180/76‑4 (MQ=255)
    aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑cgcc                                                              <  2:2546973/76‑4 (MQ=255)
    aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGCAGGGCAGATAAAGCGTT‑‑CA‑‑cgcc                                                              >  1:2673450/1‑73 (MQ=255)
           gTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTg                                                         <  2:2917401/76‑1 (MQ=37)
             tAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGc                                                       <  1:2820820/76‑1 (MQ=37)
             tAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGc                                                       <  1:2976323/76‑1 (MQ=37)
             tAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGc                                                       <  2:2979591/76‑1 (MQ=37)
               aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaa                                                     <  1:4364844/76‑1 (MQ=25)
                acccccGTACGTGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaac                                                    <  1:303693/73‑1 (MQ=21)
                aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaac                                                    <  2:4188871/76‑1 (MQ=25)
                 cGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacg                                                   >  1:1018778/1‑76 (MQ=25)
                  gagcGTTCTGGGGTACGGTTGGAATTTGCTCGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacga                                                  <  2:1866286/73‑1 (MQ=12)
                  gCCCGTACTGGCCTACGGTTGGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCaacgg                                                  <  1:2887173/76‑2 (MQ=25)
                      gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGc                                              <  1:4382326/76‑1 (MQ=25)
                      gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGc                                              <  2:3680454/76‑1 (MQ=25)
                        aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCt                                            >  2:1700362/1‑76 (MQ=25)
                          tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCttc                                          >  2:4004602/1‑75 (MQ=25)
                           ggCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGc                                         <  2:2215673/76‑1 (MQ=25)
                             ccgaCGGGTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCt                                       <  2:296946/73‑1 (MQ=21)
                             ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCt                                       <  2:2544883/76‑1 (MQ=25)
                              cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTg                                      <  2:1731013/76‑1 (MQ=25)
                              cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTg                                      <  2:2867886/76‑1 (MQ=25)
                              cTACGGTTCGAATTTGCACGAAATCGGAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTg                                      <  2:4395866/76‑1 (MQ=25)
                               tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGa                                     <  1:2742892/76‑1 (MQ=25)
                               tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGa                                     >  1:175180/1‑76 (MQ=25)
                                  ggTTGGAATTTGAACGAAATCGTAGGGAAGATAAGGCGTT‑‑CACGCCGCATCTGGNNACGAACGCCTTGCCTGACAt                                  <  1:265855/76‑1 (MQ=21)
                                  ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACAt                                  >  2:2782120/1‑76 (MQ=25)
                                  ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACAt                                  <  1:1870658/76‑1 (MQ=25)
                                  ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACAt                                  >  2:3571238/1‑76 (MQ=25)
                                  ggTTCGAATTTGCACGAAAGCGAAGGGCAGACAAAGCGTT‑‑CACGCCGCATCTCGCAACGAACGCCTTGCCTGACaa                                  >  1:2148667/1‑75 (MQ=255)
                                   gTTCGAATTTGCACTAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATa                                 <  1:851999/76‑1 (MQ=25)
                                    ttCGACTTTGCACGAAATCGTAGGGCCGATAAGCCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaa                                <  1:38666/76‑1 (MQ=25)
                                      cGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                              <  1:4356463/76‑1 (MQ=25)
                                       gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                             <  2:3179850/76‑1 (MQ=25)
                                       gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                             <  2:3914725/76‑1 (MQ=25)
                                        aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                            <  2:582900/76‑1 (MQ=25)
                                        aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                            <  2:2862630/76‑1 (MQ=25)
                                        aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTCCCTGACATAAAAgg                            >  2:299648/1‑75 (MQ=25)
                                         aTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                           >  1:2382783/1‑76 (MQ=25)
                                         aTTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                           <  2:2753211/76‑1 (MQ=25)
                                          tttGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                          <  1:2227012/76‑1 (MQ=25)
                                           ttGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTccc                         >  2:28955/1‑76 (MQ=25)
                                           ttGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                         <  1:3355120/76‑1 (MQ=25)
                                           ttGCACGAAATCGGAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                         <  1:2095496/76‑1 (MQ=25)
                                            tGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCg                        <  1:1109822/76‑1 (MQ=25)
                                               aCGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgag                     <  1:432183/76‑1 (MQ=25)
                                               aCGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgag                     <  1:1355778/76‑1 (MQ=25)
                                                cGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                    <  1:1101043/76‑1 (MQ=25)
                                                cGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                    <  2:3277835/76‑1 (MQ=25)
                                                cGAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGCGCCGgaga                    >  2:4846197/1‑76 (MQ=25)
                                                 gAAATCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                   <  2:3169628/76‑1 (MQ=25)
                                                 gAAATCGTAGGGCAGAGAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                   <  2:922704/76‑1 (MQ=25)
                                                 gAAATCGGAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                   <  2:1060853/76‑1 (MQ=25)
                                                  aaaTCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                  <  2:2171669/76‑1 (MQ=25)
                                                  aaaTCGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAat                  <  1:3322797/76‑1 (MQ=25)
                                                   aaTCGGAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                 <  2:4156898/76‑1 (MQ=25)
                                                      cGTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCNGGAGAATAtct              >  2:4150730/1‑76 (MQ=25)
                                                       gTAGGGCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc             >  1:2758650/1‑76 (MQ=37)
                                                         gggggAGATAAGGGGTTCGCG‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg           <  2:527055/75‑1 (MQ=37)
                                                          gggCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCgg          <  1:2501034/76‑1 (MQ=37)
                                                           ggCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGc         >  1:4543958/1‑76 (MQ=37)
                                                            gCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCa        <  1:4563627/76‑1 (MQ=37)
                                                            gCAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCa        >  1:3068975/1‑76 (MQ=37)
                                                             cAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAt       >  1:4082192/1‑76 (MQ=37)
                                                             cAGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAt       >  1:4076746/1‑76 (MQ=37)
                                                              aGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt      >  1:4636547/1‑76 (MQ=37)
                                                              aGATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt      >  2:1246484/1‑76 (MQ=37)
                                                               gATAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttt     <  2:1351950/76‑1 (MQ=38)
                                                                 tAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttttt   >  1:4358422/1‑76 (MQ=38)
                                                                 tAAGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttttt   <  1:701616/76‑1 (MQ=38)
                                                                  aaGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa  >  1:1222007/1‑76 (MQ=38)
                                                                  aaGGCGTT‑‑CACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa  <  2:1726284/76‑1 (MQ=38)
                                                                              ||                                                                
GGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTT‑‑CA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTA  >  NC_000913/3537329‑3537468

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