Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,268,816 (A)8→7 intergenic (‑131/‑128) rplU ← / → ispB 50S ribosomal protein L21/octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,268,8090A.91.2% 226.6 / 19.7 57intergenic (‑124/‑135)rplU/ispB50S ribosomal protein L21/octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base . (18/34);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

ATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAA  >  REL606/3268734‑3268879
                                                                           |                                                                      
ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaaga                                                                        >  1:694197/1‑76 (MQ=255)
ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaaga                                                                        >  1:694215/1‑76 (MQ=255)
ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaaga                                                                        >  1:1683145/1‑76 (MQ=255)
ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaaga                                                                        >  2:1078052/1‑76 (MQ=255)
ataAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaaga                                                                        >  1:1353807/1‑76 (MQ=255)
 taAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACgaagaa                                                                       >  1:908496/1‑76 (MQ=255)
  aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCGCGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaa                                                                      <  1:1958711/76‑1 (MQ=255)
  aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaa                                                                      >  1:761056/1‑76 (MQ=255)
  aaaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaa                                                                      <  1:1117945/76‑1 (MQ=255)
   aaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaa                                                                     >  2:10550/1‑76 (MQ=255)
   aaTATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaa                                                                     >  2:154912/1‑76 (MQ=255)
    atatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaa                                                                    <  1:1887522/76‑1 (MQ=255)
    atatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaa                                                                    <  2:1724253/76‑1 (MQ=255)
     tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaaa                                                                   >  1:1471665/1‑76 (MQ=255)
     tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaaa                                                                   >  1:1471673/1‑76 (MQ=255)
     tatTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaaa                                                                   >  1:378434/1‑76 (MQ=255)
       ttCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGaaaaaaac                                                                 <  1:1234336/76‑2 (MQ=255)
            aaTAGGGCGCGAATATTACGCAAAACGCACGCCTTTAACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAAc                                                           >  2:20369/1‑76 (MQ=255)
              tAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACt                                                          <  1:1282837/76‑1 (MQ=255)
              tAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACtt                                                         <  2:273009/76‑1 (MQ=255)
              tAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACtt                                                         <  1:719492/76‑1 (MQ=255)
              tAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACtt                                                         >  1:1581320/1‑76 (MQ=255)
                 ggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTa                                                      <  2:1825225/76‑1 (MQ=255)
                 ggCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTa                                                      <  2:1556314/76‑1 (MQ=255)
                  gcgcGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTAc                                                     <  2:880659/76‑1 (MQ=255)
                  gcacgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTAc                                                     <  2:458894/73‑1 (MQ=255)
                   cgcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACg                                                    <  1:755040/76‑1 (MQ=255)
                    gcgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACg                                                    >  1:690743/1‑76 (MQ=255)
                     cgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGt                                                  <  2:500260/76‑1 (MQ=255)
                     cgAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGt                                                  <  2:1236737/76‑1 (MQ=255)
                        atatTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAGAAACACAACTTGTACGGTAAc                                               <  2:1296808/76‑1 (MQ=255)
                          atTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGt                                             >  1:994754/1‑76 (MQ=255)
                            tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGttt                                           >  2:1626662/1‑76 (MQ=255)
                            tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGttt                                           >  2:433090/1‑76 (MQ=255)
                            tACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGttt                                           <  2:840801/76‑1 (MQ=255)
                             aCGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTa                                          <  1:348651/76‑1 (MQ=255)
                             aCGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTa                                          <  1:1662808/76‑1 (MQ=255)
                              cGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTAt                                         >  1:1337478/1‑76 (MQ=255)
                              cGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTAt                                         <  2:949829/76‑1 (MQ=255)
                               gCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATc                                        <  1:816989/76‑1 (MQ=255)
                                cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCt                                       <  2:1567615/76‑1 (MQ=255)
                                cAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCt                                       <  2:1583322/76‑1 (MQ=255)
                                 aaaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCtg                                      >  2:477609/1‑76 (MQ=255)
                                 aaaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCtg                                      <  1:1424650/76‑1 (MQ=255)
                                  aaaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCtgt                                     >  2:1695102/1‑76 (MQ=255)
                                   aaCGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCtgtg                                    <  2:1054775/76‑1 (MQ=255)
                                   aaCGCACGCCTTTGACAAGTGCTACAGTAAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCtgtg                                    <  1:576802/76‑1 (MQ=255)
                                      gcacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCa                                 >  1:31767/1‑76 (MQ=255)
                                       cacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                >  2:999972/1‑76 (MQ=255)
                                       cacgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAt                                >  2:1395474/1‑76 (MQ=255)
                                        acgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                               <  1:69987/76‑1 (MQ=255)
                                        acgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                               <  2:1879698/76‑1 (MQ=255)
                                        acgcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtt                               <  1:1801847/76‑1 (MQ=255)
                                          gcCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAtttt                             >  1:157649/1‑76 (MQ=255)
                                           ccTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCAttttt                            <  1:148466/76‑1 (MQ=255)
                                                gACAAGTGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGt                       >  2:1734194/1‑76 (MQ=255)
                                                      tGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                 >  2:256008/1‑76 (MQ=255)
                                                      tGCTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATc                 <  1:1046885/76‑1 (MQ=255)
                                                        cTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAc               >  2:1719121/1‑76 (MQ=255)
                                                        cTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAc               >  1:1687960/1‑76 (MQ=255)
                                                        cTACAGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAc               >  1:1687952/1‑76 (MQ=255)
                                                          aCAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAcc              >  1:143350/1‑76 (MQ=255)
                                                          aCAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCAcc              >  2:871335/1‑76 (MQ=255)
                                                            aGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta           <  1:1826759/76‑1 (MQ=255)
                                                            aGTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCta           <  1:1467613/76‑1 (MQ=255)
                                                             gTCAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat          <  1:1209125/76‑1 (MQ=255)
                                                              tCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCtat          <  1:147393/76‑1 (MQ=255)
                                                               cAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTatat        >  2:1785080/1‑76 (MQ=255)
                                                               cAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTatat        <  2:1645084/76‑1 (MQ=255)
                                                               cAATACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTatat        <  1:1065145/76‑1 (MQ=255)
                                                                 aTACACGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTc      <  1:280704/76‑1 (MQ=255)
                                                                   acacGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt    <  1:1437875/76‑1 (MQ=255)
                                                                    cacGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCt    <  1:1283766/76‑1 (MQ=255)
                                                                    cacGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa   <  2:586694/76‑1 (MQ=255)
                                                                    cacGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTa   <  1:213941/76‑1 (MQ=255)
                                                                     acGAAGAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTaa  <  1:1216426/76‑1 (MQ=255)
                                                                           |                                                                      
ATAAATATTCACAATAGGGCGCGAATATTACGCAAAACGCACGCCTTTGACAAGTGCTACAGTCAATACACGAAGAAAAAAAACACAACTTGTACGGTAACGTTTATCTGTGCCATTTTTTCAGTACAATCACCCTATATTCCTAA  >  REL606/3268734‑3268879

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