Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,886,204 G→A intergenic (‑34/+635) yobF ← / ← yebO hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,886,2040GA98.7% 892.8 / ‑1.6 222intergenic (‑34/+635)yobF/yebOhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (106/113);  total (109/113)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 233 total aligned reads displayed.

CGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATAC  >  REL606/1886131‑1886276
                                                                         |                                                                        
cgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTAtt                                                                        <  1:280883/76‑1 (MQ=255)
    tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATca                                                                    <  1:1365493/76‑1 (MQ=255)
    tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATca                                                                    <  2:4438506/76‑1 (MQ=255)
     ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcac                                                                   >  2:2521854/1‑76 (MQ=255)
     ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcac                                                                   >  2:1095277/1‑76 (MQ=255)
     ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcac                                                                   >  2:3114682/1‑76 (MQ=255)
      tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcaca                                                                  <  1:3941591/76‑1 (MQ=255)
      tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcaca                                                                  >  2:4041722/1‑76 (MQ=255)
       aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACAt                                                                 <  2:1875281/76‑1 (MQ=255)
       aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACAt                                                                 <  1:1431437/76‑1 (MQ=255)
       aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACAt                                                                 <  1:3235568/76‑1 (MQ=255)
        cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATa                                                                <  2:3391420/76‑1 (MQ=255)
        cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATa                                                                <  2:1967163/76‑1 (MQ=255)
         tCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATAc                                                               <  1:4942480/76‑1 (MQ=255)
         tCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATAc                                                               <  1:4815032/76‑1 (MQ=255)
         tCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATAc                                                               >  2:5853079/1‑76 (MQ=255)
          cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGACATCTATTATCACATACa                                                              >  1:897297/1‑76 (MQ=255)
          cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACa                                                              >  2:3469537/1‑76 (MQ=255)
           aGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACaa                                                             >  1:5365994/1‑76 (MQ=255)
               cttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGtt                                                         >  2:2650742/1‑76 (MQ=255)
               cttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGtt                                                         <  2:5018910/76‑1 (MQ=255)
               cttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGtt                                                         >  2:1585860/1‑76 (MQ=255)
                ttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTa                                                        <  1:5081350/76‑1 (MQ=255)
                ttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTa                                                        >  2:830153/1‑76 (MQ=255)
                 tcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTaa                                                       >  1:2366578/1‑76 (MQ=255)
                 tcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTaa                                                       <  1:1514776/76‑1 (MQ=255)
                  cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAt                                                      >  2:2427269/1‑76 (MQ=255)
                  cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAt                                                      >  1:6161717/1‑76 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACGTACAACGTTAAtt                                                     <  1:5796311/76‑1 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAtt                                                     <  1:44188/76‑1 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAtt                                                     <  2:2850983/76‑1 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAtt                                                     >  2:3726398/1‑76 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAtt                                                     >  1:6101444/1‑76 (MQ=255)
                    ttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTa                                                    >  1:3655111/1‑76 (MQ=255)
                    ttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTa                                                    >  1:3799604/1‑76 (MQ=255)
                    ttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTa                                                    >  1:2773074/1‑76 (MQ=255)
                        tGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGc                                                <  1:983969/76‑1 (MQ=255)
                        tGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGc                                                <  1:1102901/76‑1 (MQ=255)
                         gAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCa                                               >  2:5998668/1‑76 (MQ=255)
                         gAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCa                                               <  1:862901/76‑1 (MQ=255)
                           aaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAg                                             <  1:5226510/76‑1 (MQ=255)
                            aaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGc                                            <  2:2056234/76‑1 (MQ=255)
                            aaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGc                                            <  2:5081292/76‑1 (MQ=255)
                              tGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCa                                          <  1:1208879/76‑1 (MQ=255)
                              tGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCa                                          <  2:4551044/76‑1 (MQ=255)
                                ccACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAtt                                        >  2:5041642/1‑76 (MQ=255)
                                  acacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAtttt                                      <  2:1320039/76‑1 (MQ=255)
                                   cacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCAttttt                                     >  2:2969137/1‑76 (MQ=255)
                                   cacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAttttt                                     <  2:3496747/76‑1 (MQ=255)
                                    acaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTa                                    <  2:5343635/76‑1 (MQ=255)
                                    acaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTa                                    <  1:2525242/76‑1 (MQ=255)
                                    acaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTa                                    <  1:663580/76‑1 (MQ=255)
                                    acaTAAACAGAACGGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTa                                    <  2:6175498/76‑1 (MQ=255)
                                       tAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAAc                                 >  2:2303754/1‑76 (MQ=255)
                                       tAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAAc                                 >  2:3785945/1‑76 (MQ=255)
                                        aaaCAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa                                >  2:811128/1‑76 (MQ=255)
                                         aaCAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACAt                               >  1:3155808/1‑76 (MQ=255)
                                           cAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                             >  2:4210236/1‑76 (MQ=255)
                                           cAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                             >  2:5678483/1‑76 (MQ=255)
                                           cAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                             >  2:5284991/1‑76 (MQ=255)
                                             gAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtct                           <  2:651621/76‑1 (MQ=255)
                                              aaCTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtctc                          <  2:1624986/76‑1 (MQ=255)
                                               aCTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCg                         >  1:292944/1‑76 (MQ=255)
                                                cTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGa                        >  1:5267482/1‑76 (MQ=255)
                                                 tGTACCTCGTTTAACCCGAAATCTATTATCACATACaa                                                             >  1:3205941/1‑38 (MQ=255)
                                                 tGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGAt                       >  1:2232844/1‑76 (MQ=255)
                                                  gTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATc                      <  1:6255317/76‑1 (MQ=255)
                                                   tACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa                     >  2:6048428/1‑76 (MQ=255)
                                                   tACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa                     <  1:2197891/76‑1 (MQ=255)
                                                    aCCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAg                    <  2:5537797/76‑1 (MQ=255)
                                                    aCCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAg                    <  2:5711132/76‑1 (MQ=255)
                                                    aCCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAg                    <  1:4085544/76‑1 (MQ=255)
                                                     ccTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGa                   >  2:5454695/1‑76 (MQ=255)
                                                     ccTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGa                   <  1:3515276/76‑1 (MQ=255)
                                                     ccTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGa                   <  1:164447/76‑1 (MQ=255)
                                                      cTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGAc                  <  1:5927224/76‑1 (MQ=255)
                                                       tCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACg                 <  1:56697/76‑1 (MQ=255)
                                                       tCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACg                 <  2:3057407/76‑1 (MQ=255)
                                                       tCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACg                 >  2:2149928/1‑76 (MQ=255)
                                                        cGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGc                >  1:1833611/1‑76 (MQ=255)
                                                        cGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGc                <  1:2587448/76‑1 (MQ=255)
                                                         gTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGca               >  2:4688630/1‑76 (MQ=255)
                                                          tttAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGcac              <  1:3270147/76‑1 (MQ=255)
                                                           ttAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGcaca             >  1:128130/1‑76 (MQ=255)
                                                           ttAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGcaca             <  2:1469532/76‑1 (MQ=255)
                                                            tAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaa            <  2:5259234/76‑1 (MQ=255)
                                                             aaCCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaaa           <  2:5634378/76‑1 (MQ=255)
                                                             aaCCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaaa           >  1:4348114/1‑76 (MQ=255)
                                                              aCCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAt          >  2:3201490/1‑76 (MQ=255)
                                                                 cGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtctc       <  1:4710811/76‑1 (MQ=255)
                                                                  gAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCg      <  2:3852661/76‑1 (MQ=255)
                                                                   aaaTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGa     >  2:4888875/1‑76 (MQ=255)
                                                                   aaaTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGa     >  2:1953442/1‑76 (MQ=255)
                                                                    aaTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAt    <  1:2051753/76‑1 (MQ=255)
                                                                    aaTCTATTAGCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAt    <  1:5158171/76‑1 (MQ=255)
                                                                     aTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa   >  1:1897695/1‑76 (MQ=255)
                                                                     aTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa   <  2:2723074/76‑1 (MQ=255)
                                                                     aTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa   <  2:1192093/76‑1 (MQ=255)
                                                                     aTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa   >  1:4564663/1‑76 (MQ=255)
                                                                     aTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa   >  1:1013342/1‑76 (MQ=255)
                                                                         |                                                                        
CGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATAC  >  REL606/1886131‑1886276

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