Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,032,108 T→A L253F (TTA→TTT)  SAUSA300_RS05075 ← 1,4‑dihydroxy‑2‑naphthoate octaprenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,032,1080TA98.0% 137.2 / ‑4.4 50L253F (TTA→TTT) SAUSA300_RS050751,4‑dihydroxy‑2‑naphthoate octaprenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (27/21);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CTGGCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAATAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCATAAATGTTAAAG  >  NC_007793/1032017‑1032190
                                                                                           |                                                                                   
cTGGCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGcc                                                                             <  1:3540938/100‑1 (MQ=255)
cTGGCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGAAAgatcggaagag                                                                             >  1:2016372/1‑92 (MQ=255)
 ctgCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGACTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAA‑AATGAGccac                                                                           >  1:4053483/3‑100 (MQ=255)
 aggCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAAcca                                                                            >  1:4733111/2‑100 (MQ=255)
  ggCATCATTGTAGGCGGTGGATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGAAAGGAAaaaaatagatcgg                                                                           >  1:211430/1‑93 (MQ=255)
   gCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGccacc                                                                          <  1:5106882/100‑1 (MQ=255)
   gCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGccacc                                                                          >  1:2892257/1‑100 (MQ=255)
    catcatTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACACCGCAAGTAAGTAAACAAATGAAccacca                                                                         >  1:2466202/1‑100 (MQ=255)
    catcatTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGccacca                                                                         >  1:5078615/1‑100 (MQ=255)
     atcatTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGccaccgc                                                                        <  1:663751/100‑3 (MQ=255)
       catTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCATCAggc                                                                      >  1:559290/1‑99 (MQ=255)
       catTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAAGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAACGAGCCACCAgga                                                                      <  1:4651384/100‑2 (MQ=255)
         ttGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCTCCAGGttt                                                                    <  1:1222745/100‑1 (MQ=255)
         ttGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGttt                                                                    <  1:660797/100‑1 (MQ=255)
         ttGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGttt                                                                    <  1:4333858/100‑1 (MQ=255)
         ccctAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGttt                                                                    <  1:1904243/97‑1 (MQ=255)
          tGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTta                                                                   >  1:3160442/1‑100 (MQ=255)
          tGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGACCCACCAgatcgg                                                                   >  1:2063479/1‑95 (MQ=255)
            tAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCGCCAGGTTtaat                                                                 <  1:36460/100‑1 (MQ=255)
              ggCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTtaataa                                                               <  1:5636007/100‑1 (MQ=255)
                 ggTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAATAATGAGCCACCAGGTTTAATAAT‑Ga                                                           <  1:44940/100‑1 (MQ=255)
                      aTCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTa                                                      <  1:4574293/100‑1 (MQ=255)
                                 tGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAATAGATCGGAAgagcacacg                                            >  1:576739/1‑92 (MQ=255)
                                 tGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAATAGATCGGAAgagcacacg                                            >  1:663842/1‑92 (MQ=255)
                                 tGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATaaa                                           >  1:5610069/1‑100 (MQ=255)
                                     aCGTCTGATAACTTTAACAGACATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCa                                       <  1:1698878/100‑1 (MQ=255)
                                      cGTCTGATATCTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAGGGCAt                                      <  1:5097187/100‑1 (MQ=255)
                                      cGTCTGATAACTTTAACAGGCATTGGGAATGATAACCACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCAt                                      <  1:3790990/100‑1 (MQ=255)
                                      cGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCAt                                      <  1:5647549/100‑1 (MQ=255)
                                       gTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATa                                     <  1:2289886/100‑1 (MQ=255)
                                        tCTTATAACTTTAACAGGCAGTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATaa                                    <  1:4445129/100‑1 (MQ=255)
                                           gATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAGAGGCagatcgg                                 >  1:1375949/1‑94 (MQ=255)
                                           gATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAgat                                 >  1:1496432/1‑98 (MQ=255)
                                           gATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAgat                                 >  1:3288354/1‑98 (MQ=255)
                                           gATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAgat                                 >  1:1317310/1‑98 (MQ=255)
                                           gATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAAGCATAAgat                                 >  1:3822672/1‑98 (MQ=255)
                                             tAACTTTAGCAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAGGCGaa                               <  1:1222861/100‑2 (MQ=255)
                                             tAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTGCAATAAAGGCATAAGCGAt                               <  1:2637460/100‑1 (MQ=255)
                                                 tttAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAgatcggaag                           >  1:2801097/1‑92 (MQ=255)
                                                  ttAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAgatcggaaga                          >  1:995382/1‑91 (MQ=255)
                                                  ttAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGt                          >  1:1632664/1‑100 (MQ=255)
                                                  ttAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAAGCATAAgatcggaaga                          >  1:2664608/1‑91 (MQ=255)
                                                   tAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCAAAAgatcggaagag                         >  1:1573121/1‑90 (MQ=255)
                                                    aaCAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTAc                        >  1:5190831/1‑100 (MQ=255)
                                                     ccAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACGATAAAGGCATAAGCGATAAAGTACa                       >  1:2271784/2‑100 (MQ=255)
                                                       aGGCATTGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGCTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATg                     >  1:4534467/1‑100 (MQ=255)
                                                            ttGGGAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGc                >  1:5268410/1‑100 (MQ=255)
                                                               ggAATGATAACAACGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCAt             <  1:3750740/100‑1 (MQ=255)
                                                                       aacaacGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCATAAATGTTa     >  1:704323/1‑100 (MQ=255)
                                                                        acaacGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCATAAATGTTaa    <  1:5010734/100‑1 (MQ=255)
                                                                         taacGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCATAAATGTTaaa   >  1:4915565/2‑100 (MQ=255)
                                                                          aacGCAAGTAAGTAAAAAAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCATAAATGTTAaag  >  1:5719247/1‑100 (MQ=255)
                                                                                           |                                                                                   
CTGGCATCATTGTAGGCGGTGTATCATTCTTCTTGAAACGTCTGATAACTTTAACAGGCATTGGGAATGATAACAACGCAAGTAAGTAAAATAATGAGCCACCAGGTTTAATAAT‑GATCGTAAGTACAATAAAGGCATAAGCGATAAAGTACATGATTGCCATAAATGTTAAAG  >  NC_007793/1032017‑1032190

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