Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 533,863 T→C 5.9% Y525C (TAT→TGT) ‡ ybhF ← fused predicted transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported533,8630TC 5.9% 198.9 / 2.2 68Y525C (TAT→TGT) ‡ybhFfused predicted transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (37/27);  new base C (2/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

ATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTG  >  Exported/533802‑533933
                                                             |                                                                      
aTCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCt                                                                >  1:5619433/1‑70 (MQ=255)
      cGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCa                                                         >  1:2966709/1‑71 (MQ=255)
      cGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCa                                                         <  1:4462223/71‑1 (MQ=255)
        gCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCAt                                                        >  1:898690/1‑70 (MQ=255)
          gTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATaa                                                      <  1:2731463/70‑1 (MQ=255)
          gTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATa                                                       >  1:1264724/1‑69 (MQ=255)
          gTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATa                                                       >  1:5583975/1‑69 (MQ=255)
          gTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATa                                                       <  1:4489416/69‑1 (MQ=255)
             ccGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAgtg                                                  >  1:3844657/1‑71 (MQ=255)
             ccGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCAAAAAgt                                                   >  1:3957901/1‑70 (MQ=255)
                  ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAgtgggtgg                                             >  1:1589497/1‑71 (MQ=255)
                   gCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGggtggt                                            <  1:5581305/71‑1 (MQ=255)
                       ttAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTc                                                                    <  1:1840173/43‑1 (MQ=255)
                       ttAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCa                                                             >  1:4450844/1‑50 (MQ=255)
                       ttAATTTCACGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGAc                                         >  1:3068870/1‑70 (MQ=255)
                        tAATTTCCCGTGGTACACCAGGCCGATGCGGTCGCAat                                                                        >  1:533025/1‑38 (MQ=255)
                        tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAata                                                                       >  1:4129405/1‑39 (MQ=255)
                        tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAg                                                    >  1:3297072/1‑58 (MQ=255)
                        tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAg                                                    >  1:353430/1‑58 (MQ=255)
                        tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCa                                       >  1:5470136/1‑71 (MQ=255)
                         aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCg                                                                  >  1:106375/1‑43 (MQ=255)
                         aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCAt                                      >  1:1894517/1‑71 (MQ=255)
                          aTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCAt                                      >  1:546633/1‑70 (MQ=255)
                          aTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCAt                                      >  1:4060162/1‑70 (MQ=255)
                           tttCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAGGTGGGTGGTGAc                                         <  1:2905679/66‑1 (MQ=255)
                           tttCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGAcc                                        <  1:4065994/67‑1 (MQ=255)
                           tttCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCa                                    >  1:5628321/1‑71 (MQ=255)
                             tCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCAcc                                  >  1:4808392/1‑71 (MQ=255)
                             tCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCAcc                                  >  1:3003380/1‑71 (MQ=255)
                               ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                <  1:4330453/71‑1 (MQ=255)
                               ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                <  1:945889/71‑1 (MQ=255)
                               ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                <  1:1008892/71‑1 (MQ=255)
                               ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                <  1:3236423/71‑1 (MQ=255)
                                 gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCAt                                      >  1:4346371/1‑63 (MQ=255)
                                 gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCAt                                      >  1:774054/1‑63 (MQ=255)
                                 gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCACCCATAAAGTGggtggt                                            >  1:3178283/1‑57 (MQ=255)
                                  cgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGAc                             <  1:3289606/71‑1 (MQ=255)
                                  cgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGAc                             <  1:3780911/71‑1 (MQ=255)
                                     tACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGcc                          <  1:162026/71‑1 (MQ=255)
                                     tACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGcc                          <  1:2587750/71‑1 (MQ=255)
                                       cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATc                                     <  1:2085227/58‑1 (MQ=255)
                                        acCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCttt                       <  1:2252577/71‑1 (MQ=255)
                                        acCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCttt                       <  1:1976867/71‑1 (MQ=255)
                                          cAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTtct                     >  1:2367174/1‑71 (MQ=255)
                                           aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAgtgggt                                               >  1:2262723/1‑44 (MQ=255)
                                           aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTtctc                    >  1:1463739/1‑71 (MQ=255)
                                           aGGCCGATGCGGTCGCAATATTCCCTTCATCCATaa                                                      <  1:4243377/36‑1 (MQ=38)
                                            ggCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                >  1:2096646/1‑58 (MQ=255)
                                               cGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACg                            >  1:4365044/1‑59 (MQ=255)
                                                gATGCGGTCGCAACGTTCCGCTTCGTTCATAAAgtggg                                                >  1:4480777/1‑38 (MQ=21)
                                                gATGCGGTCGCAACGTTCCGCTTCGTTCATAAAgtggg                                                >  1:4929905/1‑38 (MQ=21)
                                                 aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                <  1:3458725/53‑1 (MQ=255)
                                                 aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGt                                <  1:1641699/53‑1 (MQ=255)
                                                  tGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGgagatg                                        <  1:5382951/44‑3 (MQ=19)
                                                  tGCGGTCGCAACGTTCCGCTTCGTTCATAAAGTGGGTGgagatg                                        <  1:2428844/44‑3 (MQ=19)
                                                   gCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGa                              >  1:4748028/1‑53 (MQ=255)
                                                   gCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGc            >  1:3752177/1‑71 (MQ=255)
                                                      gTCGCAATATTCCGCTTCATCCATAAAgtgggtgg                                             >  1:4195918/1‑35 (MQ=255)
                                                      gTCGCAATATTCCGCTTCATCCATAAAgtgggtgg                                             >  1:3287737/1‑35 (MQ=255)
                                                      gTCGCAATATTCCGCTTCATCCATAAAgtgggtgg                                             >  1:752131/1‑35 (MQ=255)
                                                      gTCGCAATATTCCGCTTCATCCATAAAgtgggtgg                                             >  1:2139921/1‑35 (MQ=255)
                                                       tCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGtt        <  1:1465206/71‑1 (MQ=255)
                                                       tCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGtt        <  1:1727380/71‑1 (MQ=255)
                                                       tCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGtt        >  1:256481/1‑71 (MQ=255)
                                                        cGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGtt        >  1:2831203/1‑70 (MQ=255)
                                                        cGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGtt        >  1:3238625/1‑70 (MQ=255)
                                                         gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTg       <  1:2869320/70‑1 (MQ=255)
                                                         gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTg       <  1:5446902/70‑1 (MQ=255)
                                                         gCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGa      <  1:4395727/71‑1 (MQ=255)
                                                          aaaTATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGAc                             <  1:641677/46‑1 (MQ=255)
                                                             tatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGAtgtg  <  1:1016677/71‑1 (MQ=255)
                                                             |                                                                      
ATCGTCCGGCGTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTG  >  Exported/533802‑533933

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: