Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,035,629 (C)5→6 coding (228/768 nt) wbbD ← putative glycosyltransferase WbbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,035,6241.C91.1% 217.6 / 13.6 56coding (233/768 nt)wbbDputative glycosyltransferase WbbD
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (28/23);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAATA  >  REL606/2035550‑2035692
                                                                           |                                                                    
ttCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGc                                                                     >  1:1737544/1‑76 (MQ=255)
ttCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGc                                                                     <  1:1536907/76‑1 (MQ=255)
 tCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGcc                                                                    >  2:485960/1‑76 (MQ=255)
    aTTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGccccc                                                                 >  2:687372/1‑76 (MQ=255)
    aTTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGccccc                                                                 >  1:801814/1‑76 (MQ=255)
    aTTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGccccc                                                                 >  2:1799269/1‑76 (MQ=255)
     tttATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGcccccc                                                                <  2:4159/76‑2 (MQ=255)
     tttATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGcccccc                                                                <  2:217125/76‑2 (MQ=255)
      ttATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGcccccca                                                               <  2:1040859/76‑3 (MQ=255)
      ttATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGcccccca                                                               >  1:1368123/1‑74 (MQ=255)
      ttATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGcccccca                                                               <  2:681579/76‑3 (MQ=255)
          atcatcAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAaccc                                                           <  1:1603664/76‑1 (MQ=255)
               caGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCatat                                                       >  1:1803884/1‑76 (MQ=255)
               caGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCatat                                                       >  1:1145812/1‑76 (MQ=255)
               caGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCatat                                                       >  1:602909/1‑76 (MQ=255)
                aGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATAtt                                                      >  2:847485/1‑76 (MQ=255)
                aGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATAtt                                                      >  1:510911/1‑76 (MQ=255)
                  aaTCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATAtttt                                                    >  1:1665062/1‑76 (MQ=255)
                       gCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCa                                               <  1:1038077/76‑1 (MQ=255)
                       gCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCa                                               <  2:492892/76‑1 (MQ=255)
                            tATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAgtgt                                          >  1:1818635/1‑76 (MQ=255)
                            tATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAgtgt                                          >  1:1881229/1‑76 (MQ=255)
                            tATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAgtgt                                          >  2:335678/1‑76 (MQ=255)
                             aTTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGtgtg                                         >  1:1410562/1‑76 (MQ=255)
                             aTTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGtgtg                                         >  1:1239991/1‑76 (MQ=255)
                             aTTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGtgtg                                         <  1:913835/76‑1 (MQ=255)
                               tACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGan                                       <  1:938253/76‑2 (MQ=255)
                               tACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGan                                       <  1:1385539/76‑2 (MQ=255)
                               tACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGaa                                       <  1:741643/76‑1 (MQ=255)
                                 cTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAAt                                     <  2:1235257/76‑1 (MQ=255)
                                  tAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATc                                    <  1:58237/76‑1 (MQ=255)
                                   aaaTCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCtt                                  <  2:440764/76‑1 (MQ=255)
                                      tCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCtttt                                <  1:363625/76‑1 (MQ=255)
                                       cATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTta                               >  1:649432/1‑76 (MQ=255)
                                       cATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTNta                               <  2:1333582/76‑1 (MQ=255)
                                        aTGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTtat                              <  1:750806/76‑1 (MQ=255)
                                        aTGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTtata                             >  1:1526854/1‑76 (MQ=255)
                                         tGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTtata                             >  2:787300/1‑76 (MQ=255)
                                         tGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTtata                             <  1:1713281/76‑1 (MQ=255)
                                         tGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTtata                             <  1:1665601/76‑1 (MQ=255)
                                          gCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTatat                            >  1:968744/1‑76 (MQ=255)
                                           cATACATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATAtt                           >  1:1709027/1‑76 (MQ=255)
                                              aCATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTaa                        <  2:1818635/76‑1 (MQ=255)
                                               cATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTaaa                       <  1:572351/76‑1 (MQ=255)
                                               cATTGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTaaa                       <  1:888879/76‑1 (MQ=255)
                                                  tGATTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGCGTGAAATCTTTTATATTTAAatat                    >  1:1202868/1‑76 (MQ=255)
                                                    aTTCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATAtct                  <  1:1817111/76‑1 (MQ=255)
                                                     ttCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATAtctt                 >  2:1297013/1‑76 (MQ=255)
                                                     ttCAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATAtctt                 <  2:1476726/76‑1 (MQ=255)
                                                       cAACCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATAtcttct               >  1:472493/1‑76 (MQ=255)
                                                        aaCCCAAAATTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTtctc              >  2:1497944/1‑76 (MQ=255)
                                                             aaaaTTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACt         <  1:715620/76‑1 (MQ=255)
                                                             aaaaTTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACt         <  1:1211926/76‑1 (MQ=255)
                                                              aaaTTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTa        >  1:387952/1‑76 (MQ=255)
                                                              aaaTTTAACGCTGCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTAt       >  2:1064778/1‑76 (MQ=255)
                                                               aaTTTAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTAt       <  1:805587/76‑1 (MQ=255)
                                                                 tttAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCa     >  2:1631112/1‑76 (MQ=255)
                                                                 tttAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCa     >  1:1075068/1‑76 (MQ=255)
                                                                  ttAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCaa    >  1:1747270/1‑76 (MQ=255)
                                                                  ttAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCaa    >  2:230617/1‑76 (MQ=255)
                                                                  ttAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCaa    >  1:195345/1‑76 (MQ=255)
                                                                   tAACGCTGCCCCCCAACCCCATATTTTTTTCACGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAAt   >  2:595991/1‑76 (MQ=255)
                                                                   tAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAAt   >  1:1866278/1‑76 (MQ=255)
                                                                   tAACGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAAt   <  1:94853/76‑1 (MQ=255)
                                                                   tAACGCTGCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAATa  <  2:1589369/76‑1 (MQ=255)
                                                                    aaCGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAATa  <  1:306904/76‑1 (MQ=255)
                                                                    aaCGCTGCCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAATa  >  1:1211232/1‑76 (MQ=255)
                                                                           |                                                                    
TTCTATTTATATCATCAGAATCAGCACGTATTACTAAATCATGCATACATTGATTCAACCCAAAATTTAACGCTGCCCCCAACCCCATATTTTTTTCAAGTGTGAAATCTTTTATATTTAAATATCTTCTCCAACTATCAATA  >  REL606/2035550‑2035692

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