Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,716,872 C→A 23.2% intergenic (+130/+68) yphH → / ← glyA predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,716,8720CA23.2% 163.1 / 30.9 73intergenic (+130/+68)yphH/glyApredicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base A (8/9);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCT  >  Exported/1716806‑1716940
                                                                  |                                                                    
tGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                                                  >  1:87508/1‑71 (MQ=255)
  atgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATa                                                                <  1:276209‑M2/62‑1 (MQ=255)
  atgcgcgatGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATa                                                                <  1:5482606‑M2/62‑1 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                                                  <  1:64705/69‑1 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                                                                 <  1:5506823/70‑1 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                                                                 <  1:1271959/70‑1 (MQ=255)
  aTGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGAt                                                                 <  1:4987178/70‑1 (MQ=11)
  aTGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGAt                                                                 <  1:1233144/70‑1 (MQ=11)
  aTGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATa                                                                <  1:2054670/71‑1 (MQ=11)
     cgcTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                                                             >  1:1208984/1‑71 (MQ=255)
     cgcTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAg                                                              <  1:3703517/70‑1 (MQ=11)
     cgcTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAg                                                              <  1:4010410/70‑1 (MQ=11)
       cgatGCTTATCAGGCCTACAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                                                          <  1:2714889‑M2/67‑1 (MQ=255)
         atGCTTATCAGGCCTACCAGAAGATAGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt                                                         <  1:4203362‑M2/69‑1 (MQ=255)
         aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                                                  >  1:4984271/1‑62 (MQ=255)
         aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                                                  >  1:1349458/1‑62 (MQ=255)
         aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgct                                                         >  1:4413369/1‑69 (MQ=2)
         aCGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgct                                                         >  1:1436366/1‑69 (MQ=2)
          cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                                                         <  1:5436240/70‑1 (MQ=255)
            cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                                                       >  1:4071410/1‑70 (MQ=255)
            cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                                                       >  1:2513460/1‑70 (MQ=255)
            cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc                                                      >  1:3248722/1‑71 (MQ=255)
            cTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATATGCGATGTTTTGTAGGACTGATAAGGCGCTcg                                                       >  1:5395870/1‑69 (MQ=2)
               aTCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg                                                     <  1:3798215/69‑1 (MQ=2)
                 cAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGcgccg                                                  >  1:934475/1‑70 (MQ=2)
                     ccTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                                                                 <  1:5218794/51‑1 (MQ=255)
                       tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                                                <  1:529082/50‑1 (MQ=255)
                       tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                                                <  1:1792834/50‑1 (MQ=255)
                         cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                                               >  1:1550923/1‑49 (MQ=255)
                         cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                                              >  1:1556834/1‑50 (MQ=255)
                         cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                                              >  1:4125468/1‑50 (MQ=255)
                         cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                           <  1:5050166/69‑1 (MQ=255)
                          acGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATaa                                                               >  1:5061555/1‑48 (MQ=255)
                          acGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                        >  1:5566762/1‑71 (MQ=255)
                           agaagaTTGCAATATATTGAATTTGCACTGTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                       >  1:1023432/4‑70 (MQ=11)
                            ggAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                                            >  1:3973499/1‑49 (MQ=255)
                              agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                    <  1:1869594/71‑1 (MQ=11)
                              agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                    <  1:4980081/71‑1 (MQ=11)
                              agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                    <  1:5502856/71‑1 (MQ=11)
                              agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                    <  1:1028001/71‑1 (MQ=11)
                               caTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaact                                   <  1:5329792/70‑2 (MQ=11)
                               caTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaac                                    >  1:3731757/2‑70 (MQ=11)
                                aTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATcaacaa                                  <  1:3874039/71‑1 (MQ=11)
                                aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                                                <  1:3329907/41‑1 (MQ=255)
                                  tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                 >  1:5574960/1‑70 (MQ=255)
                                  tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                 >  1:1650007/1‑70 (MQ=255)
                                    cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg                              >  1:3828834/1‑71 (MQ=255)
                                     aaTATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGc                             <  1:3284311/71‑1 (MQ=255)
                                      atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                            <  1:5076453/71‑1 (MQ=255)
                                      atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                            <  1:4537113/71‑1 (MQ=255)
                                      atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                            <  1:3494598/71‑1 (MQ=255)
                                       tctatTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATcaacaa                                  >  1:2837279/3‑64 (MQ=11)
                                       tctatTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATcaacaa                                  >  1:882821/3‑64 (MQ=11)
                                         tatTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaac                                    >  1:1107182/1‑60 (MQ=11)
                                         tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                                            >  1:2146612/1‑36 (MQ=40)
                                         tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca                          <  1:4737930/70‑1 (MQ=255)
                                         tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAt                         <  1:575186/71‑1 (MQ=255)
                                         ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCgg                                           >  1:1595834/3‑52 (MQ=0)
                                         ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCgg                                           >  1:2539797/3‑52 (MQ=0)
                                         ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCgg                                           >  1:1909925/3‑52 (MQ=0)
                                         ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGc                                          <  1:2514547/51‑1 (MQ=0)
                                           ttGAATTTGCGATGTTTTGTAGGACGGATAAGGCg                                                           <  1:4648181/35‑1 (MQ=14)
                                           ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCg                                                           >  1:1378075/1‑35 (MQ=40)
                                           ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCg                                                           >  1:2302650/1‑35 (MQ=40)
                                             gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                      >  1:3443282/1‑70 (MQ=255)
                                              gaTTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCaa                                     >  1:3873219/2‑53 (MQ=11)
                                                  tGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCa                 >  1:5128176/1‑70 (MQ=255)
                                                       aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc           >  1:1399563/1‑71 (MQ=255)
                                                       aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc           >  1:1091837/1‑71 (MQ=255)
                                                           ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt        <  1:1625196/70‑1 (MQ=255)
                                                           ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt        <  1:1857053/70‑1 (MQ=255)
                                                           ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt        <  1:4691661/70‑1 (MQ=255)
                                                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCt  <  1:810276/71‑1 (MQ=255)
                                                                  |                                                                    
TGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCT  >  Exported/1716806‑1716940

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: