Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 710,436 G→A H71Y (CAC→TAC)  fur ← DNA‑binding transcriptional dual regulator Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,4360GA98.2% 168.6 / NA 54H71Y (CAC→TAC) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/23);  minor base C (1/0);  total (31/23)
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.99e-01

TTTTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGA  >  NC_000913/710296‑710572
                                                                                                                                            |                                                                                                                                        
ttttGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGc                                                                                                                                        >  2:680470/1‑143 (MQ=255)
ttttGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGc                                                                                                                                        >  1:404223/1‑143 (MQ=255)
  ttccGGCAATTTCACGCTCACGCGCTTCGATTGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCgg                                                                                                                                      <  2:671674/140‑1 (MQ=255)
       gCAAGTTCACGCTGACTGCCTTCGATGGAATCACCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCTTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGAc                                                                                                                                 <  2:417346/143‑1 (MQ=255)
         aaTTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGa                                                                                                                               >  1:94227/1‑143 (MQ=255)
          aTTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGAt                                                                                                                              <  1:54223/143‑1 (MQ=255)
              cgcGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCa                                                                                                                          <  2:539436/141‑1 (MQ=255)
              aaCGCAGACGCGCTTCGAGGGAATCAGCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCa                                                                                                                          <  1:119632/142‑1 (MQ=255)
                gtcTGACGCGCTTCGATGGAATGTGCACTAAATTCGATAGCCTTGCCGCAGTCTAGTCATATCAGGTGATCGTGCTGATGGCGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGc                                                                                                                        <  2:1171666/141‑1 (MQ=255)
                       gcgcTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCCGTTCAGATACGGATTTACCGCCTTCAACATTGTAGCGGGCGAGGATACCAGCgccgtcc                                                                                                                 >  1:110660/1‑142 (MQ=255)
                         gcTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCaaa                                                                                                               <  1:557087/143‑1 (MQ=255)
                          cTTCGAGGCAATCATCCCTAAATGCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAAc                                                                                                              <  1:941179/143‑1 (MQ=255)
                           ttCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGCTGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCACACt                                                                                                             >  1:839721/1‑143 (MQ=255)
                               aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAGATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGtt                                                                                                         >  2:56285/1‑143 (MQ=255)
                               aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGtt                                                                                                         >  2:60464/1‑143 (MQ=255)
                                   aaTCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAACCTCGTCGCGTGTGAGGATACCAGCGTCGTCAAACTGGTTCAGt                                                                                                     >  1:232797/1‑143 (MQ=255)
                                    atcatcACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTa                                                                                                    <  1:36933/143‑1 (MQ=255)
                                     tcatcaCTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTGACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGAGTCGTCAAACTGGTTCAGTAc                                                                                                   >  1:648690/1‑143 (MQ=255)
                                           gtgacttCGAGAGCCTGGCCGCAGTCGAGGAGTATCCGGTGTTCGTTGTGACGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTTTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGata                                                                                             <  2:532184/138‑1 (MQ=255)
                                           cTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGata                                                                                             >  2:805562/1‑143 (MQ=255)
                                                ttCGATAACCTTGCCGCAGTCGAGGCAGCTCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCACACTGGTTCAGTACGCGACATACCg                                                                                        >  1:935991/1‑143 (MQ=255)
                                                    aTAGCCTTCGCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGc                                                                                    <  2:1188073/143‑1 (MQ=255)
                                                        ccTTGGCGCAGGCGAGGGAGGTCGGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGa                                                                                <  2:659444/143‑1 (MQ=255)
                                                         cTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGAc                                                                               >  1:955098/1‑143 (MQ=255)
                                                                   tCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATttcttc                                                                     <  2:969267/143‑1 (MQ=255)
                                                                       ggCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGAAACGGGTGACGATACAAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCAccc                                                                 >  1:393728/1‑143 (MQ=255)
                                                                           gATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCatat                                                             >  2:917990/1‑143 (MQ=255)
                                                                             tCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCg                                                           >  1:1238105/1‑143 (MQ=255)
                                                                                ggTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATc                                                        <  1:1265357/143‑1 (MQ=255)
                                                                                 gTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCa                                                       <  2:404223/143‑1 (MQ=255)
                                                                                  tGATCGTGGTGATGTTGCGGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAg                                                      <  1:306677/143‑1 (MQ=255)
                                                                                          gtgATTGTGCGGGGTCAGTTCAAATACGTATGTACCGCCGTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGt                                              <  2:1007938/143‑1 (MQ=255)
                                                                                           tgATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGGTCAGACGTTTGta                                             >  1:743186/1‑143 (MQ=255)
                                                                                              tGTTGCTGTGTCAATTCGAATACGGAGTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATaa                                          <  2:466644/143‑1 (MQ=255)
                                                                                               gTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATaaa                                         >  1:1045872/1‑143 (MQ=255)
                                                                                                ttGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGACGTCAAACTGATACAGTACGCGATATACCGTAGCCAGACCAATATCTTAACCCAGATAGATCAGCCGTACGTACAAAt                                        >  2:1174008/1‑143 (MQ=255)
                                                                                                 tGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATc                                       >  1:65121/1‑143 (MQ=255)
                                                                                                  gCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCt                                      <  1:722060/143‑1 (MQ=255)
                                                                                                      tgtCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCg                                  <  2:955098/143‑1 (MQ=255)
                                                                                                        tCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGGATCAATCTTCCGCa                                >  2:1210356/1‑143 (MQ=255)
                                                                                                           gTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTg                             >  1:754958/1‑143 (MQ=255)
                                                                                                                   aCGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATg                     >  1:999332/1‑143 (MQ=255)
                                                                                                                      gATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCCCCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGtt                  >  1:741822/1‑143 (MQ=255)
                                                                                                                       aTTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGCGTGGATAAATATTCCACACTGAAGCGATgtgtg                 >  2:880549/1‑139 (MQ=255)
                                                                                                                        tttACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTAAGTAGGCGATATACTGTAGCCAGACCAATTTCTTCACCCATATCGATGAGACGTTACTAAAAATCTTCCACAATGACGTAATGGTTGt                >  2:1095382/1‑143 (MQ=255)
                                                                                                                            ccGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCgg            >  1:987760/1‑143 (MQ=255)
                                                                                                                            ccGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCgg            >  2:802050/1‑143 (MQ=255)
                                                                                                                             cGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGCTCAGACGTTTGTATAAATCTTCCGCACTGACGTGCGGGTTGTCCGGc           >  2:1064758/1‑143 (MQ=255)
                                                                                                                              gCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCt          <  1:1063917/143‑1 (MQ=255)
                                                                                                                                 ttCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAAACGTTTGTATAAATCTTCCGCACCGACGTGATGGTTGTCCGGCTCCt       >  2:42828/1‑143 (MQ=255)
                                                                                                                                 ttCAAAATTGTAGCGGGAGACGATACCAGCGTCGTCAAACTGGTTCAGCACGCGATATACCGTAGCCAGCCCAATTTCTTCACCCATATCGACCAGACGCTTGCATAAACCCTCCGCACCAACGCGATCGTTGCCCGGCGCCt       >  2:505076/1‑143 (MQ=255)
                                                                                                                                     acaTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag   <  1:165197/141‑1 (MQ=255)
                                                                                                                                      aaTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaaga  <  2:475317/143‑1 (MQ=255)
                                                                                                                                      aaTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaaga  >  1:697423/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                        
TTTTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGA  >  NC_000913/710296‑710572

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: ‑