Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,586,056:1 +C 100% intergenic (+139/+76) PP_4061 → / ← PP_4063 conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,586,0561.C93.9% 105.9 / 1.2 33intergenic (+139/+76)PP_4061/PP_4063conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (14/17);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGATGCGTGCCTTGCACCAGCC  >  AE015451/4585849‑4586281
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
ccGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATgcg                                                                                                                                                                                                                <  2:1151825/229‑1 (MQ=255)
                             gggCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCt                                                                                                                                                                                                                        >  1:613334/1‑192 (MQ=255)
                              ggCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                                                                            <  1:1007033/235‑1 (MQ=255)
                                  aCCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGCC‑C‑GTGCCGGCCTTGaa                                                                                                                                                                                                                    >  1:1255544/1‑191 (MQ=255)
                                                                 ttAAGGGTAGACCACGGCACCGGCGTGAGCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGcccc                                                                                                                                                                                                      <  1:667532/174‑1 (MQ=255)
                                                                      ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                                                                                                                       <  1:736052/200‑1 (MQ=255)
                                                                      ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                                                                                                                       >  2:736049/1‑200 (MQ=255)
                                                                       gTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGcccc                                                                                                                                                                                                      >  1:3784/1‑168 (MQ=39)
                                                                           aCCACGCCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                                                                                                                       <  2:769395/195‑1 (MQ=39)
                                                                            ccACGGCACCGGCGTGACCAGCTGCTCGGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGGGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                                                                            <  2:659639/189‑1 (MQ=38)
                                                                            ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                                                                                                                 >  1:878807/1‑168 (MQ=39)
                                                                               cGGCGCTGGCGTGAACAGCTGATCTGCGCTGCTCCATGGAAAGGTGGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATgcggcg                                                                                                                                                                                                             <  2:76051/153‑1 (MQ=38)
                                                                               cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATgcgg                                                                                                                                                                                                               <  2:102930/151‑1 (MQ=39)
                                                                                  cACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCca                                                                                                                        >  1:583826/1‑235 (MQ=255)
                                                                                                    aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATgcggcg                                                                                                                                                                                                             >  2:753500/1‑132 (MQ=38)
                                                                                                                                gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCCTGCgg                                                                                                    >  1:295748/1‑209 (MQ=39)
                                                                                                                                      cGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCcac                                                                                                                       <  2:1224675/184‑1 (MQ=255)
                                                                                                                                                gCAAGGTGAGGATCTTGCATTGACCGGGCCGGCCTCTTCGCGGGTAAACCCGCTGCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGttt                                                                                                                                                                                      <  2:207943/111‑1 (MQ=25)
                                                                                                                                                   aGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACa                                                       <  1:747208/235‑1 (MQ=255)
                                                                                                                                                      tGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTg                                                                                                                                                                                                                      >  2:357414/1‑73 (MQ=18)
                                                                                                                                                                  aTTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                                                                                                                 <  2:1047957/82‑1 (MQ=18)
                                                                                                                                                                  aTTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGc                                                                                                                                   <  2:968808/144‑1 (MQ=37)
                                                                                                                                                                  aTTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGC‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCCGGTAGGAGCGGGtt                                                                                                                                                                                       >  1:757027/1‑92 (MQ=21)
                                                                                                                                                                   ttGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGCTGCAGCCGCAATCTCCGCCACACTGATCTCACg                                                                                                           >  1:230563/1‑167 (MQ=38)
                                                                                                                                                                      aTCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATAcc                                          <  2:583823/229‑1 (MQ=255)
                                                                                                                                                                         gggCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑ACC‑GGGCCGGCCTTGAATGCGGCGCTGCCCCTggaa                                                                                                                                                                                                 >  2:138594/1‑73 (MQ=17)
                                                                                                                                                                             cggccTCTTCGCGGGTAAACCCGCTCCTACCAGGG‑A‑CTGTGCCGGCCTTGAATgc                                                                                                                                                                                                                 >  2:1278851/1‑55 (MQ=11)
                                                                                                                                                                               gcctcttcgcgggTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAgt                           <  1:847979/224‑1 (MQ=255)
                                                                                                                                                                                     ccgcgggTAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATAcc                                          <  1:402578/209‑1 (MQ=255)
                                                                                                                                                                                            tAAACCCGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGa                                                         <  1:1063583/192‑1 (MQ=255)
                                                                                                                                                                                                 ccGCTCCTACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGATGCGTGCCTTGCa        >  1:960538/1‑236 (MQ=255)
                                                                                                                                                                                                        tACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGATGCGTGCCTTGCACCAGcc  >  1:771909/1‑235 (MQ=255)
                                                                                                                                                                                                        tACCAGGG‑ACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGATGCGTGCCTTGCACCAGcc  >  1:489877/1‑235 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  
CCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGG‑A‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGATGCGTGCCTTGCACCAGCC  >  AE015451/4585849‑4586281

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