Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,173 A→G 12.0% intergenic (+57/+27) chiQ → / ← fur lipoprotein ChiQ/DNA‑binding transcriptional dual regulator Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,1730AG12.0% 106.4 / 6.5 52intergenic (+57/+27)chiQ/furlipoprotein ChiQ/DNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base G (3/3);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

CAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCA  >  NC_000913/710037‑710305
                                                                                                                                         |                                                                                                                                    
cAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTgg                                                                                                                                >  1:477657/1‑143 (MQ=255)
cAACGCAAACCGGAAATAATCAGGCGGTGAAAGCCGACTGTGATACAGTCTGGCAGGAAATACGCAGTAATCACAAATAAATGAGAGCTGTAAGTCTCGCTTTGCTTATTCACCACGCACAAACACGACAAGCGGCAGGTCgg                                                                                                                                >  2:891528/1‑143 (MQ=255)
  aCGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCt                                                                                                                              <  1:113524/143‑1 (MQ=255)
  aCGCAAACCGGAAATGAGCGGGCGGTGGAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCt                                                                                                                              <  2:297776/143‑1 (MQ=255)
  aCGCAAACCGGAAATGACCAGGCCGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATCTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCt                                                                                                                              <  2:107709/143‑1 (MQ=255)
   cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCGCGCTTTGCTTATGTAACAAACAAATATTTGTCAACCCGCAGGTTGGCtt                                                                                                                             >  1:840278/1‑143 (MQ=255)
    gCAAACCGGAAATGATCAGGCGGTCAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCGTGCATAAAAAAGCCAACCTGCGGGTTGGCttt                                                                                                                            <  2:496154/143‑1 (MQ=255)
     cAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGCCtttt                                                                                                                           >  2:750682/1‑143 (MQ=255)
     cAAACCGGAAATGAGCAGTCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAAGAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCtttt                                                                                                                           <  1:176673/143‑1 (MQ=255)
      aaaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                                                          <  1:812316/143‑1 (MQ=255)
       aaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtct                                                                                                                         >  1:220415/1‑143 (MQ=255)
       aaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtct                                                                                                                         >  2:266569/1‑143 (MQ=255)
               gtGATCAGGCGGTCACAGCCGATTGTGATCAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAgg                                                                                                                 <  2:196180/142‑1 (MQ=255)
                       gCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCCGGTTGGCTTTTCTCGTTCAGGCTGGCTTa                                                                                                         >  1:650441/1‑143 (MQ=255)
                       gCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTa                                                                                                         >  1:200853/1‑143 (MQ=255)
                       gCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTGTCTTATTTCCCCTGCTTAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTa                                                                                                         >  1:735688/1‑143 (MQ=255)
                         ggTGAAAGCCGATTGTGATAAAATATGCATTTATTGGCGCTCGAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAGAACAGCCAACCCGCAGGTTCGCTTTTCTCGTTCAGGCTGGCTTact                                                                                                       >  2:925154/1‑141 (MQ=255)
                               aGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATCACAAATAAGTGAGAGGTGTAACTCTCGCTGTTCTTATTTCTCTTGCATACAAAAGACAACCCGCTGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTAGCCtt                                                                                                 >  2:347126/1‑143 (MQ=255)
                                gCCTATTGTGATTAAGTCTGGCAGGAAATACGCAGTAATAACAAATATGTGATAGCTGTAACTGTCGCGTTTCTTATTTCCCTTGCATAAAAAAGCCATCCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTc                                                                                                <  2:783589/143‑1 (MQ=255)
                                gCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTc                                                                                                >  2:427615/1‑143 (MQ=255)
                                       gtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCa                                                                                         <  1:904119/143‑1 (MQ=255)
                                       gtgATAAAGTCGGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCa                                                                                         <  2:309783/143‑1 (MQ=255)
                                         gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCGGGTTGGCTTTTCTCGTTCAGTCTGGCTTATTTGCCTTCGTGCGCATg                                                                                       >  1:704670/1‑143 (MQ=255)
                                          aTAAAGTCTGTCGGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAGTTGCCTTCGTGCGCATGt                                                                                      >  2:719687/1‑143 (MQ=255)
                                              aGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                                                                  <  2:118080/143‑1 (MQ=255)
                                                   ggCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcg                                                                             <  1:158521/143‑1 (MQ=255)
                                                   ggCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcg                                                                             <  2:719150/143‑1 (MQ=255)
                                                    gCACGAAATACGCAGCAATAACAAGTAAGGGAGAGCTGTAACTCTCGCTGTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCGTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgc                                                                            <  2:300095/143‑1 (MQ=255)
                                                     cAGGAAATACGCAGTAATAACAAATGAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGGTGGCTTTTCGCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgcg                                                                           <  1:585908/143‑1 (MQ=255)
                                                     cAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTGATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgcg                                                                           <  2:399334/143‑1 (MQ=255)
                                                          aaTACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCCGGCTGGCTTATTTGCCTTCGTGCGCATTTTCATCTTCGCTGCAAt                                                                      >  1:240598/1‑143 (MQ=255)
                                                            tACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGACTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCg                                                                    >  1:914666/1‑143 (MQ=255)
                                                                        aaCAAATGAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGcaca                                                         <  1:648039/143‑1 (MQ=255)
                                                                        aaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAttt                                                                                                      >  1:378213/1‑98 (MQ=255)
                                                                        aaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAttt                                                                                                      <  2:378213/98‑1 (MQ=255)
                                                                            actaggTGAGCCCTGTAACTCTCCCTATACTTCTTGCCCTTGCATAAAACAGCCAACCCCCAGGTTGGCGTTTCTCGTTCAGGCTCGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGa                                                     <  1:31048/138‑1 (MQ=255)
                                                                             aTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGTCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAc                                                    <  1:419530/143‑1 (MQ=255)
                                                                                   gagaGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCTCTGGTTGGCTTTTCTTGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTaaa                                              <  1:36643/143‑1 (MQ=255)
                                                                                        cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATag                                         <  1:216792/143‑1 (MQ=255)
                                                                                                 tcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTgg                                <  2:38891/143‑1 (MQ=255)
                                                                                                   gCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTACCGATAGAGTCTGTGCtt                              >  1:114241/1‑143 (MQ=255)
                                                                                                   gCTTTTCTTATTTCCCTTGCATAAAAAAGCAAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCCCAGTGACCGTACAGATCGAGACTGTGGtc                              >  2:512858/1‑142 (MQ=255)
                                                                                                       ttCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACTGTATAGATAGAGACTGTGGTTAGTc                          >  2:724764/1‑143 (MQ=255)
                                                                                                         cTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCCTCTTCGCGGCAATCGCCTTCGGCACAGTGGCCGTACAGATAGATACTGTGGTTAGTGAg                        >  2:461173/1‑143 (MQ=255)
                                                                                                          ttATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAgg                       >  1:508529/1‑143 (MQ=255)
                                                                                                            aTTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTCTTTGCCTTCGTGCGCATGTTCATCTTCGCGGAAACCGCCTTTGGCTCAGTGGTCGTAACCATACAGACTCTGCTTAGTCAGGCg                     >  1:586470/1‑143 (MQ=255)
                                                                                                                   ttGCATAACAGAGCCAACCCGCAGGTTGGCTTTGCTCGTTCAGGCTGGCTTATTTGCCTTCGGGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCa              <  2:28473/143‑1 (MQ=255)
                                                                                                                   ttGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCa              >  1:878948/1‑143 (MQ=255)
                                                                                                                        tAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTGGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGttt         >  1:13681/1‑143 (MQ=255)
                                                                                                                           aaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGc      >  1:877631/1‑143 (MQ=255)
                                                                                                                            aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCg     >  1:169855/1‑143 (MQ=255)
                                                                                                                               gCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTCGCCTTCGTGCGCATGTTCATCTTCGCGGCACTCGCCTTCGGCACAGTGACCGTAAAGATAGCGCCGGTGGTTAGTCAGGCGCATGCCACGCTTAGCGGCa  >  2:325684/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCA  >  NC_000913/710037‑710305

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