Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,070,246:1 +GA 100% intergenic (+46/+96) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,070,2461.G93.1% 92.4 / 0.8 29intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (13/14);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*AE0154511,070,2462.A93.1% 96.7 / 0.7 29intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger related‑protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (13/14);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGC  >  AE015451/1070034‑1070440
                                                                                                                                                                                                                       ||                                                                                                                                                                                                  
tGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGT‑‑Gagatc                                                                                                                                                                                              <  1:283304/219‑5 (MQ=255)
                     aCCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAgcgcgg                                                                                                                                                                                        >  1:67165/1‑204 (MQ=255)
                                              gATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGa                                                                                                                                                         >  2:502526/1‑212 (MQ=255)
                                                     gTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                                                                                            <  1:102588/202‑1 (MQ=255)
                                                       aGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                                                                                            >  2:473631/1‑216 (MQ=255)
                                                                                      tcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                                                   <  2:426703/162‑1 (MQ=255)
                                                                                          tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                                                             >  1:2339/1‑148 (MQ=255)
                                                                                                   tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                                                                                    <  1:505942/196‑1 (MQ=255)
                                                                                                       ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                                                    >  1:339440/1‑224 (MQ=255)
                                                                                                          aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGcgcc                                                                                                                                                                                         >  1:159016/1‑120 (MQ=255)
                                                                                                          aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGcgcc                                                                                                                                                                                         <  2:159012/120‑1 (MQ=255)
                                                                                                                          tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGATGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCGCAgg                                                                                                                                                               >  2:197473/1‑130 (MQ=38)
                                                                                                                                gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAg                                               <  1:502530/236‑1 (MQ=255)
                                                                                                                                   ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAc                                            <  1:229788/236‑1 (MQ=255)
                                                                                                                                                  gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                                                                              <  1:275007/187‑1 (MQ=255)
                                                                                                                                                         aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcccgcacc                                                                         >  2:190603/1‑184 (MQ=255)
                                                                                                                                                               aGGGCTGCTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                                                     <  2:93587/103‑1 (MQ=25)
                                                                                                                                                                      aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                                                                                        <  1:302578/157‑1 (MQ=255)
                                                                                                                                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTGCGGCGCAGGGGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCa                                                                                   <  1:12688/159‑1 (MQ=38)
                                                                                                                                                                                     cagacgcaTGCCTTTAGGCTTTCATCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATat                               >  1:47178/8‑198 (MQ=37)
                                                                                                                                                                                        gggTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGTGCCGGGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                                                                                        >  2:480057/1‑139 (MQ=38)
                                                                                                                                                                                        gggTATGCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAAGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGc                                                                                     >  2:297769/1‑142 (MQ=255)
                                                                                                                                                                                        gggTATGCCGTTAGGCTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                                                                                             >  2:234588/1‑118 (MQ=37)
                                                                                                                                                                                              gCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgccagc  >  2:289712/1‑219 (MQ=255)
                                                                                                                                                                                              gCCGTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgccagc  <  1:289716/219‑1 (MQ=255)
                                                                                                                                                                                                 gTTAGGTTTTCGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGgccag   <  1:341880/215‑1 (MQ=255)
                                                                                                                                                                                                         tccGGCGC‑‑CTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCCCAGGCGATGAAAGCCTCAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                                                                                    <  2:67162/92‑1 (MQ=21)
                                                                                                                                                                                                           cGGCGC‑‑CTGTGAGATCGAGCCCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg                                                                     <  1:238784/139‑1 (MQ=255)
                                                                                                                                                                                                            ggCGC‑‑CTGTGAGATCGATCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                                                                                      <  2:100592/73‑1 (MQ=18)
                                                                                                                                                                                                                       ||                                                                                                                                                                                                  
TGGGGCACCCGCTGACCAACTACCTGTGCCTGGCGTTCATCGTGCTGATCCTGGTAGTGATGTACCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑CTGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAATCGCTCCAGGCGGCCAGC  >  AE015451/1070034‑1070440

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