Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 710,882 +T intergenic (‑157/+53) uof ← / ← fldA ryhB‑regulated fur leader peptide/flavodoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,8781.T92.3% 214.4 / 11.5 65intergenic (‑153/+57)uof/fldAryhB‑regulated fur leader peptide/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base T (29/31);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

ATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTC  >  NC_000913/710779‑710970
                                                                                                    |                                                                                             
ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAatta                                                                                                >  2:1547617/1‑100 (MQ=255)
     aaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑tt‑ttt                                                                                         <  1:930272/100‑2 (MQ=255)
      aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑tt‑tttg                                                                                        <  2:1400820/100‑3 (MQ=255)
        gTGGCATAGCTTACAGTGAGTACGGCGAGGTGACTTGAACATCTGGATATTGTACAGGTACAACAGCAAAAGTCACAAATTTGCAGCAATTA‑gt‑tgt                                                                                         <  2:1313704/97‑3 (MQ=255)
         ttgCACAGCTTACATTTACAACGGCAATGCCACTTGAACACCGTGACATTGTACAGGTACTACTGCAAAAGTTACAAAGTTGTAGCAATTA‑tt‑tttgat                                                                                      <  2:458258/97‑5 (MQ=255)
         tGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑tt‑tttgatt                                                                                     <  1:1408666/100‑6 (MQ=255)
          ggCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGAAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTg                                                                                     >  1:395835/1‑100 (MQ=255)
           gCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTgg                                                                                    >  1:996960/1‑100 (MQ=255)
           gCACACCTTACATTTACAACAGCAAGGCTACGTAAACACAGGGACATTGTACAGGTACAACAGCAAAAGGTACAAATTTGTAGCAATTA‑tt‑tttgatt                                                                                     <  2:1538847/98‑6 (MQ=255)
               aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCAtt                                                                                <  2:729535/100‑1 (MQ=255)
                gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTa                                                                               >  1:1152743/1‑100 (MQ=255)
                  ttACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATc                                                                             <  1:367005/100‑1 (MQ=255)
                     cATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTAt                                                                          >  1:1119804/1‑100 (MQ=255)
                      aTTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTAtt                                                                         <  1:340082/100‑1 (MQ=255)
                       tttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTa                                                                        <  2:468227/100‑1 (MQ=255)
                       tttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTa                                                                        <  1:142912/100‑1 (MQ=255)
                        ttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGTTACAACAGCATAAGTGACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTaa                                                                       <  2:228254/100‑1 (MQ=255)
                         tACAACGGCAAGGCCACTTGAACTCCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAAt                                                                      >  2:1520956/1‑100 (MQ=255)
                         tACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAAt                                                                      >  1:1361166/1‑100 (MQ=255)
                           cAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATAc                                                                    >  2:1236592/1‑100 (MQ=255)
                             aCGGCAAGGCCACTTGAACACCGGGACATTGTACAGTTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACgg                                                                  <  2:279783/100‑1 (MQ=255)
                             aCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACgg                                                                  <  1:549988/100‑1 (MQ=255)
                               ggCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCg                                                                <  2:31531/100‑1 (MQ=255)
                               ggCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCg                                                                <  2:3821/100‑1 (MQ=255)
                                  aaGGCCACTTGAACTCCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAg                                                             <  1:1520966/100‑1 (MQ=255)
                                   aGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGa                                                            >  2:1379231/1‑100 (MQ=255)
                                    ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGAc                                                           <  2:820298/100‑1 (MQ=255)
                                      ccACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATg                                                        <  1:1300132/100‑1 (MQ=255)
                                       cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATg                                                        <  2:618941/100‑1 (MQ=255)
                                        aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGa                                                       >  1:1093515/1‑100 (MQ=255)
                                        aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                      <  1:49594/100‑1 (MQ=255)
                                         cTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                      <  1:829366/100‑1 (MQ=255)
                                         cTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                      <  1:1420797/100‑1 (MQ=255)
                                         cTTGAACACCGGGACATTGTACAGGCACAACAGCAAACGTTACCAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACCGCGTAGACATg                                                        >  1:519663/1‑98 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    <  2:1119797/100‑1 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    >  1:1498501/1‑100 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    >  2:1108296/1‑100 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    >  1:877430/1‑100 (MQ=255)
                                               cacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcg                                                >  2:1475364/1‑100 (MQ=255)
                                                acCGGGACATTGTACGGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgc                                               <  1:1547627/100‑1 (MQ=255)
                                                 ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                              >  2:597800/1‑100 (MQ=255)
                                                   gggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgc                                               >  2:1139567/1‑97 (MQ=255)
                                                   gggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                           >  2:922145/1‑100 (MQ=255)
                                                    ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCcaca                                           <  1:1355926/100‑1 (MQ=255)
                                                    ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                           <  1:1424996/100‑1 (MQ=255)
                                                    ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcat                                          >  2:901597/1‑100 (MQ=255)
                                                     gACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcat                                          >  2:1208734/1‑100 (MQ=255)
                                                         ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcaca                                      <  1:1121675/100‑1 (MQ=255)
                                                         ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcaca                                      >  2:211130/1‑100 (MQ=255)
                                                           gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcacatc                                    >  2:1341730/1‑100 (MQ=255)
                                                            tACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcacatc                                    >  2:209571/1‑99 (MQ=255)
                                                              naGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAgg                                 >  2:690683/2‑100 (MQ=255)
                                                                ggTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCa                               <  1:1130314/100‑1 (MQ=255)
                                                                 gTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAt                              >  2:184117/1‑100 (MQ=255)
                                                                 gTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TTATTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAt                              >  2:843094/1‑100 (MQ=255)
                                                                  tACAACAGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATCCGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAt                              >  1:356244/1‑99 (MQ=255)
                                                                   acaacaGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                            <  2:996952/100‑1 (MQ=255)
                                                                   acaacaGAAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                            <  2:961510/100‑1 (MQ=255)
                                                                    caacaGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga                           >  1:1029826/1‑100 (MQ=255)
                                                                     aacaGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTgag                          <  2:1404184/100‑1 (MQ=255)
                                                                        aGCAAAAGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAt                       >  2:295524/1‑100 (MQ=255)
                                                                            aaaGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                   <  1:568318/100‑1 (MQ=255)
                                                                            aaaGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                   <  1:778553/100‑1 (MQ=255)
                                                                            aaaGTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                   <  2:1014628/100‑1 (MQ=255)
                                                                               gTTACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtcg                >  2:3829/1‑100 (MQ=255)
                                                                                 tACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag              <  2:1498491/100‑1 (MQ=255)
                                                                                 tACAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag              <  1:843100/100‑1 (MQ=255)
                                                                                   cAAATTTGTAGCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGAt            <  1:716029/100‑1 (MQ=255)
                                                                                     aaTATGTAGCAATTATTT‑TTGAATGGCATTAACTATTAATTCGGCGTATATATGAGTCTATGCCTCATCACATCAGGCATTGAGAATTTCGTCGAGATGc          >  2:1110153/1‑100 (MQ=255)
                                                                                             gCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc  <  1:1522626/100‑1 (MQ=255)
                                                                                             gCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc  >  2:411163/1‑100 (MQ=255)
                                                                                             gCAATTATTT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc  >  1:1363653/1‑100 (MQ=255)
                                                                                                    |                                                                                             
ATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTA‑TT‑TTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTC  >  NC_000913/710779‑710970

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