Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,144 Δ1 bp intergenic (+146/+305) hemN → / ← glnG O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,1430C.98.8% 257.0 / ‑3.1 84intergenic (+145/+306)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (43/40);  minor base T (0/1);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC  >  NC_000913/4051069‑4051212
                                                                          |                                                                      
tGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                       <  2:3877512/76‑2 (MQ=255)
  aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                     >  1:1715395/1‑73 (MQ=255)
  aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                     <  1:2619220/76‑4 (MQ=255)
    gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                  <  1:2654158/76‑1 (MQ=255)
     cgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGAGTTGACGCCGCATCAGGCATCTGa                                                                 >  1:582541/1‑76 (MQ=25)
      gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAg                                                                >  1:3405857/1‑76 (MQ=255)
      gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAg                                                                <  2:2271789/76‑1 (MQ=255)
      gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAg                                                                <  1:965992/76‑1 (MQ=255)
      gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAg                                                                <  1:3285434/76‑1 (MQ=255)
       cTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGc                                                               >  1:2699710/1‑76 (MQ=255)
         gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCat                                                             <  1:1437702/76‑1 (MQ=255)
           gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                           <  2:723146/76‑1 (MQ=255)
           gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                           <  2:3989893/76‑1 (MQ=255)
             aacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCCTCAGGCATCTGACCatgatc                                                         >  2:168289/1‑75 (MQ=21)
              acaaAACCTGTTGAACTTATTGCTTTTGTAGGCCGGATAATGCGTTTACGCCGCATCAGGCATCTGAGCATGATGc                                                        <  2:3636233/76‑1 (MQ=25)
              acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGc                                                        >  2:3997036/1‑76 (MQ=255)
                 aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTg                                                     <  2:4831980/76‑1 (MQ=255)
                  aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                    >  1:4052777/1‑76 (MQ=255)
                  aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                    <  1:1340798/76‑1 (MQ=255)
                   agcTGTTGAGATTATTGCTTTTGTATGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGTGCATGATGCCTGAt                                                   <  2:4242455/74‑1 (MQ=21)
                   aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAt                                                   >  1:3576076/1‑76 (MQ=255)
                    ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATg                                                  <  2:2677339/76‑1 (MQ=255)
                     cTGTTGCAATTATTGCTTTTGTAGTGCGGGTAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgc                                                 <  2:3416694/76‑1 (MQ=21)
                     cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgc                                                 <  1:1892130/76‑1 (MQ=255)
                      tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCTGACCCTGCTGCATGATgcg                                                >  2:681713/1‑76 (MQ=17)
                      tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcg                                                >  2:2371839/1‑76 (MQ=255)
                      tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcg                                                >  2:654059/1‑76 (MQ=255)
                      tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcg                                                <  2:395081/76‑1 (MQ=255)
                       gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcgg                                               <  1:4415664/76‑1 (MQ=255)
                       gTTGAAATTATTGCTTGTGTAGGCCGGATAAGGCGGTTACGCCGCATCAGGCATCTGAGCATGAGGCCTGATgcgg                                               <  1:4429157/76‑1 (MQ=25)
                        ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                              <  2:4296059/76‑1 (MQ=255)
                          gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgc                                            <  1:2178364/76‑2 (MQ=255)
                          gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgc                                            <  2:1655209/76‑2 (MQ=255)
                            aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctg                                          <  2:2082580/76‑4 (MQ=21)
                                attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                      <  2:1749076/76‑1 (MQ=18)
                                attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                      <  2:2723264/76‑1 (MQ=21)
                                  tggttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                    <  2:1253639/73‑1 (MQ=18)
                                   gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                   >  1:2275095/1‑76 (MQ=21)
                                    cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                  <  1:1701747/76‑1 (MQ=21)
                                     ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGGCt                                 >  1:2443486/1‑76 (MQ=18)
                                     ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCCGCGCTCGCCGCGTCt                                 >  2:410167/1‑76 (MQ=17)
                                      tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                >  1:3874497/1‑76 (MQ=18)
                                      tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                >  1:3966757/1‑76 (MQ=18)
                                       ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                               >  1:3672564/1‑76 (MQ=18)
                                        tGTAGGTCGGCTAATGCGTTTACTCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                              <  1:1146967/76‑1 (MQ=17)
                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                              <  1:1169290/76‑1 (MQ=18)
                                          tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                            <  2:2961998/76‑1 (MQ=18)
                                           aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                           >  1:3455993/1‑76 (MQ=18)
                                           aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                           >  1:1998249/1‑76 (MQ=18)
                                           aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                           >  2:2907271/1‑76 (MQ=17)
                                            ggCGGGATGGGGCGTTGACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                          <  2:257473/76‑1 (MQ=14)
                                            ggCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                          <  1:4817138/76‑1 (MQ=17)
                                            ggCCGAATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGTGGCGCTGGTCGCGTCTTATCAgg                          <  1:4821296/76‑1 (MQ=17)
                                             gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                         <  2:3917072/76‑1 (MQ=17)
                                               cGGATAAGGCGTTTACGCCGCATCAGTCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                       <  2:1888872/76‑1 (MQ=17)
                                                 gATAGGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                     <  2:426568/76‑1 (MQ=14)
                                                    aaGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                  <  2:3029404/76‑1 (MQ=17)
                                                     aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                 >  2:2724954/1‑76 (MQ=18)
                                                         gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct             >  1:2860202/1‑76 (MQ=18)
                                                         gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct             >  2:2968455/1‑76 (MQ=18)
                                                         gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct             >  2:3002451/1‑76 (MQ=18)
                                                         gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGCCTTCTCCGCGCCACGACtttgc             >  2:1746167/1‑74 (MQ=9)
                                                          tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTtctc            >  2:444669/1‑76 (MQ=17)
                                                          tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc            <  2:2496165/76‑1 (MQ=18)
                                                          tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc            <  2:3650588/76‑1 (MQ=18)
                                                          tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCCACGACTTtctc            >  1:4753504/1‑76 (MQ=18)
                                                          tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGAGCTGGTCGCGTCTTATCAGGCCTACGACCTtctc            >  2:566171/1‑76 (MQ=18)
                                                          tttACGCCGCATCAGGCAGCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc            <  2:2760819/76‑1 (MQ=18)
                                                           ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCCGGGCTACGACTTTCTcc           >  2:878129/1‑76 (MQ=18)
                                                           ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc           >  1:1396638/1‑76 (MQ=21)
                                                            tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg          >  2:1792341/1‑76 (MQ=21)
                                                            tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg          >  2:1757836/1‑76 (MQ=21)
                                                             aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc         >  2:2821285/1‑76 (MQ=21)
                                                             aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc         <  2:1594107/76‑1 (MQ=18)
                                                              cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa        >  2:3980126/1‑76 (MQ=21)
                                                              cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa        <  1:673766/76‑1 (MQ=21)
                                                               gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt       >  1:64683/1‑76 (MQ=21)
                                                               gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt       >  1:4782650/1‑76 (MQ=21)
                                                               gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt       >  2:76887/1‑76 (MQ=21)
                                                                 cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc     <  1:3355365/76‑1 (MQ=21)
                                                                  gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg    >  1:647711/1‑76 (MQ=25)
                                                                   cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCCTcccg   >  2:284885/1‑74 (MQ=21)
                                                                   cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg   >  2:277183/1‑76 (MQ=25)
                                                                    aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc  >  2:3559681/1‑76 (MQ=25)
                                                                    aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc  >  2:3289364/1‑76 (MQ=25)
                                                                    aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc  >  2:1436396/1‑76 (MQ=25)
                                                                    aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc  >  1:1465514/1‑76 (MQ=25)
                                                                          |                                                                      
TGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC  >  NC_000913/4051069‑4051212

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