Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,807,677 T→A 5.7% intergenic (+496/‑611) tyrB → / → aphA tyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/acid phosphatase/phosphotransferase, class B, non‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,807,6770TA 5.7% 261.4 / 6.1 88intergenic (+496/‑611)tyrB/aphAtyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/acid phosphatase/phosphotransferase, class B, non‑specific
Reads supporting (aligned to +/- strand):  ref base T (28/55);  new base A (3/2);  total (31/57)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 106 total aligned reads displayed.

TTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCC  >  Exported/2807612‑2807740
                                                                 |                                                               
tttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTaa                                                            >  1:1184996/1‑71 (MQ=255)
tttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTaa                                                            >  1:2091152/1‑71 (MQ=255)
tttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTATTACCCCTGATTTTGTTa                                                             >  1:842982/1‑70 (MQ=255)
 ttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACTTAAAGTCTTTCTTACCCCTGATTTTGTTAAc                                                           >  1:3079201/1‑71 (MQ=255)
 ttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTaa                                                            >  1:4833462/1‑70 (MQ=255)
 ttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAAc                                                           >  1:884539/1‑71 (MQ=255)
 ttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAAc                                                           >  1:860081/1‑71 (MQ=255)
  tAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAAc                                                           >  1:1307409/1‑70 (MQ=255)
        cctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATc                                                     >  1:1042096/1‑70 (MQ=255)
         ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCa                                                    <  1:3461027/70‑1 (MQ=255)
         ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCa                                                    <  1:3880453/70‑1 (MQ=255)
            ggCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:2954011/70‑1 (MQ=255)
            ggCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:2247133/70‑1 (MQ=255)
            ggCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:4797662/70‑1 (MQ=255)
            ggCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 >  1:4991946/1‑70 (MQ=255)
             gCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATaaa                                               >  1:2541822/1‑71 (MQ=255)
              cAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTAGTTAACAAAATCag                                                   <  1:1559923/66‑2 (MQ=39)
              cAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTAGTTAACAAAATCag                                                   <  1:1860764/66‑2 (MQ=39)
               aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAAt                                              <  1:3369173/70‑1 (MQ=255)
               aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAAt                                              <  1:2355216/70‑1 (MQ=255)
               aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAAt                                              <  1:222998/70‑1 (MQ=255)
               aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATg                                             <  1:1123420/71‑1 (MQ=255)
               aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATg                                             <  1:2239347/71‑1 (MQ=255)
               aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATg                                             >  1:1421595/1‑71 (MQ=255)
                aaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTAGTTAACAAAATCaggg                                                 >  1:4900492/1‑63 (MQ=38)
                aaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTAGTTAACAAAATCaggg                                                 >  1:3993938/1‑63 (MQ=38)
                aaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTAGTTAACAAAATCaggg                                                 >  1:3884995/1‑63 (MQ=38)
                      cAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:4619558/71‑1 (MQ=255)
                       aaaGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:596546/70‑1 (MQ=255)
                       aaaGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:726659/70‑1 (MQ=255)
                        aaGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaaaa                                    >  1:3855125/1‑71 (MQ=255)
                        aaGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaaaa                                    >  1:237995/1‑71 (MQ=255)
                        aaGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaaaa                                    >  1:3352523/1‑71 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:40198/53‑1 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:413336/53‑1 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:5135548/53‑1 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:3326187/53‑1 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:1530057/53‑1 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:96929/53‑1 (MQ=255)
                             tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATa                                                 <  1:1705323/53‑1 (MQ=255)
                                  atatACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACttt                                <  1:2152533/65‑1 (MQ=255)
                                  atatACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACttt                                <  1:4091490/65‑1 (MQ=255)
                                  atatACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACttt                                <  1:1511756/65‑1 (MQ=255)
                                  atatACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACttt                                <  1:1001483/65‑1 (MQ=255)
                                  atatACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACttt                                <  1:5169670/65‑1 (MQ=255)
                                   tataCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAAt                                              <  1:3967903/50‑1 (MQ=255)
                                   tataCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGc                          <  1:5113069/70‑1 (MQ=255)
                                   tataCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGc                          >  1:2693855/1‑70 (MQ=255)
                                    ataCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:76746/1‑71 (MQ=255)
                                     taCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:3133547/1‑70 (MQ=255)
                                     taCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:4576981/1‑70 (MQ=255)
                                        aTAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGa                    <  1:3073258/71‑1 (MQ=255)
                                           aaGTCTTTCTTACCCCTGATTTTGTTAACTAAATCaa                                                   >  1:1564703/1‑37 (MQ=255)
                                           aaGTCTTTCTTACCCCTGATTTTGTTAACTAAATCaa                                                   >  1:1858582/1‑37 (MQ=255)
                                           aaGTCTTTCTTACCCCTGATTTTGTTAACTAAATCaa                                                   >  1:4986055/1‑37 (MQ=255)
                                           aaGTCTTTCTTACCCCTGATTTTGTTAACTAAATCaa                                                   >  1:1314260/1‑37 (MQ=255)
                                           aaGTCTTTCTTACCCCTGATTTTGTTAACTAAATCaa                                                   >  1:1239406/1‑37 (MQ=255)
                                           aaGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACt                 >  1:2484057/1‑71 (MQ=255)
                                            aGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACtt                <  1:1486589/71‑1 (MQ=255)
                                            aGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACtt                <  1:4959571/71‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATTCGTGaa                                       <  1:671792/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:4642771/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:1983521/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:4686378/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:1933556/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:4864264/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:2390636/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:2113102/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:816723/47‑1 (MQ=255)
                                             gTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaa                                       <  1:3981786/47‑1 (MQ=255)
                                              tctttcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              <  1:3648997/71‑1 (MQ=255)
                                              tctttcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              <  1:1645810/71‑1 (MQ=255)
                                                tttcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaaaa                                    >  1:176205/1‑47 (MQ=255)
                                                 ttcttTCCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACtt                <  1:4198450/66‑1 (MQ=255)
                                                      aCCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATActtc        <  1:2192259/69‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:1907757/35‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:4918473/35‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:4024857/35‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:4089831/35‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:4055274/35‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGaaa                                      <  1:4046447/35‑1 (MQ=255)
                                                          cTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTtctc   >  1:823157/1‑70 (MQ=255)
                                                           tGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTa                               >  1:2506620/1‑41 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTAc                              >  1:28427/1‑38 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTAc                              >  1:2945526/1‑38 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTAc                              >  1:3867339/1‑38 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTAc                              >  1:3159225/1‑38 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTAc                              >  1:536646/1‑38 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTAc                              >  1:5138751/1‑38 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:40076/1‑44 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:1390589/1‑44 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:2164906/1‑44 (MQ=255)
                                                               tttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCaa                        >  1:4212356/1‑44 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGTGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:2343507/1‑53 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:526947/1‑53 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:3480723/1‑53 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:3109206/1‑53 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:4097187/1‑53 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:1572684/1‑53 (MQ=255)
                                                                ttGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa              >  1:3501528/1‑53 (MQ=255)
                                                                 |                                                               
TTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCC  >  Exported/2807612‑2807740

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