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,2040GA88.9% 175.4 / 7.9 54intergenic (‑34/+635)yobF/yebOhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (25/23);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCG  >  REL606/1886131‑1886272
                                                                         |                                                                    
cgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTAtt                                                                    >  2:494730/1‑76 (MQ=255)
 gtgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTa                                                                   >  1:951588/1‑76 (MQ=255)
  tgtTTACTCAGGANTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCtattat                                                                  >  1:854508/1‑76 (MQ=255)
   gtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATc                                                                 >  2:62271/1‑76 (MQ=255)
      tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcaca                                                              >  1:344017/1‑76 (MQ=255)
      tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcaca                                                              >  1:290209/1‑76 (MQ=255)
       aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACAt                                                             >  1:824698/1‑76 (MQ=255)
          cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACa                                                          <  1:918379/76‑1 (MQ=255)
          cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACa                                                          <  2:961770/76‑1 (MQ=255)
            ggACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAAc                                                        <  2:848035/76‑1 (MQ=255)
            ggACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAAc                                                        <  2:70974/76‑1 (MQ=255)
               cttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGtt                                                     <  1:389511/76‑1 (MQ=255)
               cttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGtt                                                     <  1:674427/76‑1 (MQ=255)
                 tcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTaa                                                   <  1:804104/76‑1 (MQ=255)
                  cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAt                                                  <  2:205490/76‑1 (MQ=255)
                  cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAt                                                  >  2:412645/1‑76 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACTTCGTTTAACCCGAAATCTATTATCACATACAACGTTAAtt                                                 >  2:1141209/1‑76 (MQ=255)
                   tttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAAtt                                                 <  2:165176/76‑1 (MQ=255)
                    ttACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTa                                                >  1:175688/1‑76 (MQ=255)
                     tACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATCCAACGCTCATTAt                                               >  1:151128/1‑76 (MQ=255)
                     tACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTAt                                               >  1:1050029/1‑76 (MQ=255)
                     tACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTAt                                               <  2:350017/76‑1 (MQ=255)
                     tACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTAt                                               <  2:1066702/76‑1 (MQ=255)
                      aCTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATg                                              >  1:638244/1‑76 (MQ=255)
                      aCTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATg                                              >  1:1191176/1‑76 (MQ=255)
                      aCTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATg                                              <  1:417789/76‑1 (MQ=255)
                       cTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATgg                                             <  1:1004347/76‑1 (MQ=255)
                        tGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGc                                            <  2:42706/76‑1 (MQ=255)
                         gAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCa                                           >  2:667937/1‑76 (MQ=255)
                          aaaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAATGTTAATTATGGCaa                                          <  1:678789/76‑1 (MQ=255)
                          aaaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCaa                                          >  1:360391/1‑76 (MQ=255)
                           aaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAg                                         <  1:54356/76‑1 (MQ=255)
                               gCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAt                                     >  1:149802/1‑76 (MQ=255)
                                   cacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAttttt                                 >  2:151453/1‑76 (MQ=255)
                                   cacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAttttt                                 >  2:402098/1‑76 (MQ=255)
                                     caTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa                               <  1:1126884/76‑1 (MQ=255)
                                      aTAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTaaa                              >  2:53017/1‑76 (MQ=255)
                                           cAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGt                         <  2:32883/76‑1 (MQ=255)
                                             gAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtct                       <  2:560310/76‑1 (MQ=255)
                                               aCTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCg                     >  2:834724/1‑76 (MQ=255)
                                               aCTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCg                     <  1:192771/76‑1 (MQ=255)
                                                cTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGa                    >  2:281182/1‑76 (MQ=255)
                                                 tGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGAt                   >  1:539678/1‑76 (MQ=255)
                                                  gTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATc                  >  1:592566/1‑76 (MQ=255)
                                                   tACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa                 <  2:635781/76‑1 (MQ=255)
                                                   tACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa                 <  1:310870/76‑1 (MQ=255)
                                                    aCCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAg                >  1:1088362/1‑76 (MQ=255)
                                                      cTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGAc              >  2:1056741/1‑76 (MQ=255)
                                                       ccGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACg             <  2:394845/75‑1 (MQ=255)
                                                        cGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGc            >  1:995875/1‑76 (MQ=255)
                                                             aaCCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaaa       >  1:1181976/1‑76 (MQ=255)
                                                               cccGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtc     <  1:1105474/76‑1 (MQ=255)
                                                               cccGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtc     <  1:802364/76‑1 (MQ=255)
                                                                 cGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtctc   >  2:1170635/1‑76 (MQ=255)
                                                                  gAAATCTATTATCACATACAACGTTAATTATGGCATGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCg  >  2:1133784/1‑76 (MQ=255)
                                                                         |                                                                    
CGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCG  >  REL606/1886131‑1886272

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: