Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,294,017 T→A 22.3% intergenic (+44/‑231) gldA → / → yijF glycerol dehydrogenase, NAD/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,294,0170TA22.3% 95.5 / 22.7 49intergenic (+44/‑231)gldA/yijFglycerol dehydrogenase, NAD/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base A (3/8);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

AGCGTTTCCTGCAAGAGTGGGAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTG  >  Exported/2293948‑2294083
                                                                     |                                                                  
aGCGTTTCCTGCAAGAGTGGGAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAAt                                                                    >  1:3363786/1‑70 (MQ=255)
             agagTGGGAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTa                                                      >  1:5656924/1‑71 (MQ=255)
                 tGGGAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCt                                                  <  1:3310428/71‑1 (MQ=255)
                   ggAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCtt                                                 <  1:3205349/70‑1 (MQ=255)
                   ggAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTa                                                <  1:1404335/71‑1 (MQ=255)
                   ggAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTa                                                <  1:3821979/71‑1 (MQ=255)
                   ggAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTa                                                >  1:4039048/1‑71 (MQ=255)
                        aaCCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAgg                                           >  1:2207878/1‑71 (MQ=255)
                         aCCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                          >  1:791442/1‑71 (MQ=255)
                             aCTCCAATCTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                      >  1:5379214/1‑71 (MQ=255)
                                 cAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCaa                                  >  1:2816222/1‑71 (MQ=255)
                                    aCTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAtct                               >  1:1456124/1‑71 (MQ=255)
                                     cTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                   <  1:3503152/66‑1 (MQ=255)
                                     cTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                   <  1:5360329/66‑1 (MQ=255)
                                     cTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCa                                   <  1:147098/66‑1 (MQ=255)
                                          ggCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCa                          >  1:4320852/1‑70 (MQ=255)
                                           gCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTa                                                <  1:3322978/47‑1 (MQ=255)
                                             ttGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCaa                         <  1:4947334/68‑1 (MQ=255)
                                             ttGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAata                       >  1:4860091/1‑70 (MQ=255)
                                             ttGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAata                       >  1:1397808/1‑70 (MQ=255)
                                              tGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACCCTTATCAg                                            <  1:5208743/48‑1 (MQ=255)
                                               gTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatatat                    >  1:2165227/1‑71 (MQ=255)
                                                   gggAGTTTGAAGGCAAAATTGCCTGATGCGCTACGCTTTTCAGGCCt                                        <  1:1375929/47‑1 (MQ=17)
                                                   gggAGTTTGAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTg                  <  1:2340155/69‑1 (MQ=255)
                                                    ggAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTg                  <  1:650623/68‑1 (MQ=255)
                                                      aGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAt               <  1:3568451/69‑1 (MQ=255)
                                                       gTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTg            >  1:2191214/1‑71 (MQ=255)
                                                        tttGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAatata                     <  1:327036/61‑1 (MQ=255)
                                                         ttGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAt                                 >  1:3972747/1‑48 (MQ=255)
                                                         ttGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAAt                                 >  1:5078741/1‑48 (MQ=255)
                                                          tGAACGCAAAATTGCCTGATGCGCTACGCTTATCa                                             >  1:1166492/1‑35 (MQ=37)
                                                          tGAACGCAAAATTGCCTGATGCGCTACGCTTATCa                                             >  1:2506994/1‑35 (MQ=37)
                                                              cGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCtt     >  1:3491545/1‑71 (MQ=255)
                                                              cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCaa                         <  1:3392733/51‑1 (MQ=14)
                                                              cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCtt     <  1:2258950/71‑1 (MQ=255)
                                                              cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCtt     <  1:3711802/71‑1 (MQ=255)
                                                              cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCtt     <  1:3338977/71‑1 (MQ=255)
                                                              cGCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCtt     <  1:4396510/71‑1 (MQ=255)
                                                               gCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAat                        >  1:4557979/1‑51 (MQ=14)
                                                               gCAAAAATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCttt    <  1:4666928/71‑1 (MQ=255)
                                                                cAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGc           >  1:3312632/1‑63 (MQ=255)
                                                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTg  >  1:5264793‑M2/1‑70 (MQ=255)
                                                                 aaaaATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATctct                             >  1:1819627/1‑44 (MQ=14)
                                                                 aaaaATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGATTTTGCGTGCTTTTg  >  1:5198063‑M2/1‑70 (MQ=255)
                                                                 aaaaATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg  <  1:5487399‑M2/71‑2 (MQ=255)
                                                                 aaaaATGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg  <  1:1393426‑M2/71‑2 (MQ=255)
                                                                  aaaaTGCCTGATGCGCTACGCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTg  <  1:44739‑M2/70‑2 (MQ=255)
                                                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt             <  1:4161666/59‑1 (MQ=255)
                                                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAAttt             <  1:4476929/59‑1 (MQ=255)
                                                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTg            <  1:5504369/60‑1 (MQ=255)
                                                                   aaTTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGa                 <  1:4782298/54‑1 (MQ=255)
                                                                     |                                                                  
AGCGTTTCCTGCAAGAGTGGGAATAACCTACTCCAAACTCCCGGCTTGTCCGGGAGTTTGAACGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTG  >  Exported/2293948‑2294083

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