Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,036,374 (A)9→8 coding (585/1113 nt) wbbC ← putative glycosyltransferase WbbC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,036,3660A.96.4% 256.1 / 3.5 55coding (593/1113 nt)wbbCputative glycosyltransferase WbbC
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (29/24);  total (30/25)
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.52e-01

CTAATTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTACAATAA  >  REL606/2036292‑2036441
                                                                          |                                                                           
cTAATTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaa                                                                            >  2:796490/1‑76 (MQ=255)
 tAATTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaa                                                                           <  1:1055720/76‑1 (MQ=255)
 tAATTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaa                                                                           >  2:1414481/1‑76 (MQ=255)
  aaTTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaa                                                                          >  1:48185/1‑76 (MQ=255)
  aaTTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaa                                                                          >  1:1552793/1‑76 (MQ=255)
    ttttttGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaaan                                                                        <  1:30568/76‑2 (MQ=255)
      ttttGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaaaaaa                                                                      <  1:1167162/76‑1 (MQ=255)
       tttGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCACNGCACCCACCGAGaaaaaaaag                                                                     >  2:289823/1‑75 (MQ=255)
        ttGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaaaaaagg                                                                    >  1:1765982/1‑74 (MQ=255)
         tGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaaaaaaggt                                                                   >  2:1078439/1‑73 (MQ=255)
         tGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGaaaaaaaaggt                                                                   >  2:1065994/1‑73 (MQ=255)
            aCTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTca                                                               <  2:1061142/76‑1 (MQ=255)
              tttAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTcaca                                                             <  2:363743/76‑1 (MQ=255)
                tAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAAt                                                           >  2:50475/1‑76 (MQ=255)
                 aGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATc                                                          <  1:1315938/76‑1 (MQ=255)
                 aGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATc                                                          >  1:960708/1‑76 (MQ=255)
                   aaaaGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCat                                                        >  2:591243/1‑76 (MQ=255)
                   aaaaGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCat                                                        <  1:118327/76‑1 (MQ=255)
                   aaaaGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCat                                                        >  2:885258/1‑76 (MQ=255)
                    aaaGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCata                                                       >  1:1000464/1‑76 (MQ=255)
                            tATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGCAAAAAAAGGTTCACAATCATATTGGTGta                                               >  2:318006/1‑76 (MQ=255)
                            tATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGCAAAAAAAGGTTCACAATCATATTGGTGta                                               >  2:317997/1‑76 (MQ=255)
                            tATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGta                                               >  2:729847/1‑76 (MQ=255)
                             atcatcAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGtat                                              <  2:1681429/76‑1 (MQ=255)
                              tcatcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGtata                                             >  2:1277825/1‑76 (MQ=255)
                              tcatcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGtata                                             <  1:643449/76‑1 (MQ=255)
                              tcatcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGtata                                             <  1:643454/76‑1 (MQ=255)
                               catcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATag                                            <  1:474803/76‑1 (MQ=255)
                               catcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATag                                            <  2:220189/76‑1 (MQ=255)
                               catcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATag                                            <  2:1557318/76‑1 (MQ=255)
                                atcaAATTGGTTTGGTTTTTTTCATTGAAAGGCCCCCCCGGAAAAAAAAAAGGTTCACAATCATATTGGTGTATag                                            <  2:220173/76‑1 (MQ=255)
                                 tcaAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATagag                                          <  2:1719230/76‑1 (MQ=255)
                                  caAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAgg                                         >  1:482727/1‑76 (MQ=255)
                                     aTCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAtt                                      <  2:863514/76‑1 (MQ=255)
                                       cGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCGCAATCATATTGGTGTATAGAGGAtttt                                    >  1:1247036/1‑76 (MQ=255)
                                       cGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCGCAATCATATTGGTGTATAGAGGAtttt                                    >  1:1247047/1‑76 (MQ=255)
                                       cGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAtttt                                    <  2:389613/76‑1 (MQ=255)
                                        gTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAttttt                                   <  2:1773854/76‑1 (MQ=255)
                                         tATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAtttttt                                  >  1:623861/1‑76 (MQ=255)
                                         tATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAtttttt                                  <  1:487538/76‑1 (MQ=255)
                                         tATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAtttttt                                  <  1:308050/76‑1 (MQ=255)
                                         tATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGAtttttt                                  >  1:1188595/1‑76 (MQ=255)
                                          aTTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTa                                 >  1:670320/1‑76 (MQ=255)
                                          aTTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTa                                 >  1:249903/1‑76 (MQ=255)
                                           ttGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTaa                                <  2:28749/76‑1 (MQ=255)
                                           ttGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTaa                                >  2:436790/1‑76 (MQ=255)
                                           ttGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTaa                                >  2:1523864/1‑76 (MQ=255)
                                               ttATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATc                            >  2:412999/1‑76 (MQ=255)
                                               ttATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATc                            <  1:1432517/76‑1 (MQ=255)
                                                tATCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCa                           <  2:821325/76‑1 (MQ=255)
                                                 aTCTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCa                           >  2:1441883/1‑76 (MQ=255)
                                                  tCTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCAcc                         >  1:1469931/1‑76 (MQ=255)
                                                   cTTCATGCAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCt                        <  1:973958/76‑1 (MQ=255)
                                                          cAAAGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGTTGTATAGAGGATTTTTTAAAATCACCTAatatat                 <  1:372661/76‑1 (MQ=255)
                                                           aaaGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATAtt                >  2:1378799/1‑76 (MQ=255)
                                                           aaaGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATAtt                >  2:912091/1‑76 (MQ=255)
                                                           aaaGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATAtt                <  1:1065098/76‑1 (MQ=255)
                                                            aaGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTg               >  2:947634/1‑76 (MQ=255)
                                                            aaGCACCCACCGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTg               <  2:147566/76‑1 (MQ=255)
                                                               cacccaccGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTaa            <  2:1056455/76‑1 (MQ=255)
                                                                acccaccGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAg           >  1:1745719/1‑76 (MQ=255)
                                                                   caccGAGAAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGtt         <  2:1575419/76‑1 (MQ=255)
                                                                   caccGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTa        >  2:1806106/1‑76 (MQ=255)
                                                                   caccGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTa        >  2:1726930/1‑76 (MQ=255)
                                                                   caccGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTa        >  2:117380/1‑76 (MQ=255)
                                                                    accGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTAc       >  2:1278535/1‑76 (MQ=255)
                                                                     ccGAGAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTACa      >  1:1574900/1‑76 (MQ=255)
                                                                          aaaaaaaaaGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTACAATaa  >  1:114537/1‑76 (MQ=255)
                                                                          |                                                                           
CTAATTTTTTGAACTTTAGAAAAGCCTTTATCATCAAATCGTATTGTTTATCTTCATGCAAAGCACCCACCGAGAAAAAAAAAGGTTCACAATCATATTGGTGTATAGAGGATTTTTTAAAATCACCTAATATATTGTAAGTTACAATAA  >  REL606/2036292‑2036441

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