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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,886,6080GA97.3% 1318.3 / 13.5 336intergenic (‑438/+231)yobF/yebOhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (138/189);  total (142/194)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 339 total aligned reads displayed.

AGTCAATCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACGTTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGAAAAT  >  REL606/1886535‑1886680
                                                                         |                                                                        
 gTCAATCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACAttt                                                                       >  1:5467099/1‑76 (MQ=255)
   cAATCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCg                                                                     >  1:4022304/1‑76 (MQ=255)
   cAATCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCg                                                                     >  1:4227758/1‑76 (MQ=255)
   cAATCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCg                                                                     >  2:6644107/1‑76 (MQ=255)
     aTCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCt                                                                   <  1:34576/76‑1 (MQ=255)
     aTCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCt                                                                   >  2:3198591/1‑76 (MQ=255)
     aTCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTCCGCt                                                                   <  1:8198490/76‑1 (MQ=255)
      tCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACGTTTCGCta                                                                  >  2:8396638/1‑76 (MQ=255)
       ccccAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCtaa                                                                 >  1:4436334/1‑76 (MQ=255)
       ccccAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCtaa                                                                 >  1:1848615/1‑76 (MQ=255)
       ccccAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCtaa                                                                 <  2:8036400/76‑1 (MQ=255)
        cccAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCtaat                                                                >  2:3345682/1‑76 (MQ=255)
           aaaGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTaataaa                                                             <  2:4451156/76‑1 (MQ=255)
              gggATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAataaataa                                                          >  1:3673377/1‑76 (MQ=255)
                 aTTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATAc                                                       <  2:2444402/76‑1 (MQ=255)
                  ttGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATAcc                                                      >  1:5040020/1‑76 (MQ=255)
                  ttGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATAcc                                                      <  1:7328740/76‑1 (MQ=255)
                   tGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCa                                                     <  1:8044531/76‑1 (MQ=255)
                   tGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCa                                                     <  2:5011090/76‑1 (MQ=255)
                   tGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCa                                                     >  1:4137427/1‑76 (MQ=255)
                     aaaaTTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACg                                                   <  2:9232322/76‑1 (MQ=255)
                            aaaaTAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATc                                            <  1:5926699/76‑1 (MQ=255)
                            aaaaTAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATc                                            >  1:5857720/1‑76 (MQ=255)
                            aaaaTAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATc                                            >  1:7237515/1‑76 (MQ=255)
                                 aaGAAAAATTGATGAATGAGCAAAAAATCAATAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACaaa                                       >  1:1273685/1‑76 (MQ=255)
                                  aGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAAt                                      >  1:3291255/1‑76 (MQ=255)
                                    aaaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATcc                                    <  1:9387472/76‑1 (MQ=255)
                                    aaaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATcc                                    <  1:895827/76‑1 (MQ=255)
                                    aaaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATcc                                    <  1:6557958/76‑1 (MQ=255)
                                     aaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCg                                   >  2:6439453/1‑76 (MQ=255)
                                     aaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCg                                   >  2:765343/1‑76 (MQ=255)
                                     aaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCg                                   >  2:1236415/1‑76 (MQ=255)
                                     aaaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCg                                   >  1:8417383/1‑76 (MQ=255)
                                      aaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGa                                  >  2:5627987/1‑76 (MQ=255)
                                      aaaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGa                                  >  1:2207229/1‑76 (MQ=255)
                                       aaTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGAt                                 >  2:1994450/1‑76 (MQ=255)
                                        aTTGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATg                                <  1:1923212/76‑1 (MQ=255)
                                         ttGATGAATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGt                               <  1:3861199/76‑1 (MQ=255)
                                          tgatgaATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTa                              >  1:1248751/1‑76 (MQ=255)
                                           gatgaATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTaa                             <  2:2198243/76‑1 (MQ=255)
                                           gatgaATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTaa                             >  1:7682923/1‑76 (MQ=255)
                                           gatgaATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTaa                             <  2:2645490/76‑1 (MQ=255)
                                           gatgaATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTaa                             <  1:3170240/76‑1 (MQ=255)
                                           gatgaATGAGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTaa                             >  1:4889547/1‑76 (MQ=255)
                                            atgaatgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTaaa                            <  2:8626231/76‑1 (MQ=255)
                                             tgaatgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAAc                           <  2:8229416/76‑1 (MQ=255)
                                             tgaatgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAAc                           <  1:2540333/76‑1 (MQ=255)
                                              gaatgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACa                          <  2:3463077/76‑1 (MQ=255)
                                              gaatgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACa                          <  2:1073181/76‑1 (MQ=255)
                                               aatgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACAt                         >  1:2319304/1‑76 (MQ=255)
                                                atgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATc                        >  1:2355647/1‑76 (MQ=255)
                                                 tgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATcc                       >  2:737487/1‑76 (MQ=255)
                                                 tgaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATcc                       >  2:1492042/1‑76 (MQ=255)
                                                  gaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCt                      <  1:1412450/76‑1 (MQ=255)
                                                  gaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCt                      >  2:6052134/1‑76 (MQ=255)
                                                  gaGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCt                      <  2:1828640/76‑1 (MQ=255)
                                                   aGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCta                     >  2:2934939/1‑76 (MQ=255)
                                                   aGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCta                     >  2:9419016/1‑76 (MQ=255)
                                                   aGCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCta                     <  1:5636340/76‑1 (MQ=255)
                                                    gCAAAAAATCAAGAGAGAAACGTTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaa                    <  1:7575658/76‑1 (MQ=255)
                                                    gCAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaa                    <  2:1775850/76‑1 (MQ=255)
                                                     cAAAAAATCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaat                   <  1:5244266/76‑1 (MQ=255)
                                                      aaaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaata                  <  1:7088485/76‑1 (MQ=255)
                                                      aaaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaata                  >  2:31560/1‑76 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:6824305/76‑1 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:1580328/76‑1 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 >  2:6927453/1‑76 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:2116043/76‑1 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:6124607/76‑1 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:466205/76‑1 (MQ=255)
                                                       aaaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:5320357/76‑1 (MQ=255)
                                                       aaaaaCCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCtaataa                 <  1:698793/76‑1 (MQ=255)
                                                        aaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATaaa                >  1:9102335/1‑76 (MQ=255)
                                                        aaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATaaa                <  2:8512875/76‑1 (MQ=255)
                                                        aaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATaaa                >  1:7276886/1‑76 (MQ=255)
                                                        aaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATaaa                >  1:5067677/1‑76 (MQ=255)
                                                        aaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATaaa                <  2:6312359/76‑1 (MQ=255)
                                                        aaaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATaaa                <  2:5462450/76‑1 (MQ=255)
                                                         aaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAAt               <  1:6420954/76‑1 (MQ=255)
                                                         aaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAAt               >  2:3690412/1‑76 (MQ=255)
                                                         aaaTCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAAt               >  1:6769620/1‑76 (MQ=255)
                                                           aTCAAGAGAGAAACGTTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTa             <  1:876498/76‑1 (MQ=255)
                                                            tCAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTaa            >  2:4231406/1‑76 (MQ=255)
                                                             cAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAAt           <  1:6288323/76‑1 (MQ=255)
                                                             cAAGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAAt           <  1:204599/76‑1 (MQ=255)
                                                              aaGAGAGAAACATTTCTCTAATAAATAATACCACGCTTAATCACAAATCCGATGGAAACATCCTAATAAATTAATg          <  2:4436333/76‑1 (MQ=255)
                                                              aaGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATg          >  2:4844793/1‑76 (MQ=255)
                                                              aaGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATg          <  1:785615/76‑1 (MQ=255)
                                                              aaGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATg          <  1:1033015/76‑1 (MQ=255)
                                                              aaGAGAGAAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATg          >  2:7216700/1‑76 (MQ=255)
                                                                gagagaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAg        <  1:8009029/76‑1 (MQ=255)
                                                                  gagaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGa      <  1:2918129/76‑1 (MQ=255)
                                                                  gagaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGa      <  1:9499821/76‑1 (MQ=255)
                                                                   agaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGaa     <  1:2974188/76‑1 (MQ=255)
                                                                    gaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGaaa    <  2:52719/76‑1 (MQ=255)
                                                                    gaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGaaa    >  1:957446/1‑76 (MQ=255)
                                                                    gaAACATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGaaa    >  1:5733984/1‑76 (MQ=255)
                                                                     aaaCATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGaaaa   <  1:2592664/76‑1 (MQ=255)
                                                                      aaCGTTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGAAAat  <  1:1695097/76‑1 (MQ=255)
                                                                      aaCATTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGAAAat  >  1:4380767/1‑76 (MQ=255)
                                                                         |                                                                        
AGTCAATCCCCAAAGGGATTGAAAATTTAAAATAAGAAAAATTGATGAATGAGCAAAAAATCAAGAGAGAAACGTTTCGCTAATAAATAATACCACGCTTAATCACAAATCCGATGTAAACATCCTAATAAATTAATGAGGAAAAT  >  REL606/1886535‑1886680

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