Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,173 A→G 7.7% 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,1730AG 7.7% 139.0 / 3.2 56intergenic (+57/+27)chiQ/furlipoprotein ChiQ/DNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base A (33/18);  new base G (3/2);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

TGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCA  >  NC_000913/710031‑710311
                                                                                                                                              |                                                                                                                                          
tGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTCTTTTCCTTGCCTAAAAAAGCCAACCCGCa                                                                                                                                            >  2:223298/1‑143 (MQ=255)
 gCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTCGCATAAAACACCCAACCCGCAg                                                                                                                                           >  2:749620/1‑143 (MQ=255)
 gCGTGCAAAGCACACCGGAAAGGCTCAGGCGGTGAAGGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGGGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAg                                                                                                                                           <  1:681818/143‑1 (MQ=255)
    tGCAACGCAAACCGAAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAAAAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGtt                                                                                                                                        <  2:492017/143‑1 (MQ=255)
     gCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAACAAAGCCAACCCGCAGGTTg                                                                                                                                       >  1:872270/1‑143 (MQ=255)
       aaCGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCAGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGCTTGGc                                                                                                                                     >  2:250332/1‑143 (MQ=255)
       aaCGCAAACCGGAAATGATCAGGCGATGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                     >  2:866640/1‑143 (MQ=255)
            gaaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                                                                <  1:776190/142‑1 (MQ=255)
            aaaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                                                                <  1:19805/143‑1 (MQ=255)
              acaggAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtctc                                                                                                                              <  1:546674/140‑1 (MQ=255)
               ccGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGGTGGCTCTTCTCg                                                                                                                             >  1:315474/1‑143 (MQ=255)
                ctgAAATGATCAGGCCGTCAAAGCCGATTGTGATAAAGTCTGGCAGTAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGt                                                                                                                            <  1:93957/141‑1 (MQ=255)
                cGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCCACCCGCAGGTTGGCTTTTCTCGt                                                                                                                            >  1:771997/1‑143 (MQ=255)
                 ggAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTCTCTCATTTCCCTCACATAAACACGCCAACCCGCAGGCTGGCGTTtctcagt                                                                                                                           >  2:793465/1‑140 (MQ=255)
                  gAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTc                                                                                                                          >  1:606602/1‑143 (MQ=255)
                       gATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTGCAggct                                                                                                                     >  2:778878/1‑143 (MQ=255)
                           aGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctggct                                                                                                                 <  2:606665/143‑1 (MQ=255)
                                gTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAttt                                                                                                            <  1:348063/143‑1 (MQ=255)
                                 tGAAAGCCGATTGTGATAAAGTCTGGCATGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACGCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTg                                                                                                           <  2:121364/143‑1 (MQ=255)
                                          aTTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTgc                                                                                                  >  1:618367/1‑143 (MQ=255)
                                                aTAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTCTTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGt                                                                                            >  1:137298/1‑143 (MQ=255)
                                                 tAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGtt                                                                                           >  2:527400/1‑143 (MQ=255)
                                                      tCTGACAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCACGCTTTTCTTATTTCCCTTGCATACAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTCCCTTCTTGCGCATGTTGATCt                                                                                      >  1:635433/1‑143 (MQ=255)
                                                            aGGAAATACGCAGTACTAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                                                                >  2:866500/1‑143 (MQ=255)
                                                                    cGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTGCATCTTCGCGGcgggcgcc                                                                        >  1:25513/1‑136 (MQ=255)
                                                                              aCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAACCCAACCCGCACGTTGACTTTTCTCGTTCCCGCTGGCTTCTTTCCCTTCGTCCTCATGCTCATCTTCTCTGCACTCGCCTTCgctactg                                                              >  2:3530/1‑137 (MQ=255)
                                                                                aaaTAAGTGAGAGCTGTAACTCTCGCTCATCTTCTTTCACTTGCATAAAAAATCCAACCCGCAGGTTGACTGTTCTCGTTCATCCTGCCTTATTTGCCTTCGTGCGCGTGTGCATCTTCTCCGCTCTCGCCTTACGCCCAGTg                                                            >  2:668645/1‑143 (MQ=255)
                                                                                       tGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCAGCGCAGTAACCGtca                                                     >  2:661098/1‑141 (MQ=255)
                                                                                         agagCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAAGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCGTTCTCCTTCGGCCCAGTGACCGTagag                                                   >  1:235733/1‑143 (MQ=255)
                                                                                           agCTGTACCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGAt                                                 <  2:361552/143‑1 (MQ=255)
                                                                                             cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTGGCGGCAATCGCCTGCGGCACAGTGACCGTacgtagag                                               >  1:494312/1‑136 (MQ=255)
                                                                                             cTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATag                                               <  1:866640/143‑1 (MQ=255)
                                                                                             cTGTAACGCTCGCTTTTCTTATTCCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATag                                               <  1:250332/143‑1 (MQ=255)
                                                                                              tGTACCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga                                              <  2:89142/143‑1 (MQ=255)
                                                                                                 aaCTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGAc                                           >  1:387299/1‑143 (MQ=255)
                                                                                                    tctcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCGGACTGGCTTATTGGCCTTCGTGCGCATGTTCATCTTCGCGGGACTCGCCTGCGGCACAGTGACCGTACAGATACAGACtgc                                        >  2:304183/1‑142 (MQ=255)
                                                                                                      tcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTgg                                      >  1:507200/1‑143 (MQ=255)
                                                                                                      tcGCTTGTCTTATTTCCCTTGCATACAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATCTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTgg                                      <  2:847283/143‑1 (MQ=255)
                                                                                                             tCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCa                               <  1:714080/143‑1 (MQ=255)
                                                                                                               ttATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAgg                             >  1:100517/1‑143 (MQ=255)
                                                                                                                   ttCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGaa                         <  2:229389/143‑1 (MQ=255)
                                                                                                                   ttCCCTGGCATAAAAAAGCCAACCCGCAGGGTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGaa                         <  1:700360/143‑1 (MQ=255)
                                                                                                                        ttGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCATTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCa                    >  2:9069/1‑143 (MQ=255)
                                                                                                                          gCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCagg                  >  1:668717/1‑141 (MQ=255)
                                                                                                                            aTAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTACAGATAGAGCCAGTGCTTAGTCAGGCGAATGCCAtcat                >  2:37555/1‑140 (MQ=255)
                                                                                                                             tAAGAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGttt               <  1:832304/143‑1 (MQ=255)
                                                                                                                             gagaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGttt               >  1:144621/4‑143 (MQ=255)
                                                                                                                              agaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTCTTTGCCTTCGTGCGCATGTTCATATTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGACAGAGACTGGGGTTAGTCAGGCGAATGCCAtcgtgt              >  2:592501/3‑138 (MQ=255)
                                                                                                                              aaaaaaGCCAACCGGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTGCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGCTAGAGACTGTGGTTAGTCACGCGAATGCCATGtctt              >  2:821690/1‑143 (MQ=255)
                                                                                                                              aaaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGACTGCCATGtttt              >  2:491627/1‑143 (MQ=255)
                                                                                                                               gaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATg                       <  2:170263/132‑1 (MQ=255)
                                                                                                                               gaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATg                       >  1:170263/2‑133 (MQ=255)
                                                                                                                                 aaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCg           >  1:99847/1‑143 (MQ=255)
                                                                                                                                  aaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTGCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGCTAGAGACTGTGGTTAGTAAGAAGAAGCGCATGCTTTGCgg          >  2:804264/1‑143 (MQ=255)
                                                                                                                                      cAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTACAGCTAGAGGCTGTGGTTAGTCAGGCGACTTCCTCGCTGTGCGGCAAt      >  1:43453/1‑143 (MQ=255)
                                                                                                                                          ccGCGGGTTGGCTTTTCGCGTGCAGGCGGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCGTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCa  <  2:488596/143‑1 (MQ=255)
                                                                                                                                          ccGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTGAAGATAGAGACTGTGGTTAGTCGGGCGACTGCCATGTTTTGCGGCAATTTCa  >  1:736710/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                          
TGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCA  >  NC_000913/710031‑710311

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