Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,261,861 G→A F103F (TTC→TTT hflB ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,261,8610GA99.2% 500.8 / ‑3.6 123F103F (TTC→TTThflBprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (69/53);  total (70/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

Only 100 of 134 total aligned reads displayed.

CACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTA  >  REL606/3261788‑3261933
                                                                          |                                                                        
ccccgccgCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATaaa                                                                        <  1:2900732/74‑1 (MQ=255)
 aCCGCCGCCCTGCATTTGCCGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAg                                                                       <  2:4723678/76‑1 (MQ=255)
 aCCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAg                                                                       <  1:3095056/76‑1 (MQ=255)
 aCCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAg                                                                       >  1:2669590/1‑76 (MQ=255)
 aCCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAg                                                                       >  2:4250170/1‑76 (MQ=255)
 aCCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAg                                                                       >  1:4062453/1‑76 (MQ=255)
  ccgccgCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGa                                                                      >  1:3960541/1‑76 (MQ=255)
   cgccgcCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGAt                                                                     <  1:1820431/76‑1 (MQ=255)
   cgccgcCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGAt                                                                     <  1:2379830/76‑1 (MQ=255)
    gccgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATa                                                                    <  1:3526704/76‑1 (MQ=255)
    gccgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATa                                                                    <  2:1972522/76‑1 (MQ=255)
     ccgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAg                                                                   >  1:596768/1‑76 (MQ=255)
     ccgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAg                                                                   <  1:2198609/76‑1 (MQ=255)
      cgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGa                                                                  >  2:2537577/1‑76 (MQ=255)
      cgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGa                                                                  >  1:3366181/1‑76 (MQ=255)
       gccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGaa                                                                 <  1:2232675/76‑1 (MQ=255)
       gccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGaa                                                                 <  1:3352636/76‑1 (MQ=255)
        cccTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAg                                                                >  1:1381567/1‑76 (MQ=255)
          cTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGcc                                                              >  2:1622395/1‑76 (MQ=255)
            gCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCcag                                                            >  2:1995825/1‑76 (MQ=255)
            gCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCcag                                                            <  1:2845628/76‑1 (MQ=255)
              aTTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCcagca                                                          >  2:271916/1‑76 (MQ=255)
              aTTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCcagca                                                          >  2:1749194/1‑76 (MQ=255)
               tttGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCcagcag                                                         >  1:897088/1‑76 (MQ=255)
                ttGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAgg                                                        <  1:3826857/76‑1 (MQ=255)
                 tGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGc                                                       >  2:2544337/1‑76 (MQ=255)
                  gACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCt                                                      >  1:4776869/1‑76 (MQ=255)
                  gACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCt                                                      >  1:1239506/1‑76 (MQ=255)
                    cGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTg                                                    >  1:2973783/1‑76 (MQ=255)
                    cGCAGGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTg                                                    <  2:1807290/76‑1 (MQ=255)
                     gCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTgg                                                   <  1:1935572/76‑1 (MQ=255)
                       aTGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGtt                                                 >  2:4477978/1‑76 (MQ=255)
                       aTGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGtt                                                 >  2:3770964/1‑76 (MQ=255)
                        tGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGttc                                                >  1:2601628/1‑76 (MQ=255)
                        tGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGttc                                                <  2:3633649/76‑1 (MQ=255)
                         gaagaaGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGttct                                               >  1:626034/1‑76 (MQ=255)
                           agaagaTCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGttcttc                                             >  2:2292418/1‑76 (MQ=255)
                            gaagaTCCCGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCa                                            <  1:4102643/76‑1 (MQ=255)
                            gaaaaTCCAGAACCCATGAAGAAACAACGTCGGGAACCAGGAGATAAAGATAAAAGCCAGCAGGCTTGGTTCTTCa                                            <  2:3806922/76‑1 (MQ=255)
                             aagaTCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAg                                           <  2:118664/76‑1 (MQ=255)
                                  ccAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGt                                      >  2:4004612/1‑76 (MQ=255)
                                  ccAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGt                                      >  2:1346603/1‑76 (MQ=255)
                                   cAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGtt                                     >  1:1600789/1‑76 (MQ=255)
                                   cAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGtt                                     >  1:1247008/1‑76 (MQ=255)
                                    aGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTc                                    <  1:1983444/76‑1 (MQ=255)
                                    aGACACCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTc                                    <  1:3160451/76‑1 (MQ=255)
                                      acacCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCAc                                  <  2:4833893/76‑1 (MQ=255)
                                       cacCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCAcc                                 <  2:2396125/76‑1 (MQ=255)
                                        acCAATCAGCAACAGCATGGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCg                                <  2:300504/76‑1 (MQ=255)
                                        acCAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCg                                <  1:2310037/76‑1 (MQ=255)
                                         ccAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGa                               <  2:3751051/76‑1 (MQ=255)
                                         ccAATCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGa                               >  2:1717041/1‑76 (MQ=255)
                                           aaccAGAAACAGCATCGGAAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACa                             <  1:3112936/73‑1 (MQ=255)
                                           aaTCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACa                             <  1:1983064/76‑1 (MQ=255)
                                            aTCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTTACCGACaa                            >  1:2126851/1‑76 (MQ=255)
                                            aTCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACaa                            >  2:2890677/1‑76 (MQ=255)
                                            aTCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACaa                            <  2:3058911/76‑1 (MQ=255)
                                            aTCAGCACAGCCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACaa                            <  2:330590/76‑1 (MQ=255)
                                             tCAGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAAc                           <  2:182359/76‑1 (MQ=255)
                                              cAGAAACAGCATCGGGAACCAGGAGATAAGGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAAcc                          <  2:1227766/76‑1 (MQ=255)
                                               agaaaCAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCt                         <  1:4180391/73‑1 (MQ=255)
                                               aGCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCt                         >  2:4057387/1‑76 (MQ=255)
                                                gaaaCAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCtt                        <  1:687129/74‑1 (MQ=255)
                                                gCAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCtt                        <  1:3500335/76‑1 (MQ=255)
                                                 cAACAGCATGGGAAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTg                       <  1:3733280/76‑1 (MQ=255)
                                                 cAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTg                       >  1:878312/1‑76 (MQ=255)
                                                 cAACAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTg                       >  1:4039332/1‑76 (MQ=255)
                                                  aaCAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGa                      >  2:1042140/1‑76 (MQ=255)
                                                  aaCAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGa                      >  2:4709794/1‑76 (MQ=255)
                                                  aaCAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGa                      >  1:4302741/1‑76 (MQ=255)
                                                  aaCAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGa                      >  2:243120/1‑76 (MQ=255)
                                                    cAGCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACg                    <  1:4925670/76‑1 (MQ=255)
                                                      gCATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGtt                  <  2:4515775/76‑1 (MQ=255)
                                                        cATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTc                 >  1:2409924/1‑76 (MQ=255)
                                                        cATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTc                 >  1:3001986/1‑76 (MQ=255)
                                                        cATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTc                 >  2:2804707/1‑76 (MQ=255)
                                                        cATCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTc                 >  1:1563132/1‑76 (MQ=255)
                                                         aTCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACGTTGACGTTCt                >  2:2990946/1‑76 (MQ=255)
                                                         aTCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCt                >  2:4943691/1‑76 (MQ=255)
                                                         aTCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCt                >  2:4416675/1‑76 (MQ=255)
                                                          tCGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCtt               <  2:3468848/76‑1 (MQ=255)
                                                           cGGGAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTg              >  2:2048945/1‑76 (MQ=255)
                                                            gggAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTgg             <  1:756331/76‑1 (MQ=255)
                                                            gggAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTgg             <  2:4560931/76‑1 (MQ=255)
                                                            gggAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTgg             >  2:2588261/1‑76 (MQ=255)
                                                              gAACCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTc           <  1:4444881/76‑1 (MQ=255)
                                                                aCCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCaa         >  1:4286605/1‑76 (MQ=255)
                                                                aCCAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCaa         <  2:2836900/76‑1 (MQ=255)
                                                                 ccAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAAc        <  1:417157/76‑1 (MQ=255)
                                                                 ccAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAAc        <  1:1727947/76‑1 (MQ=255)
                                                                  cAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACa       >  1:1747405/1‑76 (MQ=255)
                                                                  cAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACa       >  2:2901920/1‑76 (MQ=255)
                                                                  cAGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACa       <  2:1560887/76‑1 (MQ=255)
                                                                   aGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAg      <  1:358308/76‑1 (MQ=255)
                                                                   aGGAGATAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAg      >  1:3686244/1‑76 (MQ=255)
                                                                     gagaTAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGt    >  2:2698375/1‑76 (MQ=255)
                                                                      agaTAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGtt   <  1:3410492/76‑1 (MQ=255)
                                                                      agaTAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGtt   >  1:1170147/1‑76 (MQ=255)
                                                                      agaTAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGtt   >  1:4295539/1‑76 (MQ=255)
                                                                       gaTAAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTa  >  1:729243/1‑76 (MQ=255)
                                                                          |                                                                        
CACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTA  >  REL606/3261788‑3261933

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