Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,505,141 A→G intergenic (+175/‑7) ydcD → / → ECB_01416 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,505,1410AG98.7% 617.9 / ‑1.4 148intergenic (+175/‑7)ydcD/ECB_01416hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (71/75);  total (73/75)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 150 total aligned reads displayed.

ATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGTT  >  REL606/1505068‑1505213
                                                                         |                                                                        
aTTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCg                                                                        <  2:4576660/76‑1 (MQ=255)
aTTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCg                                                                        >  1:3953128/1‑76 (MQ=255)
   ttCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCACCa                                                                     <  1:2515280/76‑1 (MQ=255)
    tCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAg                                                                    >  1:75007/1‑76 (MQ=255)
     cATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAgg                                                                   >  1:1028516/1‑76 (MQ=255)
      aTCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGt                                                                  >  2:550515/1‑76 (MQ=255)
      aTCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGt                                                                  <  2:5370362/76‑1 (MQ=255)
      aTCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGt                                                                  <  1:1649925/76‑1 (MQ=255)
      aTCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGt                                                                  >  2:2619785/1‑76 (MQ=255)
       tCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTg                                                                 <  2:1485902/76‑1 (MQ=255)
         tttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAg                                                               >  1:5925786/1‑76 (MQ=255)
          ttGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAga                                                              >  1:5435471/1‑75 (MQ=255)
             cTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTAtt                                                           <  1:1213073/76‑1 (MQ=255)
             cTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTAtt                                                           >  2:4136711/1‑76 (MQ=255)
              tttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTc                                                          <  2:2387732/76‑1 (MQ=255)
              tttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTc                                                          >  2:3103868/1‑76 (MQ=255)
               ttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCa                                                         >  2:2719792/1‑76 (MQ=255)
               ttAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCa                                                         >  1:4534607/1‑76 (MQ=255)
                tAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCaa                                                        <  1:4690928/76‑1 (MQ=255)
                 aaCGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCaaa                                                       <  2:4420612/76‑1 (MQ=255)
                 aaCGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCaaa                                                       <  2:5435471/76‑1 (MQ=255)
                   cGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGt                                                     >  1:1224961/1‑76 (MQ=255)
                     ccTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGttt                                                   <  2:4058389/76‑1 (MQ=255)
                     acaaCAGGTGATCAAAATGTAACCGTTGGACTAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATACAAAGttt                                                   <  2:4058390/73‑1 (MQ=255)
                       tACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGc                                                 >  1:2655044/1‑76 (MQ=255)
                       tACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGc                                                 <  1:5370891/76‑1 (MQ=255)
                        aCAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCt                                                >  1:2217756/1‑76 (MQ=255)
                        aCAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCt                                                >  1:5448217/1‑76 (MQ=255)
                        aCAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCt                                                >  2:5604245/1‑76 (MQ=255)
                        aCAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCt                                                >  2:590874/1‑76 (MQ=255)
                         cAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCtt                                               >  2:1095330/1‑76 (MQ=255)
                          aGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTg                                              <  1:202883/76‑1 (MQ=255)
                          aGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTg                                              >  2:507024/1‑76 (MQ=255)
                          aGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTg                                              >  1:5238031/1‑76 (MQ=255)
                            gTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGAt                                            >  1:1244372/1‑76 (MQ=255)
                             tGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGAtt                                           >  1:1004638/1‑76 (MQ=255)
                             tGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGAtt                                           <  1:1089663/76‑1 (MQ=255)
                              gATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATta                                          <  1:4621954/76‑1 (MQ=255)
                               aTCAAAATTTCCCCTTTGGAATAACTAAAAAGATACAAACGGACGCCAGGCGAGACACCGAAGTTTGCTTGATtat                                         >  2:2634520/1‑76 (MQ=255)
                               aTCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATtat                                         <  2:6102972/76‑1 (MQ=255)
                               aTCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATtat                                         >  1:6116563/1‑76 (MQ=255)
                                  aaaaTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATAtt                                      <  2:1999890/76‑1 (MQ=255)
                                  aaaaTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATAtt                                      >  1:2644592/1‑76 (MQ=255)
                                   aaaTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATAttt                                     <  2:5360686/76‑1 (MQ=255)
                                    aaTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTc                                    <  1:3863192/76‑1 (MQ=255)
                                    aaTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTc                                    >  1:5090688/1‑76 (MQ=255)
                                     aTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCa                                   >  1:3527752/1‑76 (MQ=255)
                                     aTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCa                                   >  1:1674131/1‑76 (MQ=255)
                                     aTTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCa                                   <  2:4496979/76‑1 (MQ=255)
                                      tttAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAg                                  >  1:3169202/1‑76 (MQ=255)
                                       ttAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGt                                 <  2:4754330/76‑1 (MQ=255)
                                       ttAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGt                                 <  2:1111007/76‑1 (MQ=255)
                                        tAACCTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTg                                <  1:4050556/76‑1 (MQ=255)
                                           ccTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTGTATTTCAGTGAcc                             >  1:2949366/1‑76 (MQ=255)
                                            cTTTGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGAccc                            <  2:585563/76‑1 (MQ=255)
                                              ttGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCt                          >  1:1251932/1‑76 (MQ=255)
                                               tGGAATAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTg                         <  1:3672751/76‑1 (MQ=255)
                                                   aTAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGatat                     <  2:1032342/76‑1 (MQ=255)
                                                   aTAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGatat                     <  1:3251435/76‑1 (MQ=255)
                                                    tAACTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGAtata                    >  2:870312/1‑76 (MQ=255)
                                                     aaCTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATAc                   <  1:6236973/76‑1 (MQ=255)
                                                     aaCTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATAc                   >  2:2915531/1‑76 (MQ=255)
                                                      aCTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACg                  <  1:1501489/76‑1 (MQ=255)
                                                      aCTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACg                  <  1:4673977/76‑1 (MQ=255)
                                                      aCTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACg                  <  2:2540778/76‑1 (MQ=255)
                                                       cTAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGa                 <  1:69333/76‑1 (MQ=255)
                                                        tAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGac                >  2:5397058/1‑76 (MQ=255)
                                                        tAAAAAGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGac                >  1:2208787/1‑76 (MQ=255)
                                                        tAAAAAAAAAAAAAAAGGGGCCAGGGGAGTATTCAAAATTTGCTTGATTATATTTTAGTGACCCCCGATATACGaa                >  2:4917912/1‑75 (MQ=255)
                                                          aaaaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaa              >  2:979745/1‑76 (MQ=255)
                                                          aaaaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaa              >  1:246001/1‑76 (MQ=255)
                                                          aaaaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaa              <  2:1256973/76‑1 (MQ=255)
                                                            aaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca            <  1:1973800/76‑1 (MQ=255)
                                                            aaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca            >  2:4285672/1‑76 (MQ=255)
                                                            aaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca            <  1:3492802/76‑1 (MQ=255)
                                                            aaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca            >  2:549169/1‑76 (MQ=255)
                                                            aaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca            <  1:320496/76‑1 (MQ=255)
                                                            aaGATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGacaaca            <  2:3898399/76‑1 (MQ=255)
                                                              gATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAg          <  1:5121445/76‑1 (MQ=255)
                                                              gATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAg          <  1:686038/76‑1 (MQ=255)
                                                              gATAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAg          <  2:5595542/76‑1 (MQ=255)
                                                               aTAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAgg         <  1:1473514/76‑1 (MQ=255)
                                                               aTAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAgg         <  2:1058757/76‑1 (MQ=255)
                                                               aTAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAgg         <  2:3479833/76‑1 (MQ=255)
                                                               aTAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAgg         <  2:3167998/76‑1 (MQ=255)
                                                                tAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGa        >  2:781252/1‑76 (MQ=255)
                                                                tAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGa        >  2:2958065/1‑76 (MQ=255)
                                                                tAAAAAAGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGa        >  2:4623533/1‑76 (MQ=255)
                                                                 aaaaaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGaa       <  1:1195823/76‑1 (MQ=255)
                                                                 aaaaaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGaa       >  1:997687/1‑76 (MQ=255)
                                                                 aaaaaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGaa       <  2:4674632/76‑1 (MQ=255)
                                                                 aaaaaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGaa       <  2:1840615/76‑1 (MQ=255)
                                                                   aaaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAg     >  2:6072998/1‑76 (MQ=255)
                                                                   aaaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAg     >  2:951650/1‑76 (MQ=255)
                                                                    aaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAgg    <  2:2766406/76‑1 (MQ=255)
                                                                    aaaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAgg    <  2:194115/76‑1 (MQ=255)
                                                                    aaaGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAgg    >  2:4398861/1‑76 (MQ=255)
                                                                     aaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGt   >  2:5082721/1‑76 (MQ=255)
                                                                     aaGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGt   >  1:5909931/1‑76 (MQ=255)
                                                                      aGGGCGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGtt  <  2:1985135/76‑1 (MQ=255)
                                                                         |                                                                        
ATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGCTTGATTATATTTCAGTGACCCCTGATATACGACAACAAGGAAAGGTT  >  REL606/1505068‑1505213

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: