Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 430,835 C→T 100% intergenic (‑48/‑108) insL‑2 ← / → lon putative transposase insL for insertion sequence IS186/DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606430,8350CT92.0% 66.6 / ‑0.9 25intergenic (‑48/‑108)insL‑2/lonputative transposase insL for insertion sequence IS186/DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (13/10);  total (14/11)
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.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGCCGAAGTATCCAGTTCTTCGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCC  >  REL606/430712‑430978
                                                                                                                           |                                                                                                                                               
gTGCCGAAGTATCCAGTTCTTCGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATAc                                                                                                                       <  2:479958/150‑1 (MQ=255)
     gAAGTATCCAGTTCTTCGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACg                                                                                                                  <  1:553545/150‑1 (MQ=255)
          aTCCAGTTCTTCGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACAt                                                                                                             >  2:259884/1‑150 (MQ=255)
          aTCCAGTTCTTCGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACAt                                                                                                             >  2:681586/1‑150 (MQ=255)
                ttcttcGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAAt                                                                                                       <  2:568687/150‑1 (MQ=255)
                gtcttcGTGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACGGACGTACATGTTAAt                                                                                                       <  1:679588/149‑1 (MQ=255)
                  cttcTGGCATACCCATATGGGCCAGAATTGCTGGCAAGATATCGTGAGAGTAATTGATCGGCACGTTACATCATATCAGGCGCAAGACCACAACCCTTAAGTTAGTGCTTATGGGGAAACGTCCCCATATACTGACGTACATGTTGATAg                                                                                                     <  1:137/150‑1 (MQ=255)
                           tACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                                                            >  2:590525/1‑150 (MQ=255)
                           tACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                                                            >  2:672561/1‑150 (MQ=255)
                           tACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                                                            >  2:514070/1‑150 (MQ=255)
                           tACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                                                            >  2:380774/1‑150 (MQ=255)
                              cAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGc                                                                                         <  1:379670/150‑1 (MQ=255)
                                        cAGAATTGCTGACGAGTTAGCGTGAGAGTATTTCATCGGAACGTTAAAACATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACAGCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt                                                                               <  1:609850/150‑1 (MQ=21)
                                          gAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCa                                                                             <  1:259143/150‑1 (MQ=255)
                                          gAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCAAGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCa                                                                             <  1:588789/150‑1 (MQ=255)
                                                   gccAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                                                                    <  2:546648/149‑1 (MQ=255)
                                                    ccAGTTATCGTGAGAGTAATTCATCGGCACGTTTAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATTGCGTGAAGCACAGTCGTGTCATATGATTAcc                                                                   >  1:545068/1‑150 (MQ=255)
                                                                                   ttAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTAt                                    <  1:670581/150‑1 (MQ=255)
                                                                                                ggCGTAATACCACAACCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGt                       >  2:565231/1‑150 (MQ=255)
                                                                                                       caCCACAGCCCGTCAGTTAGGGCTGCGGGTGCAACATCCTCATGTACTGACGTAGATGAGAATAGATGGCGTGAAGCACGGTCGTGTCATCTGCGTACCTGGCGGAAATTAAGCTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAAc                <  1:172603/149‑1 (MQ=255)
                                                                                                              aCCCTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCCATGAATCCTGAGCGTTCTGAACGGATTGa         >  1:243392/1‑150 (MQ=255)
                                                                                                               cccTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGaa        >  2:563776/1‑150 (MQ=255)
                                                                                                               cccTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGaa        >  2:178005/1‑150 (MQ=255)
                                                                                                                ccTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGaaa       >  2:274470/1‑150 (MQ=255)
                                                                                                                ccTTAAGTTAGTGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGaaa       >  2:282064/1‑150 (MQ=255)
                                                                                                                     aGTTAGTGCTTATGGGGAAACATCCGCATATACTGACGTACATGTTAATAGATGGCGTGCAGCACAGTCGTGTCATCTGATTACCTTGCGGAAATTAGACTCAGATAGAGCTTCAAAAATCCTGAGGGATCTTAAAGCGTTTAAATccca  >  1:682105/1‑149 (MQ=255)
                                                                                                                           |                                                                                                                                               
GTGCCGAAGTATCCAGTTCTTCGGGCTTACCAATATGGGCCAGAATTGCTGACCAGTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCC  >  REL606/430712‑430978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: