Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,467,908 T→A 26.3% intergenic (+67/+61) ubiG → / ← yfaL bifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,467,9080TA26.3% 128.8 / 36.4 59intergenic (+67/+61)ubiG/yfaLbifunctional 3‑demethylubiquinone‑9 3‑methyltransferase and 2‑octaprenyl‑6‑hydroxy phenol methylase/adhesin
Reads supporting (aligned to +/- strand):  ref base T (26/17);  new base A (8/7);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTT  >  Exported/1467838‑1467977
                                                                          |                                                                       
gTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCt                                                                         <  1:1940781/71‑1 (MQ=255)
     gTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTcc‑‑cc                                                                   <  1:843942/70‑1 (MQ=255)
       tGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGt                                                                 >  1:2447871/1‑70 (MQ=255)
          tGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGc                                                             >  1:3144410/1‑71 (MQ=255)
          tGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGc                                                             <  1:2023096/71‑1 (MQ=255)
              tGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGa                                                         >  1:3473915/1‑71 (MQ=255)
                  cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAg                                                     <  1:3733454/71‑1 (MQ=255)
                  cgcgCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAg                                                     <  1:468940/71‑1 (MQ=255)
                    cgcCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGc                                                   <  1:150269/71‑1 (MQ=255)
                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGt                                                 >  1:5322616/1‑71 (MQ=255)
                      ccAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGt                                                 >  1:1534616/1‑71 (MQ=255)
                       cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt                                                <  1:1838129/71‑1 (MQ=255)
                       cAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt                                                <  1:647468/71‑1 (MQ=255)
                        aaTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGtt                                                >  1:125423/1‑70 (MQ=255)
                            ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGt                                                                 >  1:5569488/1‑49 (MQ=255)
                              tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg                                                                  <  1:2097987/46‑1 (MQ=255)
                                aTGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       >  1:2453231/1‑71 (MQ=255)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                            <  1:5374231/71‑1 (MQ=0)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                            >  1:4963371/1‑71 (MQ=0)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcg                                            >  1:4890722/1‑71 (MQ=0)
                                aTGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAc                                             <  1:4719260/70‑1 (MQ=0)
                                 tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg                                                                  >  1:37603/1‑43 (MQ=255)
                                 tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCg                                                                  >  1:1700515/1‑43 (MQ=255)
                                 tGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       <  1:4427823/70‑1 (MQ=255)
                                 tGCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTAcgc                                           >  1:1345902/1‑71 (MQ=0)
                                  gCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                          >  1:29656/1‑71 (MQ=0)
                                  gCGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                          >  1:4485858/1‑71 (MQ=0)
                                   cGACGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTAcgcc                                          <  1:4567717/70‑1 (MQ=0)
                                     aCGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                  >  1:1858644/1‑71 (MQ=255)
                                     aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       >  1:1204526/1‑71 (MQ=0)
                                     aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       >  1:3608400/1‑71 (MQ=0)
                                     aCGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       <  1:3533806/71‑1 (MQ=0)
                                      cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                      <  1:4720923/71‑1 (MQ=0)
                                      cGCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                      <  1:5289091/71‑1 (MQ=0)
                                       gCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                 >  1:3375370/1‑70 (MQ=255)
                                       gCTTGCGCGTCTTATCAGGCCTACAAATTCCCGCACCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATc                                     <  1:445947/71‑1 (MQ=0)
                                           tgcgcgTCTTATCAGGCCTACAAATTCCCGCACCCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                  >  1:3564206/3‑71 (MQ=2)
                                            acgcgTCTTATCAGGCCTACAAATTCCCGCGCCATCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                   <  1:2023343/68‑1 (MQ=0)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                         >  1:338519/1‑71 (MQ=255)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                         >  1:4990484/1‑71 (MQ=255)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                         >  1:5225554/1‑71 (MQ=255)
                                              gcgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGc                         >  1:4659247/1‑71 (MQ=255)
                                               cgTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGcc                        <  1:3636583/71‑1 (MQ=255)
                                                   ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGa                      >  1:4731841/1‑69 (MQ=255)
                                                   ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                     >  1:2112006/1‑70 (MQ=255)
                                                   ttATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGAc                     >  1:5660861/1‑70 (MQ=255)
                                                             tACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGAt          >  1:545890/1‑71 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAc                              >  1:2670849/1‑49 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa        >  1:1430356/1‑71 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa        >  1:5600342/1‑71 (MQ=255)
                                                               cAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTa        >  1:1377418/1‑71 (MQ=255)
                                                                aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       <  1:3969461/39‑1 (MQ=255)
                                                                aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCa                                       <  1:108964/39‑1 (MQ=255)
                                                                aaaT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTAc       >  1:695619/1‑71 (MQ=255)
                                                                  aT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAACGTGCCGACTAGACAGTGATTACCa     >  1:5178548/1‑71 (MQ=255)
                                                                   t‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAt    <  1:2224276/71‑1 (MQ=255)
                                                                         cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt  <  1:2905643/71‑1 (MQ=255)
                                                                         cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt  <  1:3436856/71‑1 (MQ=255)
                                                                         cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt  <  1:3545106/71‑1 (MQ=255)
                                                                         cTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCAttt  <  1:3545542/71‑1 (MQ=255)
                                                                          |                                                                       
GTGAGGTTGATGTTTGGCCGCGCCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAT‑‑‑‑GCTCC‑‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACCGTGCCGACTAGACAGTGATTACCATTT  >  Exported/1467838‑1467977

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: