Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,155,390:1 +A 100% intergenic (‑181/+43) trnY‑Lys ← / ← purA tRNA‑Lys/adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,155,3901.A97.4% 128.3 / ‑3.0 39intergenic (‑181/+43)trnY‑Lys/purAtRNA‑Lys/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (15/23);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-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.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAAC‑AAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAG  >  BS168_reference_genome/4155259‑4155519
                                                                                                                                    |                                                                                                                                  
aTCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCaaac‑aaa‑aaa                                                                                                                             >  3:273736/1‑137 (MQ=255)
     gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCaaac‑aaa‑aaacaa                                                                                                                          <  5:30659/136‑5 (MQ=255)
     gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCaaac‑aaa‑aaacaa                                                                                                                          <  8:54182/136‑5 (MQ=255)
     gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCaaac‑aaa‑aaacaa                                                                                                                          <  7:261499/136‑5 (MQ=255)
     gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGca                                                                                                                        >  1:221457/1‑138 (MQ=255)
     gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGc                                                                                                                         <  4:273736/138‑1 (MQ=255)
     gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATAACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCaaac‑aaa‑aaacaa                                                                                                                          <  7:172374/136‑5 (MQ=255)
     gTAATCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCaaac‑aaa‑aaacaa                                                                                                                          <  5:34552/136‑5 (MQ=255)
          ctATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCtt                                                                                                                    >  6:136524/1‑138 (MQ=255)
           tATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTa                                                                                                                   >  8:244958/1‑138 (MQ=255)
            aTTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTaaa                                                                                                                 >  6:9819/1‑139 (MQ=255)
               gAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATa                                                                                                               <  2:101572/138‑1 (MQ=255)
               gAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATa                                                                                                               <  8:128020/138‑1 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAATTAGAAGTGGCTTTAACGTTTTTGGTTTTTTTTCATTTTTTGTTGTCAGGCCGACTTATATGCGGCTTTGAAAC‑AAATAACTAGCGCGTTAATGCGGGGCtg                                                                                                        >  2:193242/1‑138 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGcccc                                                                                                                      <  3:17503/125‑1 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGcccc                                                                                                                      >  4:17503/1‑125 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCtt                                                                                                        >  8:169457/1‑139 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCtt                                                                                                        >  3:165635/1‑139 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCtt                                                                                                        >  4:52641/1‑139 (MQ=255)
                     gTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCtt                                                                                                        >  1:49959/1‑139 (MQ=255)
                      tttCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTg                                                                                                       <  6:149510/139‑1 (MQ=255)
                         ccGTAACAACGAGATTTATAATATCACGAACTAATAGAAGTGTCAATTACTTTTTTGTTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATCCGGCTATAAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCa                                                                                                     <  3:171204/138‑1 (MQ=255)
                           gTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGAc                                                                                                  <  1:167128/139‑1 (MQ=255)
                                aacGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATAtt                                                                                             <  6:179242/139‑1 (MQ=255)
                                  cGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGAc                                                                                                  <  4:104994/132‑1 (MQ=255)
                                  cGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGAc                                                                                                  >  3:104994/1‑132 (MQ=255)
                                       tttATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTAtt                                                                                        >  4:250425/1‑137 (MQ=255)
                                           taataTCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTc                                                                                  <  7:35398/139‑1 (MQ=255)
                                             atatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc                                                                                <  6:193375/139‑1 (MQ=255)
                                             atatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc                                                                                <  7:169457/139‑1 (MQ=255)
                                             atatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc                                                                                <  7:244958/139‑1 (MQ=255)
                                             atatAACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGc                                                                                <  6:141224/139‑1 (MQ=255)
                                              tatCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCa                                                                               <  5:9819/139‑1 (MQ=255)
                                                tCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACatat                                                                                              <  8:269239/122‑1 (MQ=255)
                                                tCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACatat                                                                                              >  7:269239/1‑122 (MQ=255)
                                                tCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCAc                                                                              <  5:281723/138‑1 (MQ=255)
                                                tCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCAc                                                                              <  3:52641/138‑1 (MQ=255)
                                                       cAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTacaca                                                                      <  1:194898/139‑1 (MQ=255)
                                                             gAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCg                                                                 <  6:293808/138‑1 (MQ=255)
                                                                               ttGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGa                                               <  8:12333/138‑1 (MQ=255)
                                                                                tGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAg                                              <  4:60728/138‑1 (MQ=255)
                                                                                                           aCCAACTTATAGGCGGCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATGTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGa                    >  8:90815/1‑137 (MQ=255)
                                                                                                                          gCTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGt     >  6:18463/1‑137 (MQ=255)
                                                                                                                           cTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAg  >  6:14159/1‑139 (MQ=255)
                                                                                                                           cTTTCAAACAAAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAg  <  5:14159/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
ATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAAC‑AAA‑AACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAG  >  BS168_reference_genome/4155259‑4155519

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