Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 711,161 C→T 13.7% intergenic (‑93/+196) fur ← / ← fldA DNA‑binding transcriptional dual regulator/flavodoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048711,1610CT13.7% 138.7 / 11.6 51intergenic (‑93/+196)fur/fldADNA‑binding transcriptional dual regulator/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base C (18/26);  new base T (2/5);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAAC  >  AP009048/711097‑711231
                                                                |                                                                      
cACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTa                                                                   <  1:351494/70‑1 (MQ=255)
  cTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAAt                                                                >  1:824989/1‑71 (MQ=255)
     cATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATg                                                               <  1:4208824/69‑1 (MQ=255)
     cATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGa                                                              >  1:22351/1‑70 (MQ=255)
      aTTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaa                                                            <  1:4465391/71‑1 (MQ=255)
       ttAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaaa                                                           >  1:3656884/1‑71 (MQ=255)
         agaagaAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAAc                                                          >  1:4514626/1‑70 (MQ=255)
                  gACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGa                                                              <  1:4322025/57‑1 (MQ=255)
                   aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTccc                                                <  1:429856/70‑1 (MQ=255)
                   aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaaa                                                                 >  1:2855099/1‑53 (MQ=255)
                    cAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTccc                                                <  1:1280913/69‑1 (MQ=255)
                      gAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTgg                                             >  1:4266469/1‑70 (MQ=255)
                        aTTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc                                                                        >  1:3439759/1‑41 (MQ=255)
                           tGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTc                                                  >  1:904217/1‑60 (MQ=255)
                             ctctTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAatat                                      <  1:3779806/70‑1 (MQ=255)
                              tctTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTccc                                                <  1:3029215/59‑1 (MQ=255)
                               ctTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCt                                               <  1:3918144/59‑1 (MQ=255)
                                ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaa                                                            <  1:435707/45‑1 (MQ=255)
                                ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaa                                                            <  1:2554206/45‑1 (MQ=255)
                                ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACcac                                                       <  1:957054/50‑1 (MQ=255)
                                ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTg                                    <  1:1876226/69‑1 (MQ=255)
                                ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGt                                   >  1:644259/1‑70 (MQ=255)
                                ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGt                                   >  1:880753/1‑70 (MQ=255)
                                  gagaTAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACcac                                                       >  1:1382428/1‑48 (MQ=255)
                                  gagaTAATGCGTATCATTATAGAATTGCCACGCCTAAAt                                                                >  1:217369/1‑39 (MQ=255)
                                   agaTAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaa                                >  1:3430533/1‑70 (MQ=255)
                                   agaTAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaa                                >  1:2647644/1‑70 (MQ=255)
                                     aTAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaataa                             <  1:602390/71‑1 (MQ=255)
                                     aTAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaataa                             <  1:3081043/71‑1 (MQ=255)
                                      tAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAatat                                      <  1:4205212/61‑1 (MQ=255)
                                      tAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaa                            <  1:649565/71‑1 (MQ=255)
                                      tAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaa                            <  1:2255389/71‑1 (MQ=255)
                                       aaTGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACaa                                                     >  1:3078717/1‑45 (MQ=255)
                                       aaTGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAata                                       <  1:553777/59‑1 (MQ=255)
                                         tGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAAc                         <  1:3057775/71‑1 (MQ=255)
                                              aTCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTgg                     >  1:3026248/1‑70 (MQ=255)
                                               tCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGca                   <  1:2072511/71‑1 (MQ=255)
                                               tCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGca                   <  1:2243379/71‑1 (MQ=255)
                                                cATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATAtt                                     <  1:1938202/52‑1 (MQ=255)
                                                   tataGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTg                      >  1:549432/1‑64 (MQ=255)
                                                    ataGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCt              >  1:1049716/1‑71 (MQ=255)
                                                       gAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTa            <  1:105370/70‑1 (MQ=255)
                                                       gAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTAc           <  1:3855468/71‑1 (MQ=255)
                                                        aaTTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATAtt                                     >  1:2433384/1‑44 (MQ=255)
                                                         aTTGCCACGCCTAAATGAAAACCACAAGCCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACa          <  1:1979981/70‑1 (MQ=255)
                                                          ttGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAata                                       >  1:1108196/1‑40 (MQ=255)
                                                            gCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACAttt       <  1:1660609/70‑1 (MQ=255)
                                                             ccATGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTAc     <  1:3507686/71‑1 (MQ=255)
                                                             ccACGCCTAAATGAAAACCACTAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTAc     <  1:3110184/71‑1 (MQ=255)
                                                             ccACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGc                    <  1:4254604/56‑1 (MQ=255)
                                                               aCGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTta                                 >  1:1680846/1‑41 (MQ=255)
                                                               aCGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACaa   <  1:2071607/71‑1 (MQ=255)
                                                                cGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaataa                             >  1:3144645/1‑44 (MQ=255)
                                                                cGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAAc  >  1:4055013/1‑71 (MQ=255)
                                                                |                                                                      
CACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAAC  >  AP009048/711097‑711231

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