Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,504,870 T→C S129S (AGT→AGC)  ydcD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,504,8700TC93.0% 432.5 / 15.7 116S129S (AGT→AGC) ydcDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base C (64/43);  total (65/51)
Fisher's exact test for biased strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 118 total aligned reads displayed.

GTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGTTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGTTTTT  >  REL606/1504797‑1504942
                                                                         |                                                                        
gTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGgctagcta                                                                        >  2:2221134/1‑76 (MQ=255)
gTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGgctagcta                                                                        <  2:3855635/76‑1 (MQ=255)
 tAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtat                                                                       >  1:2590391/1‑76 (MQ=255)
 tAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtat                                                                       <  1:1438544/76‑1 (MQ=255)
 tAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtat                                                                       >  1:1698/1‑76 (MQ=255)
 tAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtat                                                                       >  1:1967189/1‑76 (MQ=255)
  aagcTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtata                                                                      <  2:5599120/73‑1 (MQ=255)
  aaaCTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtata                                                                      >  2:2201957/1‑76 (MQ=255)
  aaaCTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCtata                                                                      <  1:2439447/76‑1 (MQ=255)
    aCTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATaaa                                                                    >  1:4396609/1‑76 (MQ=255)
    aCTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATaaa                                                                    <  1:3353100/76‑1 (MQ=255)
     cTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATaaaa                                                                   >  1:4489139/1‑76 (MQ=255)
     cTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATaaaa                                                                   >  1:5529878/1‑76 (MQ=255)
     cTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATaaaa                                                                   >  1:896674/1‑76 (MQ=255)
      tGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaat                                                                  >  1:5343261/1‑76 (MQ=255)
      tGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaat                                                                  <  2:852498/76‑1 (MQ=255)
      tGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaat                                                                  <  2:896674/76‑1 (MQ=255)
       ggAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaata                                                                 >  2:6177563/1‑76 (MQ=255)
       ggAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaata                                                                 >  1:1227917/1‑76 (MQ=255)
       ggAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaata                                                                 <  1:964382/76‑1 (MQ=255)
       ggAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGATATAAaata                                                                 <  1:5823343/76‑1 (MQ=255)
         agaATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaataat                                                               <  2:4306360/76‑1 (MQ=255)
         agaATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAaataat                                                               >  2:3399177/1‑76 (MQ=255)
          gaATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAAtaatt                                                              <  2:6056029/76‑1 (MQ=255)
          gaATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAAtaatt                                                              <  1:4783670/76‑1 (MQ=255)
           aaTATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAAtaatta                                                             <  2:716642/76‑1 (MQ=255)
            atatGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAAtaattaa                                                            <  2:5521419/76‑1 (MQ=255)
              ctGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGTTATAAAATaattaatt                                                          <  1:6207627/75‑1 (MQ=255)
              atGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATaattaatt                                                          <  2:2009777/76‑1 (MQ=255)
              atGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATaattaatt                                                          >  1:3402983/1‑76 (MQ=255)
               tGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAattaatta                                                         <  1:5486000/76‑1 (MQ=255)
               tGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAattaatta                                                         <  2:4829364/76‑1 (MQ=255)
               tGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAattaatta                                                         >  1:4010136/1‑76 (MQ=255)
                 gAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATg                                                       >  2:5287436/1‑76 (MQ=255)
                   gTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGTTATAAAATAATTAATTATGga                                                     <  1:416209/76‑1 (MQ=255)
                   gTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGga                                                     <  2:2926044/76‑1 (MQ=255)
                    ttCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGgag                                                    >  2:2228823/1‑76 (MQ=255)
                        gcgcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACta                                                >  1:1187174/1‑76 (MQ=255)
                        gcgcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACta                                                >  1:4791971/1‑76 (MQ=255)
                         cgcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtat                                               >  2:4073272/1‑76 (MQ=255)
                          gcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatt                                              >  2:4476139/1‑76 (MQ=255)
                          gcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatt                                              >  2:4412142/1‑76 (MQ=255)
                          gcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatt                                              >  1:2492806/1‑76 (MQ=255)
                          gcTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatt                                              >  1:2895907/1‑76 (MQ=255)
                           cTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatta                                             >  1:5048309/1‑76 (MQ=255)
                           cTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatta                                             >  2:2600759/1‑76 (MQ=255)
                           cTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtatta                                             >  1:6296404/1‑76 (MQ=255)
                            tACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtattat                                            <  1:4491670/76‑1 (MQ=255)
                            tACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtattat                                            >  2:4690928/1‑76 (MQ=255)
                            tACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACtattat                                            <  1:18051/76‑1 (MQ=255)
                             aCAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTatg                                           >  1:2902850/1‑76 (MQ=255)
                             aCAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTatg                                           <  1:3033497/76‑1 (MQ=255)
                               aaaaCAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTatgat                                         >  2:4869813/1‑76 (MQ=255)
                                 aacaacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGt                                       >  2:2540595/1‑76 (MQ=255)
                                 aacaacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGt                                       >  1:5976237/1‑76 (MQ=255)
                                 aacaacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGt                                       >  1:585563/1‑76 (MQ=255)
                                 aacaacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGt                                       >  2:4739001/1‑76 (MQ=255)
                                  acaacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGtt                                      >  1:4171744/1‑76 (MQ=255)
                                  acaacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGtt                                      <  1:4414753/76‑1 (MQ=255)
                                   caacGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTg                                     <  2:5217055/76‑1 (MQ=255)
                                      cGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGatt                                  <  1:4186879/76‑1 (MQ=255)
                                       gggTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGatta                                 <  1:132050/76‑1 (MQ=255)
                                        ggTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGattat                                >  1:115025/1‑76 (MQ=255)
                                        ggTGTATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGattat                                <  2:5529878/76‑1 (MQ=255)
                                         gtgtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGattatt                               <  1:5066127/76‑1 (MQ=255)
                                         gtgtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGattatt                               <  1:3613142/76‑1 (MQ=255)
                                          tgtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGAttatta                              >  2:2184511/1‑76 (MQ=255)
                                          tgtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGAttatta                              >  1:4420612/1‑76 (MQ=255)
                                          tgtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGAttatta                              >  1:2578482/1‑76 (MQ=255)
                                           gtATTTTAGTGATAAAAACAAGAAGGCTAGTTATAAAATAATTAATTATGGAGACTATTATGATGTTGATtattat                             <  1:2071368/76‑1 (MQ=255)
                                           gtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATtattat                             <  2:4891417/76‑1 (MQ=255)
                                           gtATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATtattat                             >  2:4265018/1‑76 (MQ=255)
                                            tATTTTAGTGATAAAAACAAGAAGGCTAGTTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATg                            >  2:2308077/1‑76 (MQ=255)
                                            tATTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATg                            <  1:6072449/76‑1 (MQ=255)
                                             aTTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGa                           <  1:6037789/76‑1 (MQ=255)
                                             aTTTTAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGa                           >  1:1485902/1‑76 (MQ=255)
                                              ttttAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGAt                          <  1:4817053/76‑1 (MQ=255)
                                              ttttAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGAt                          <  2:1510575/76‑1 (MQ=255)
                                                ttAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATaa                        >  2:1329472/1‑76 (MQ=255)
                                                ttAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATaa                        <  2:5235679/76‑1 (MQ=255)
                                                 tAGTGATAAAAACAAGAAGGCTAGTTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAAc                       <  1:5442127/76‑1 (MQ=255)
                                                 tAGTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAAc                       >  2:2331414/1‑76 (MQ=255)
                                                   gTGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACaa                     <  2:2182212/76‑1 (MQ=255)
                                                    tGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAAt                    >  1:1236945/1‑76 (MQ=255)
                                                    tGATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAAt                    >  2:686038/1‑76 (MQ=255)
                                                     gATAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAAtt                   >  2:2515280/1‑76 (MQ=255)
                                                       tAAAAACAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTg                 <  2:1052981/76‑1 (MQ=255)
                                                        aaaaaCAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGa                >  1:720645/1‑76 (MQ=255)
                                                         aaaaCAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGaa               >  2:5912851/1‑76 (MQ=255)
                                                         aaaaCAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGaa               >  1:1999890/1‑76 (MQ=255)
                                                           aaCAAGAAGGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGaaaa             >  2:1501489/1‑76 (MQ=255)
                                                              aagaagGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATg          >  1:1111007/1‑76 (MQ=255)
                                                                gaagGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGaa        <  2:1197878/76‑1 (MQ=255)
                                                                gaagGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGaa        >  2:3876587/1‑76 (MQ=255)
                                                                 aagGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAg       >  1:3446924/1‑76 (MQ=255)
                                                                  agGCTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGt      >  1:3928325/1‑76 (MQ=255)
                                                                     cTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGtttt   >  2:3865641/1‑76 (MQ=255)
                                                                     cTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGtttt   >  1:5385005/1‑76 (MQ=255)
                                                                     cTAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGtttt   <  2:1364552/76‑1 (MQ=255)
                                                                      tAGCTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGttttt  >  1:5360686/1‑76 (MQ=255)
                                                                         |                                                                        
GTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGGGTGTATTTTAGTGATAAAAACAAGAAGGCTAGTTATAAAATAATTAATTATGGAGACTATTATGATGTTGATTATTATGATAACAATTTGAAAAATGAAGTTTTT  >  REL606/1504797‑1504942

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