Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,044,167 A→G L28P (CTG→CCG)  rfbA ← glucose‑1‑phosphate thymidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,044,1670AG99.4% 694.6 / ‑4.0 167L28P (CTG→CCG) rfbAglucose‑1‑phosphate thymidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (70/96);  total (71/96)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

Only 100 of 196 total aligned reads displayed.

AATATCGCGAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCAGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCTTTAC  >  REL606/2044094‑2044239
                                                                         |                                                                        
  tatCGCGAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGgttg                                                                      >  1:5387452/1‑76 (MQ=255)
   atCGCGAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGgttgt                                                                     >  1:296993/1‑76 (MQ=255)
   atCGCGAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGgttgt                                                                     >  1:2082245/1‑76 (MQ=255)
   atCGCGAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGgttgt                                                                     >  1:6960830/1‑76 (MQ=255)
      gcgAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTtta                                                                  <  1:775141/76‑1 (MQ=255)
        gAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTttact                                                                >  1:6325064/1‑76 (MQ=255)
        gAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTttact                                                                >  2:3356405/1‑76 (MQ=255)
        gAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTttact                                                                <  1:5630522/76‑1 (MQ=255)
         aaTACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTttactt                                                               <  2:4352656/76‑1 (MQ=255)
          aTACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTttactta                                                              <  1:124216/76‑1 (MQ=255)
            aCCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGGTGTTtaatttct                                                            >  2:1502994/1‑76 (MQ=255)
             cccGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGTGGTTGTTtacttactc                                                           >  2:6643906/1‑75 (MQ=255)
             cccGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTttacttatca                                                           <  2:6457612/76‑4 (MQ=255)
             cccGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTttcttactt                                                           >  2:240043/1‑75 (MQ=255)
             cccGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTg                                                           <  2:2357678/76‑1 (MQ=255)
              ccGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGGTTGATTTCTGc                                                          >  2:4570988/1‑76 (MQ=255)
               cGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTggc                                                         >  1:5824540/1‑74 (MQ=255)
               cGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGcc                                                         >  1:5410634/1‑76 (MQ=255)
               cGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGcc                                                         >  1:3345031/1‑76 (MQ=255)
               cGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGcc                                                         <  1:3379448/76‑1 (MQ=255)
               cGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGcc                                                         >  1:6040223/1‑76 (MQ=255)
                gCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCa                                                        <  2:509857/76‑1 (MQ=255)
                gCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCa                                                        >  1:6817057/1‑76 (MQ=255)
                 cTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCAt                                                       >  1:2844854/1‑76 (MQ=255)
                  tAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATc                                                      <  2:3299400/76‑1 (MQ=255)
                  tAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATc                                                      >  1:4265595/1‑76 (MQ=255)
                  tAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATc                                                      <  2:2512639/76‑1 (MQ=255)
                  tAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATc                                                      <  1:1602729/76‑1 (MQ=255)
                   aaCATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTTCTGCCATCg                                                     >  2:2466332/1‑76 (MQ=255)
                   aaCATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTTCCATCg                                                     >  2:7091444/1‑76 (MQ=255)
                   aaCATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCg                                                     >  1:7122484/1‑76 (MQ=255)
                   aaCATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCg                                                     >  1:120775/1‑76 (MQ=255)
                   aaCATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGGTGTTTACTTTCTGCCATCg                                                     >  2:3850943/1‑76 (MQ=255)
                    aCATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGt                                                    >  1:5956742/1‑76 (MQ=255)
                     cATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCCTCGTc                                                   >  2:3541213/1‑76 (MQ=255)
                      aTAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTca                                                  <  1:3934090/76‑1 (MQ=255)
                         aGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAg                                               <  1:2878730/76‑1 (MQ=255)
                         aGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAg                                               <  1:6200731/76‑1 (MQ=255)
                          gCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAgg                                              <  1:999036/76‑1 (MQ=255)
                          gCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAgg                                              <  1:4117576/76‑1 (MQ=255)
                            gTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGAt                                            <  1:4677064/76‑1 (MQ=255)
                             ttGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATa                                           <  2:3480676/76‑1 (MQ=255)
                             ttGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATa                                           <  1:6840141/76‑1 (MQ=255)
                               gAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATaaa                                         <  1:4299993/76‑1 (MQ=255)
                                aaaGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATaaaa                                        >  2:5812762/1‑75 (MQ=255)
                                aaaGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAg                                        >  1:3644568/1‑76 (MQ=255)
                                aaaGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAg                                        <  1:1725183/76‑1 (MQ=255)
                                aaaGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAg                                        >  1:2438143/1‑76 (MQ=255)
                                 aaGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGa                                       <  1:2796720/76‑1 (MQ=255)
                                  aGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGAc                                      <  2:6208835/76‑1 (MQ=255)
                                   gCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACg                                     >  1:1335554/1‑76 (MQ=255)
                                   gCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACg                                     <  1:3189531/76‑1 (MQ=255)
                                       ataataAATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTg                                 <  2:5122828/76‑1 (MQ=255)
                                          ataaATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCg                              <  2:5953044/76‑1 (MQ=255)
                                           taaATCATCGGCTTATCATAAATCGGCAGCGGTTGTTTACTTGCTGCCATCGTCACAGGATAAAGACGAGTGCCgg                             <  1:6980125/76‑1 (MQ=255)
                                              atcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTTCTGCCATCGTCACAGGATAAAGACGGGTGCCGGCac                          >  1:2256416/1‑76 (MQ=255)
                                              atcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAac                          <  2:3272362/76‑1 (MQ=255)
                                              atcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAac                          >  1:5211984/1‑76 (MQ=255)
                                              atcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAac                          >  1:5002906/1‑76 (MQ=255)
                                              atcatcGGCTTATCATAAATCGGCAGCGGGTGTTTACTTACTGCCATCGTCACAGGATAAAAACGAGGGGCGGGac                          >  2:1997915/1‑76 (MQ=255)
                                               tcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCCTCGTCACAGGATAAAGACGAGTGCCGGaaac                         >  2:5231866/1‑74 (MQ=255)
                                               tcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAacc                         <  2:2971344/76‑1 (MQ=255)
                                               tcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAacc                         >  2:2578832/1‑76 (MQ=255)
                                               tcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAacc                         <  2:844633/76‑1 (MQ=255)
                                               tcatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAacc                         <  1:4570462/76‑1 (MQ=255)
                                               acatcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACAGCCATCGTCACAGGATAAAGACGAGTGCCGGAacc                         <  2:6678293/75‑1 (MQ=255)
                                                catcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAacca                        <  2:4137877/76‑1 (MQ=255)
                                                  tcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAaccacc                      <  2:2744012/76‑1 (MQ=255)
                                                  tcGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAaccacc                      <  2:6135982/76‑1 (MQ=255)
                                                   cGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAAccacca                     <  2:4705684/76‑1 (MQ=255)
                                                   cGGCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAAccacca                     >  1:3719838/1‑76 (MQ=255)
                                                    ggCTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCAccag                    <  2:6928358/76‑1 (MQ=255)
                                                      cTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCAccagcc                  <  2:1891448/76‑1 (MQ=255)
                                                      cTTATCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCAccagcc                  <  2:1609604/76‑1 (MQ=255)
                                                          tCATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaa              >  1:5321446/1‑76 (MQ=255)
                                                           cATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaat             >  1:3597870/1‑76 (MQ=255)
                                                           cATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaat             >  1:6639102/1‑76 (MQ=255)
                                                           cATAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaat             >  1:5141503/1‑76 (MQ=255)
                                                            aaaaaTCGGCAGCGGTTGTTTCCTTGCTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaata            <  2:4103921/74‑1 (MQ=255)
                                                            aTAAATCGGCAGCGTTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaata            <  1:5569339/76‑1 (MQ=255)
                                                            aTAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaata            <  1:4486042/76‑1 (MQ=255)
                                                             tAAATCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaataa           <  2:3812308/76‑1 (MQ=255)
                                                              aaaTCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGTATAAAGACGAGTGCCGGAACCACCAGCCAGaataat          >  1:5663623/1‑76 (MQ=255)
                                                              aaaTCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaataat          <  1:2415571/76‑1 (MQ=255)
                                                              aaaTCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaataat          >  1:4935232/1‑76 (MQ=255)
                                                              aaaTCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGaataat          <  1:7106908/76‑1 (MQ=255)
                                                               aaTCGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCCAAataata         >  2:4795273/1‑76 (MQ=255)
                                                                  cGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCt      >  2:1953963/1‑76 (MQ=255)
                                                                  cGGCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCt      >  1:429118/1‑76 (MQ=255)
                                                                   ggCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCtt     >  2:5482462/1‑76 (MQ=255)
                                                                   ggCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCtt     <  2:3904462/76‑1 (MQ=255)
                                                                   ggCAGCGGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCtt     <  2:2934094/76‑1 (MQ=255)
                                                                    gcagcgGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCttt    <  2:516447/76‑1 (MQ=255)
                                                                    gcagcgGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCttt    >  2:4436097/1‑76 (MQ=255)
                                                                    gcagcgGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCttt    >  1:7013509/1‑76 (MQ=255)
                                                                    gcagcgGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCttt    <  1:1984059/76‑1 (MQ=255)
                                                                    gcagcgGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCttt    <  1:1111509/76‑1 (MQ=255)
                                                                     cagcggTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCTTTa   >  1:5738633/1‑76 (MQ=255)
                                                                     cagcggTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCTTTa   >  1:2658098/1‑76 (MQ=255)
                                                                     cagcggTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCTTTa   <  1:1031103/76‑1 (MQ=255)
                                                                         |                                                                        
AATATCGCGAATACCCGCTAACATAAGCGTTGAAAGCGGATAATAAATCATCGGCTTATCATAAATCGGCAGCAGTTGTTTACTTACTGCCATCGTCACAGGATAAAGACGAGTGCCGGAACCACCAGCCAGAATAATACCTTTAC  >  REL606/2044094‑2044239

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