Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,170 +C intergenic (+172/+279) 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,1701.C99.1% 277.0 / ‑3.5 114intergenic (+172/+279)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (63/50);  total (64/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 128 total aligned reads displayed.

TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC  >  NC_000913/4051094‑4051239
                                                                             |                                                                      
tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                                         <  1:344628/76‑1 (MQ=255)
tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                                         >  2:2102151/1‑76 (MQ=255)
 gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgc                                                                        <  2:2440166/76‑2 (MQ=255)
 gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgc                                                                        <  1:957522/76‑2 (MQ=255)
  gaaTTGTTGTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgct                                                                       <  2:5684933/75‑3 (MQ=18)
  aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgct                                                                       <  1:1738880/76‑3 (MQ=255)
  aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgct                                                                       >  2:753531/1‑74 (MQ=255)
         tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                               >  1:3452700/1‑76 (MQ=18)
         tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                               <  1:460886/76‑1 (MQ=18)
          ggttttGTGGGGCGGATAAGGCGTTTGCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                              <  2:5323602/74‑1 (MQ=17)
             tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                           <  2:3611446/76‑1 (MQ=18)
              ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtca                                                          >  1:3056857/1‑74 (MQ=18)
              ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                          <  2:4424094/76‑1 (MQ=18)
               tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                         >  1:3637006/1‑76 (MQ=18)
                gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                        >  2:4395508/1‑76 (MQ=18)
                 tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                       <  2:2256642/76‑1 (MQ=17)
                  aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                      >  2:5162777/1‑76 (MQ=17)
                  aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                      >  1:1153518/1‑76 (MQ=17)
                    gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCCGGTCGCGTCCTCTCCGGc                                                    >  2:680959/1‑76 (MQ=17)
                     ccGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                   >  2:5678996/1‑76 (MQ=17)
                      cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                  <  2:5187965/76‑1 (MQ=17)
                       ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                 >  1:790874/1‑76 (MQ=17)
                       ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                 <  1:2170883/76‑1 (MQ=17)
                       ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGACCCGGCGCCGGGCGCGGCTTCTCCGGccca                                                 >  2:17670/1‑74 (MQ=2)
                        gATGAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                <  2:1057782/76‑1 (MQ=14)
                        gATAGGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                <  2:5229387/76‑1 (MQ=14)
                        gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                >  1:12309/1‑76 (MQ=17)
                        gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                >  1:1976523/1‑76 (MQ=17)
                          tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGCTGCCTGATGCGGCGCTGGGCGCGTCTTATCAGGCCCACGa                                              >  1:1242545/1‑76 (MQ=17)
                          tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                              <  1:2467049/76‑1 (MQ=17)
                           aaGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                             <  2:4892949/76‑1 (MQ=17)
                            aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                            >  1:381072/1‑76 (MQ=17)
                              gCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                          >  2:5241175/1‑76 (MQ=18)
                              gCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                          <  1:5298661/76‑1 (MQ=18)
                                gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                >  1:1429169/1‑68 (MQ=17)
                                gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                <  2:1429169/68‑1 (MQ=17)
                                gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                        >  2:1574908/1‑76 (MQ=18)
                                  ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                      >  1:4845867/1‑76 (MQ=18)
                                  ttACGCAGCATAAGGCGTCTGAGCATGATGCCTGATGAGGCGCTGGTCGCGTATTATCAGGCCTAAGACTTTCTcc                                      <  2:1153518/76‑1 (MQ=16)
                                   tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                     <  1:1410198/76‑1 (MQ=18)
                                   tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                     <  2:4042498/76‑1 (MQ=18)
                                   tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                     >  1:4161388/1‑76 (MQ=21)
                                    aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGTGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                                    >  2:2705991/1‑76 (MQ=21)
                                    aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                                    >  1:3222830/1‑76 (MQ=21)
                                     cgccggATCAGGCATCTGAGCATGATGCCTGAGGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                   <  2:3002479/76‑1 (MQ=21)
                                     cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                   >  2:3212946/1‑76 (MQ=21)
                                      gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  >  2:915616/1‑76 (MQ=21)
                                      gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  <  1:4542714/76‑1 (MQ=21)
                                      gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  >  1:3170321/1‑76 (MQ=21)
                                      cccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                  <  1:1223799/75‑1 (MQ=21)
                                       ccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                 >  2:4045179/1‑76 (MQ=25)
                                        cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                >  1:5790573/1‑76 (MQ=25)
                                         gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                               >  2:1887012/1‑76 (MQ=25)
                                         gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                               <  1:2894823/76‑1 (MQ=25)
                                          cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                              <  1:5678998/76‑1 (MQ=25)
                                          cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                              >  2:110897/1‑76 (MQ=25)
                                           aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             >  2:5038566/1‑76 (MQ=25)
                                           aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGATCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             <  2:4335093/76‑1 (MQ=18)
                                           aTCAGGCATCAGAGCATGATGCCTGATGCGGCGCTGGGCGCGTATTATCAGGCCTACGACTTTCTCCGACTCCGGc                             >  2:5045465/1‑76 (MQ=18)
                                           aTCAGGCAGCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             <  2:1867130/76‑1 (MQ=25)
                                            tCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCa                            <  1:2421496/76‑1 (MQ=25)
                                               gggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                          <  2:183800/73‑1 (MQ=25)
                                               aggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                          <  2:381072/73‑1 (MQ=25)
                                               aggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTGTCTCCGCATCCGGCATa                          >  1:4207133/4‑76 (MQ=25)
                                                ggcATCTGAGCTTGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                         >  1:7522/3‑76 (MQ=25)
                                                ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                         <  2:4643181/74‑1 (MQ=37)
                                                ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                         <  2:2430103/74‑1 (MQ=25)
                                                ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                         <  1:3802799/74‑1 (MQ=25)
                                                 gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa                        >  2:2693411/2‑76 (MQ=255)
                                                  cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCCACTCCGGCATAaac                       >  2:1119013/1‑76 (MQ=21)
                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                     >  1:980994/1‑76 (MQ=255)
                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                     >  2:2912078/1‑76 (MQ=255)
                                                    tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                     >  1:4798750/1‑76 (MQ=255)
                                                     cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                    <  2:3033947/76‑1 (MQ=255)
                                                     cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                    <  1:3255270/76‑1 (MQ=255)
                                                     cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                    >  2:2407061/1‑76 (MQ=255)
                                                     cTGAGCATCACGCCTGATGCGGCGCTGGGCGCGACTTATCAGGCCTACGACAATCTCCGCATCCGGCATAaaccac                    >  2:3173901/1‑76 (MQ=16)
                                                      tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                   >  1:3399480/1‑76 (MQ=255)
                                                      tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                   <  1:2005735/76‑1 (MQ=255)
                                                        aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg                 <  1:753531/76‑1 (MQ=255)
                                                         gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGTATCCGGCATAAACAACGAGc                <  1:5235989/76‑1 (MQ=25)
                                                           atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt              >  1:4197494/1‑76 (MQ=255)
                                                           atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt              >  1:2812979/1‑76 (MQ=255)
                                                            tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATACACAACGAGCTTc             >  2:3095213/1‑76 (MQ=37)
                                                             gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa            >  2:4528259/1‑76 (MQ=255)
                                                             gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa            >  2:267844/1‑76 (MQ=255)
                                                                gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt         <  2:1493833/76‑1 (MQ=255)
                                                                gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt         >  2:1813033/1‑76 (MQ=255)
                                                                 gcTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa        <  2:3309992/75‑1 (MQ=255)
                                                                  cggATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaa       <  2:3231825/74‑1 (MQ=38)
                                                                  cTGAGGCGGCGAGGGCGACGCCTTATCACGCCTGCGGCTCTCGCCGCAGCCGGCATAAACAACGAGCGTCAGCCaa       >  2:318237/1‑76 (MQ=9)
                                                                   tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa      >  2:3431960/1‑76 (MQ=255)
                                                                    gATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc     >  2:4664341/1‑76 (MQ=255)
                                                                    gATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc     <  1:414984/76‑1 (MQ=255)
                                                                     aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             <  2:1747733/51‑1 (MQ=25)
                                                                     aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                             >  1:1747733/1‑51 (MQ=37)
                                                                     aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa    >  2:761517/1‑76 (MQ=255)
                                                                     aTGCGGCGCGGGTCGCGTCTTATCGGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa    <  2:1247476/76‑1 (MQ=38)
                                                                       gcggcgcTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc  <  2:3452700/76‑1 (MQ=255)
                                                                       gcggcgcTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCACAAACAACGAGCTTCAGCTAACCAGc  >  2:5494614/1‑76 (MQ=38)
                                                                             |                                                                      
TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC  >  NC_000913/4051094‑4051239

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