Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,160,666 T→G intergenic (+130/‑190) rplL → / → rpoB 50S ribosomal protein L7/L12/DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,160,6660TG93.3% 173.7 / ‑1.9 45intergenic (+130/‑190)rplL/rpoB50S ribosomal protein L7/L12/DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (22/20);  total (24/21)
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.97e-01

ACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAG  >  REL606/4160595‑4160738
                                                                       |                                                                        
cncNGCCTTTNTGCGCTGTAAGGCGCCAGTAGCGTATCACACTNTATGACTACTGCTATNCCTTTCACTGCGTGtt                                                                      <  1:112372/74‑1 (MQ=255)
ccgaGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGtt                                                                      <  2:885255/73‑1 (MQ=255)
  cAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTc                                                                    <  2:1367839/76‑1 (MQ=255)
   gcccTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCt                                                                   <  1:710439/74‑1 (MQ=255)
   aGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATTCGTGTTTCt                                                                   >  1:43896/1‑76 (MQ=255)
      cTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATc                                                                >  2:900763/1‑76 (MQ=255)
      cTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATc                                                                >  2:501979/1‑76 (MQ=255)
       tttttGCGCTGTAAGGCACCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATcg                                                               >  2:481154/1‑76 (MQ=255)
         tttGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATcgac                                                             >  1:561314/1‑76 (MQ=255)
          ttGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATcgacg                                                            >  2:61708/1‑76 (MQ=255)
           tGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATcgacga                                                           >  1:611200/1‑76 (MQ=255)
                 gTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAat                                                     <  1:1224871/76‑1 (MQ=255)
                  tAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAata                                                    >  1:385089/1‑76 (MQ=255)
                  tAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAata                                                    >  2:371205/1‑76 (MQ=255)
                    aGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAAtata                                                  >  2:1333873/1‑76 (MQ=255)
                        gcNAGTCGCGTCTCACACTGTTTGAGTACTGCTGTGCCTTTCAATGCGTGTTTCTATCTACGACTTAATATTCTGc                                              <  2:1442482/76‑1 (MQ=255)
                        gcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGc                                              <  1:308969/76‑1 (MQ=255)
                          cAGTAGCGTTTCACACTGTTTGACTACTGCAGTGCCTTTCAATGCGTGTTTCTATCGACGACNTAATATACTGCGa                                            >  2:1401539/1‑76 (MQ=255)
                           aGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGAc                                           >  1:615870/1‑76 (MQ=255)
                             tAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAg                                         >  1:1273279/1‑76 (MQ=255)
                                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATNCGTGTTTACANCGCCAACTTAATATACTGCGCCANGACg                                     >  1:175403/1‑76 (MQ=255)
                                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACg                                     >  1:1040401/1‑76 (MQ=255)
                                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACg                                     >  1:1040391/1‑76 (MQ=255)
                                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGt                                    >  2:144819/1‑76 (MQ=255)
                                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGt                                    >  1:1005276/1‑76 (MQ=255)
                                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGt                                    >  2:1091698/1‑76 (MQ=255)
                                    tCACACTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATAAACTGCGACAGGACGTcc                                  >  2:1075065/1‑76 (MQ=255)
                                     cacaCTGTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTNTCGACGACTTAATATACTGCGACAGGACGTCCg                                 <  2:1295857/76‑1 (MQ=255)
                                     cacaCTGTTTGACTACTGCTGTGCCTTTCAACGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCg                                 <  2:245317/76‑1 (MQ=255)
                                           gTTTGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACNACTTAATATACTGCGACAGGACGTCCGTTCtgt                           <  2:1304262/76‑1 (MQ=255)
                                             ttGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTgtgt                         >  2:1337310/1‑76 (MQ=255)
                                             ttGACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTgtgt                         >  1:382725/1‑76 (MQ=255)
                                               gACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTaa                       <  2:936050/76‑1 (MQ=255)
                                               gACTACTGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTaa                       <  2:478229/76‑1 (MQ=255)
                                                actactGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTaaa                      >  2:1510790/1‑76 (MQ=255)
                                                 ctactGCTGTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAAt                     <  1:402887/76‑1 (MQ=255)
                                                    ctgctgTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATAGACTGCGACAGGACGTCCGTTCTGGGTAAATCGc                  <  2:1313829/76‑1 (MQ=255)
                                                    ctgctgTGCCTTTCAATGCGTGTTTCTATAGACGACATAATATACTGCGACAGGACGNCCGTTCTGTGTAAATCGc                  <  1:1454175/76‑1 (MQ=255)
                                                     tgctgTGCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCa                 <  1:501980/76‑1 (MQ=255)
                                                         gtgCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGa             >  1:1496580/1‑76 (MQ=255)
                                                          tgCCTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGaa            <  1:1194136/76‑1 (MQ=255)
                                                            ccTTTCAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAAt          <  2:814386/76‑1 (MQ=255)
                                                                 cAATGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGttt     <  1:1091699/76‑1 (MQ=255)
                                                                  aaTGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTa    <  1:563604/76‑1 (MQ=255)
                                                                    tGCGTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAg  <  2:797608/76‑1 (MQ=255)
                                                                       |                                                                        
ACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAG  >  REL606/4160595‑4160738

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: