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

ATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTC  >  NC_000913/710779‑710970
                                                                                                    |                                                                                             
ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAatta                                                                                                >  2:1547617/1‑100 (MQ=255)
     aaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAttttt                                                                                         <  1:930272/100‑2 (MQ=255)
      aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttttg                                                                                        <  2:1400820/100‑3 (MQ=255)
        gTGGCATAGCTTACAGTGAGTACGGCGAGGTGACTTGAACATCTGGATATTGTACAGGTACAACAGCAAAAGTCACAAATTTGCAGCAATTAgttgt                                                                                         <  2:1313704/97‑3 (MQ=255)
         ttgCACAGCTTACATTTACAACGGCAATGCCACTTGAACACCGTGACATTGTACAGGTACTACTGCAAAAGTTACAAAGTTGTAGCAATTAtttttgat                                                                                      <  2:458258/97‑5 (MQ=255)
         tGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTAtttttgatt                                                                                     <  1:1408666/100‑6 (MQ=255)
          ggCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGAAAAAGTTACAAATTTGTAGCAATTATTTTTGATTg                                                                                     >  1:395835/1‑100 (MQ=255)
           gCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTgg                                                                                    >  1:996960/1‑100 (MQ=255)
           gCACACCTTACATTTACAACAGCAAGGCTACGTAAACACAGGGACATTGTACAGGTACAACAGCAAAAGGTACAAATTTGTAGCAATTAtttttgatt                                                                                     <  2:1538847/98‑6 (MQ=255)
               aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCAtt                                                                                <  2:729535/100‑1 (MQ=255)
                gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTa                                                                               >  1:1152743/1‑100 (MQ=255)
                  ttACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATc                                                                             <  1:367005/100‑1 (MQ=255)
                     cATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTAt                                                                          >  1:1119804/1‑100 (MQ=255)
                      aTTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTAtt                                                                         <  1:340082/100‑1 (MQ=255)
                       tttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTa                                                                        <  2:468227/100‑1 (MQ=255)
                       tttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTa                                                                        <  1:142912/100‑1 (MQ=255)
                        ttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGTTACAACAGCATAAGTGACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTaa                                                                       <  2:228254/100‑1 (MQ=255)
                         tACAACGGCAAGGCCACTTGAACTCCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAAt                                                                      >  2:1520956/1‑100 (MQ=255)
                         tACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAAt                                                                      >  1:1361166/1‑100 (MQ=255)
                           cAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATAc                                                                    >  2:1236592/1‑100 (MQ=255)
                             aCGGCAAGGCCACTTGAACACCGGGACATTGTACAGTTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACgg                                                                  <  2:279783/100‑1 (MQ=255)
                             aCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACgg                                                                  <  1:549988/100‑1 (MQ=255)
                               ggCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCg                                                                <  2:31531/100‑1 (MQ=255)
                               ggCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCg                                                                <  2:3821/100‑1 (MQ=255)
                                  aaGGCCACTTGAACTCCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAg                                                             <  1:1520966/100‑1 (MQ=255)
                                   aGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGa                                                            >  2:1379231/1‑100 (MQ=255)
                                    ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGAc                                                           <  2:820298/100‑1 (MQ=255)
                                      ccACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATg                                                        <  1:1300132/100‑1 (MQ=255)
                                       cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATg                                                        <  2:618941/100‑1 (MQ=255)
                                        aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGa                                                       >  1:1093515/1‑100 (MQ=255)
                                        aCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                      <  1:49594/100‑1 (MQ=255)
                                         cTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                      <  1:829366/100‑1 (MQ=255)
                                         cTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                      <  1:1420797/100‑1 (MQ=255)
                                         cTTGAACACCGGGACATTGTACAGGCACAACAGCAAACGTTACCAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACCGCGTAGACATg                                                        >  1:519663/1‑98 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    <  2:1119797/100‑1 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    >  1:1498501/1‑100 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    >  2:1108296/1‑100 (MQ=255)
                                           tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTc                                                    >  1:877430/1‑100 (MQ=255)
                                               cacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcg                                                >  2:1475364/1‑100 (MQ=255)
                                                acCGGGACATTGTACGGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgc                                               <  1:1547627/100‑1 (MQ=255)
                                                 ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                              >  2:597800/1‑100 (MQ=255)
                                                   gggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgc                                               >  2:1139567/1‑97 (MQ=255)
                                                   gggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                           >  2:922145/1‑100 (MQ=255)
                                                    ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCcaca                                           <  1:1355926/100‑1 (MQ=255)
                                                    ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                           <  1:1424996/100‑1 (MQ=255)
                                                    ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcat                                          >  2:901597/1‑100 (MQ=255)
                                                     gACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcat                                          >  2:1208734/1‑100 (MQ=255)
                                                         ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcaca                                      <  1:1121675/100‑1 (MQ=255)
                                                         ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcaca                                      >  2:211130/1‑100 (MQ=255)
                                                           gTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcacatc                                    >  2:1341730/1‑100 (MQ=255)
                                                            tACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcacatc                                    >  2:209571/1‑99 (MQ=255)
                                                              naGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAgg                                 >  2:690683/2‑100 (MQ=255)
                                                                ggTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCa                               <  1:1130314/100‑1 (MQ=255)
                                                                 gTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAt                              >  2:184117/1‑100 (MQ=255)
                                                                 gTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTATTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAt                              >  2:843094/1‑100 (MQ=255)
                                                                  tACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATCCGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAt                              >  1:356244/1‑99 (MQ=255)
                                                                   acaacaGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                            <  2:996952/100‑1 (MQ=255)
                                                                   acaacaGAAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                            <  2:961510/100‑1 (MQ=255)
                                                                    caacaGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga                           >  1:1029826/1‑100 (MQ=255)
                                                                     aacaGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTgag                          <  2:1404184/100‑1 (MQ=255)
                                                                        aGCAAAAGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAt                       >  2:295524/1‑100 (MQ=255)
                                                                            aaaGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                   <  1:568318/100‑1 (MQ=255)
                                                                            aaaGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                   <  1:778553/100‑1 (MQ=255)
                                                                            aaaGTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                   <  2:1014628/100‑1 (MQ=255)
                                                                               gTTACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtcg                >  2:3829/1‑100 (MQ=255)
                                                                                 tACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag              <  2:1498491/100‑1 (MQ=255)
                                                                                 tACAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag              <  1:843100/100‑1 (MQ=255)
                                                                                   cAAATTTGTAGCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGAt            <  1:716029/100‑1 (MQ=255)
                                                                                     aaTATGTAGCAATTATTTTTGAATGGCATTAACTATTAATTCGGCGTATATATGAGTCTATGCCTCATCACATCAGGCATTGAGAATTTCGTCGAGATGc          >  2:1110153/1‑100 (MQ=255)
                                                                                             gCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc  <  1:1522626/100‑1 (MQ=255)
                                                                                             gCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc  >  2:411163/1‑100 (MQ=255)
                                                                                             gCAATTATTTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc  >  1:1363653/1‑100 (MQ=255)
                                                                                                    |                                                                                             
ATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTC  >  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: