Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,499,506:1 +C 100% intergenic (+83/‑81) trpS → / → zapE Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,499,5061.C97.1% 100.9 / ‑2.9 34intergenic (+83/‑81)trpS/zapETryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (9/24);  total (9/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.16e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGGCACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACA  >  AE015451/1499273‑1499713
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                
tGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAcg‑c                                                                                                                                                                                                                 <  1:487566/236‑1 (MQ=255)
                               aGCTCCGGCCTGTGGCTCACTGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAcg‑c                                                                                                                                                                                                                 <  1:101707/205‑1 (MQ=38)
                                                cAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                                                                                               <  1:309814/222‑1 (MQ=255)
                                                   gggCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGGGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGAc                                                                                                                                                                                                   <  2:792600/198‑1 (MQ=37)
                                                                 ccGACGGGCGGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACGGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAcg‑c                                                                                                                                                                                                                 <  2:95702/171‑1 (MQ=37)
                                                                                                       gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                                                                                           <  1:728461/235‑1 (MQ=39)
                                                                                                             aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCCTTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc                                                                                                                                               >  1:220505/1‑193 (MQ=37)
                                                                                                                 agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                                                                                                                            <  2:220503/192‑1 (MQ=37)
                                                                                                                  gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc                                                                                                  <  2:311518/233‑1 (MQ=255)
                                                                                                                     ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                                                                                                  <  2:688947/182‑1 (MQ=37)
                                                                                                                     ccTGAACGAAGCCGTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAcg‑c                                                                                                                                                                                                                 <  2:241790/119‑1 (MQ=21)
                                                                                                                        gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                                                                                                                >  1:667224/1‑165 (MQ=25)
                                                                                                                            gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTGTTGTTGCCTCGCTAACGa                                                                                                                                    >  2:703628/1‑189 (MQ=37)
                                                                                                                                         cacaGCAGGACGGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTGTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                                                                                     <  2:23683/191‑1 (MQ=37)
                                                                                                                                               agGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                                                                                       >  1:7882/1‑199 (MQ=37)
                                                                                                                                                   cTGAGGCTGCCCAACGGTGGGGGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAcg‑c                                                                                                                                                                                                                 <  2:256/89‑1 (MQ=17)
                                                                                                                                                   cTGAGGCGGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                                                                                <  2:217445/170‑1 (MQ=25)
                                                                                                                                                     gcggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACt                                                                                                                                                                                                  <  2:46456/100‑1 (MQ=11)
                                                                                                                                                      cggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCGCCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                                                                                     <  1:144636/177‑1 (MQ=37)
                                                                                                                                                      cggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCACG‑CCAGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                                                                                <  2:307070/166‑1 (MQ=25)
                                                                                                                                                        gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                                                                                                  <  1:709338/163‑1 (MQ=25)
                                                                                                                                                               aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGt                                                    >  2:294428/1‑234 (MQ=38)
                                                                                                                                                                   gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                                                                                               >  1:278405/1‑107 (MQ=14)
                                                                                                                                                                   gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                                                                                               <  2:278403/107‑1 (MQ=14)
                                                                                                                                                                          cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                                                                                     <  1:120452/158‑1 (MQ=37)
                                                                                                                                                                                 gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTg                                 >  1:712685/1‑235 (MQ=255)
                                                                                                                                                                                  cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTa                               <  1:799684/236‑1 (MQ=255)
                                                                                                                                                                                              gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCGCGCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                                                                                                  <  2:231280/77‑1 (MQ=11)
                                                                                                                                                                                              gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTgg                    <  1:82425/235‑1 (MQ=255)
                                                                                                                                                                                                  cgcgGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTcc                                         <  1:391548/210‑1 (MQ=39)
                                                                                                                                                                                                   gcgGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                                                                                                  <  2:369062/120‑1 (MQ=25)
                                                                                                                                                                                                       tGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTc                                          <  2:667222/204‑1 (MQ=39)
                                                                                                                                                                                                                ccGGCTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAacaaca  <  1:294430/234‑1 (MQ=255)
                                                                                                                                                                                                                    cTTTGCCGGTGTTCG‑CGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAATCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGa                                                             >  1:137167/1‑172 (MQ=255)
                                                                                                                                                                                                                    cTTTGCCGGTGTT‑CGCGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGttc                              >  1:422057/1‑201 (MQ=255)
                                                                                                                                                                                                                                ttcgcGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAgc                                      <  2:7882/180‑1 (MQ=255)
                                                                                                                                                                                                                                 tcgcGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCAcgcg            <  1:599355/206‑1 (MQ=255)
                                                                                                                                                                                                                                   gcGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGc                 <  1:171153/201‑1 (MQ=255)
                                                                                                                                                                                                                                    cGGGCACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCg                                                      <  2:115712/164‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                
TGGCGAGGCCGATGTCCGCATCGAAGGCCTGAGCTTCGGCCTGTGGCTCAATGGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGGCACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGCAGGAAACTGCGGTGCGTCACCTCCAGCGCCTGTACGACGACCTGGTGCACGCGCAGAACAACA  >  AE015451/1499273‑1499713

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